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(ABSTRACT) 

The  influence  of  freestream  turbulence  representative  of  the  flow  downstream  of  a 
modem  gas  turbine  combustor  and  first  stage  vane  on  turbine  blade  heat  transfer  has  been 
measured  and  analytically  modeled  in  a  linear,  transonic  turbine  cascade.  Measurements 
were  performed  on  a  high  turning,  transonic  turbine  blade.  The  facility  is  capable  of 
heated  flow  with  inlet  total  temperature  of  120°C  and  inlet  total  pressure  of  10  psig.  The 
Reynolds  number  based  on  blade  chord  and  exit  conditions  (5x10^)  and  the  inlet  and  exit 
Mach  numbers  (0.4  and  1 .2,  respectively)  are  representative  of  conditions  in  a  modem 
gas  turbine  engine.  High  intensity,  large  length-scale  freestream  turbulence  was 
generated  using  a  passive  turbulence-generating  grid  to  simulate  the  turbulence  generated 
in  modem  combustors  after  it  has  passed  through  the  first  stage  vane  row.  The  grid 
produced  freestream  turbulence  with  intensity  of  approximately  10-12%  and  an  integral 
length  scale  of  2  cm  near  the  entrance  of  the  cascade  passages,  which  is  believed  to  be 
representative  of  the  core  flow  entering  a  first  stage  gas  turbine  rotor  blade  row.  Mean 
heat  transfer  results  showed  an  increase  in  heat  transfer  coefficient  of  approximately  8% 
on  the  suction  surface  of  the  blade,  with  increases  on  the  pressure  surface  on  the  order  of 
two  times  higher  than  on  the  suction  surface  (approximately  17%).  This  corresponds  to 
increases  in  blade  surface  temperature  of  5-10%,  which  can  significantly  reduce  the  life 
of  a  turbine  blade.  The  heat  transfer  data  were  compared  with  correlations  from 
published  literature  with  good  agreement. 


Time-resolved  surface  heat  transfer  and  passage  velocity  measurements  were 
performed  to  investigate  and  quantify  the  effects  of  the  turbulence  on  heat  transfer  and  to 
correlate  velocity  fluctuations  with  heat  transfer  fluctuations.  The  data  demonstrates 
strong  coherence  in  velocity  and  heat  flux  at  a  frequency  correlating  with  the  most 
energetic  eddies  in  the  turbulence  flow  field  (the  integral  length-scale).  An  analytical 
model  was  developed  to  predict  increases  in  surface  heat  transfer  due  to  freestream 
turbulence  based  on  local  measurements  of  turbulent  velocity  fluctuations  (u'rms)  and 
length-scale  (A*).  The  model  was  shown  to'  predict  measured  increases  in  heat  flux  on 
both  blade  surfaces  in  the  current  data.  The  model  also  successfully  predicted  the 
increases  in  heat  transfer  measured  in  other  work  in  the  literature,  encompassing  different 
geometries  (flat  plate,  cylinder,  turbine  vane  and  turbine  blade)  as  well  as  both  laminar 
and  turbulent  boundary  layers,  but  demonstrated  limitations  in  predicting  early  transition 
and  heat  transfer  in  turbulent  boundary  layers.  Model  analyses  in  the  frequency  domain 
provided  valuable  insight  into  the  scales  of  turbulence  that  are  most  effective  at 
increasing  surface  heat  transfer. 
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Chapter  1 


Introduction 


1.1  Turbine  Blade  Heat  Transfer  -  Problem  Statement 

Gas  turbine  engine  designers  continually  strive  to  improve  engine  performance  by 
increasing  overall  engine  efficiency  and  thrust-to-weight  ratio.  Improvements  in  overall 
engine  efficiency  are  driven  in  large  part  by  increasing  the  combustor  exit  temperature 
and  turbine  inlet  temperature.  Increasing  thrust-to-weight  ratio  necessitates  the  use  of 
lighter  materials,  which  may  operate  near  their  thermal  limits. 

The  implementation  of  blade  cooling  techniques  in  conjunction  with  advances  in 
blade  materials  has  increased  the  envelope  of  turbine  inlet  temperatures.  With  the  blades 
operating  near  their  thermal  limits,  it  is  critical  to  the  engine  designer  to  have  accurate 
predictions  of  hot  mainstream  gas  to  blade  surface  heat  transfer  rates.  The  heat  transfer 
to  the  turbine  blades  is  composed  of  steady  heat  transfer  based  on  mean-flow  conditions, 
as  well  as  unsteady  heat  transfer  due  to  fluctuations  in  the  mean  flow.  The  unsteadiness 
of  the  flowfield  is  a  product  of  wake  passing  and  shock  waves  from  blade  row 
interactions  and  freestream  turbulence  generated  in  the  combustor.  The  effects  of  these 
unsteady  features  of  the  gas  turbine  flowfield  on  blade  heat  transfer  need  to  be  quantified 
to  provide  more  accurate  predictions  of  blade  surface  heat  transfer  and  facilitate  the 

development  of  more  effective  cooling  schemes. 

Analysis  of  the  effects  of  flow  unsteadiness  on  heat  transfer  is  often  performed  by 
decoupling  the  effects  of  the  unsteady  flow  phenomena  and  studying  each  independently. 
Studies  range  from  simple  flat  plate  and  cylinder  stagnation  point  analyses  to  cascade 
testing  and  full  rotating  rigs.  Rotating  rigs  simulate  turbine-operating  conditions  more 
accurately,  however  they  often  make  it  difficult  to  separate  the  effects  of  various 
unsteady  flow  phenomena.  Cascade  testing  provides  the  ability  to  study  heat  transfer 
matching  many  engine  conditions  while  decoupling  the  effects  of  various  unsteady  flow 
phenomena. 
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The  focus  of  the  current  research  work  is  to  investigate  the  effects  of  freestream 
turbulence  on  turbine  blade  heat  transfer  in  a  transonic  cascade  with  high  turning  turbine 
blades.  The  steady  flowfield  matches  engine  non-dimensional  conditions  (Mach  number 
and  Reynolds  number)  and  the  generated  turbulence  field,  intensity  and  length  scale,  is 
consistent  with  recent  reports  from  modem  gas  turbine  engine  combustors  (Van  Fossen 
and  Bunker,  2002)  and  reports  of  turbulence  decay  through  the  first  stage  vane  row 
(Radomsky  and  Thole,  1998,  2000).  The  turbulence  field  simulates  the  core  flow 
turbulence  in  the  turbine,  decoupling  this  turbulence  from  the  turbulence  in  the  wake. 
The  use  of  high  frequency  instrumentation  allows  for  the  analysis  of  both  time-averaged 
and  time-resolved  heat  transfer. 

Following  is  a  review  of  the  literature  of  investigations  of  the  effects  of 
turbulence  on  heat  transfer.  The  review  attempts  to  cover  the  broad  spectrum  of  studies 
which  have  been  performed  on  the  topic  and  also  reviews  turbulence  generation 
techniques  and  time-resolved  measurements  of  heat  flux. 

1.2  Literature  Review 

1.2.1  Turbulent  Heat  Transfer  Studies  and  Correlations 

Many  studies  have  been  performed  in  the  last  25  years  focusing  on  quantifying 
the  effects  of  freestream  turbulence  on  heat  transfer.  Much  of  the  previous  work  has  been 
experimental  analyses  on  simplified  geometries,  namely  the  flat  plate  and  cylinder  in 
cross-flow.  Flat  plate  experiments  present  the  most  basic  geometry  for  experiments,  with 
the  absence  of  curvature  and  pressure  gradient  effects.  Cylinder  experiments  are  of 
particular  importance  as  the  stagnation  point  is  typically  the  region  of  highest  heat 
transfer  on  turbine  vanes  and  blades  except  for  regions  of  transitional  and  turbulent  flow. 
Cylinder  experiments  also  provide  data  on  the  effects  of  curvature  and  pressure  gradient. 
In  more  recent  years,  studies  have  been  performed  in  cascades  of  turbine  blades  and 
vanes.  These  studies  present  data  sets  that  are  generally  more  applicable  to  real  gas 
turbine  engines.  The  following  review  of  previous  work  will  be  presented  by  grouping 
experiments  together  by  classification  of  the  experiments,  rather  than  presenting  a 
chronological  succession  of  studies. 
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1.2. 1. 1  Flat  Plate  Experiments 

Several  research  groups  have  performed  flat  plate  studies  of  turbulent  heat 
transfer.  Bradshaw  (1974)  and  Hancock  and  Bradshaw  (1983)  studied  the  effects  of 
freestream  turbulence  on  turbulent  flat  plate  boundary  layers.  Data  were  collected  with 
and  without  turbulence  and  comparisons  between  the  data  sets  were  made  at  similar 
Reynolds  numbers  based  on  boundary  layer  thickness  (Ree  and  Res).  The  level  of 
freestream  turbulence  in  these  experiments  was  rather  low,  with  reported  intensity  at  7% 
and  below.  A -correlation  was  developed  which  correlates  increases  in  skin  friction  with 
turbulence  intensity,  dissipation  length  scale  (Lu*)  and  boundary  layer  thickness  (5). 
Reynolds  number  was  reported  to  have  no  effect  on  how  freestream  turbulence  affected 
skin  friction  enhancement. 

Blair  (1983)  correlated  skin  friction  and  heat  transfer  increases  due  to  freestream 
turbulence  on  turbulent  flat  plates  and  demonstrated  that  increases  in  freestream 
turbulence  induce  upstream  boundary  layer  transition.  His  work  showed  that  there  was  a 
low  Reynolds  number  effect  that  attenuated  skin  friction  enhancement  and  modified  the 
correlation  of  Hancock  and  Bradshaw  to  include  low  Reynolds  number  effects. 
Measured  heat  transfer  enhancement  was  shown  to  increase  at  a  higher  rate  than  skin 
fnction  with  increasing  turbulence  intensity.  The  Reynolds  analogy  factor  was 
demonstrated  to  increase  linearly  with  increasing  freestream  turbulence  levels.  The  work 
of  Blair  was  later  corrected  based  on  additional  low  Reynolds  number  data  by  Castro 
(1984).  It  was  shown  that  the  low  Reynolds  number  effect  was  less  significant  at  higher 
turbulence  intensities.  Moss  and  Oldfield  (1992)  combined  the  data  sets  of  Blair  and 
Castro  with  their  own  data,  which  included  higher  turbulence  intensity  (up  to  16%)  and 
modified  the  correlation  of  Blair.  Work  by  MacMullin,  et  al.  (1989),  Maciejewski  and 
Moffat  (1989,  1992),  Ames  and  Moffat  (1990)  and  Sahm  and  Moffat  (1992)  also 
extended  the  range  of  the  correlation  developed  by  Blair  and  proposed  a  new  parameter 
(TLR)  which  uses  integral  quantities  rather  than  boundary  layer  thickness,  since  defining 
the  edge  of  the  boundary  layer  at  high  levels  of  freestream  turbulence  can  be  difficult. 

Maciejewski  and  Moffat  (1989)  presented  flat  plate  turbulent  boundary  layer  data 
with  turbulence  intensities  as  high  as  60%  using  a  free  jet.  They  proposed  a  simpler 
method  of  analyzing  high  turbulence  heat  transfer  data.  They  observed  a  linear 
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relationship  between  measured  heat  transfer  coefficient  (h)  and  turbulent  fluctuations  in 
the  freestream  (u’).  Their  correlation  was  independent  of  geometry  or  length  scale,  and 
was  simply  a  function  of  near  wall  turbulence  level.  Although  the  range  of  turbulence 
intensity  is  much  higher  than  what  is  believed  to  be  representative  of  turbine  conditions, 
these  results  are  of  interest.  Later  work  Jjy  Maciejewski  and  Moffat  (1992),  Maciejewski 
and  Anderson  (1996)  and  Denninger  and  Anderson  (1999)  revised  the  correlation  of 
Maciejewski  and  Moffat  (1989)  using  the  value  of  uVms,  max  in  the  inner  wall  region  as  an 
input  parameter  to  improve  the  correlation  and  extend  it’s  applicability  to  lower  levels  of 
freestream  turbulence. 

Thole  and  Bogard  (1995)  performed  a  comprehensive  analysis  of  all  of  these 
worics  and  reported  that  the  correlation  using  integral  quantities  is  more  successful  in 
scaling  heat  transfer  enhancement,  while  Blair’s  correlation  is  more  appropriate  for  skin 
friction  enhancement  predictions.  They  also  concluded  that  the  correlation  of 
Maciejewski  and  Moffat  (1992)  is  good  for  cases  of  very  high  freestream  turbulence 
(Tu>20%).  Relative  to  the  differences  in  measured  heat  transfer  and  skin  friction 
enhancements  in  the  presence  of  high  freestream  turbulence,  they  conclude  that  the 
difference  is  due  to  the  fact  that  at  higher  levels  of  turbulence,  large  scale  eddies 
penetrate  into  the  boundary  layer.  These  large  scale  eddies  are  effective  at  removing 
heat,  but  are  non-stress  producing,  so  at  higher  turbulence  levels,  they  continue  to 
enhance  heat  transfer,  but  their  affect  on  shear  stress  is  attenuated. 

Barrett  and  Hollingsworth  (2003)  studied  the  effects  of  very  small  scale  (ratio  of 
length  scale  to  boundary  layer  momentum  thickness  of  4  to  32),  moderate  intensity  (0.1  - 
8%  Tu)  turbulence  on  skin  friction.  They  developed  a  new  boundary  layer  model  to 
establish  a  skin  friction  correlation  that  does  not  rely  on  knowledge  of  the  momentum 
thickness  Reynolds  number. 

1.2.1.2  Cylinder  Experiments 

Many  of  the  original  works  studying  the  effects  of  freestream  turbulence  on  heat 
transfer  were  cylinder  and  leading  edge  studies.  Since  the  flow  in  the  stagnation  region 
can  be  simulated  by  studying  a  cylinder  in  cross-flow,  these  studies  represent  a  simple 
but  effective  experimental  analysis  which  has  a  valid  application  to  turbine  vane  and 
blade  studies.  As  mentioned,  the  heat  transfer  in  the  leading  edge  region  of  turbine 
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blades  and  vanes  is  typically  the  region  of  highest  heat  transfer  with  the  exception  of 
regions  of  transitional  and  turbulent  flow. 

Smith  and  Keuthe  (1966),  Kestin  and  Wood  (1971)  and  Lowery  and  Vachon 
(1975)  developed  and  modified  correlations  that  correlated  Nu/Re'^  (Frossling  number) 
against  TuRe'^.  The  correlations  were  only  a  function  of  the  turbulence  intensity  (Tu) 
and  cylinder  Reynolds  number  (Re),  length  scale  was  not  included  in  these  studies  as  it 
was  in  later  work  on  cylinders  and  in  most  of  the  work  on  flat  plates.  Mehendale,  et  al 
(1991)  later  modified  the  correlation  for  a  larger  range  of  turbulence  levels  up  to  15%. 

Van  Fossen  and  Simoneau  (1985)  performed  a  streamwise  traverse  of  a  hot-wire 
along  the  stagnation  streamline  for  a  cylinder  in  cross-flow.  From  measurement  of 
velocity  fluctuations  approaching  the  leading  edge  of  the  cylinder  it  was  determined  that 
vorticity  increases  as  it  approaches  the  stagnation  region.  This  has  been  attributed  to 
vortex  stretching  around  the  leading  edge. 

Van  Fossen,  et.  al  (1995)  studied  the  influence  of  turbulence  on  heat  transfer  on 
cylinders  and  cylinders  with  wedge  shaped  afterbodies.  In  these  studies,  it  was 
determined  that  augmentation  of  stagnation  region  heat  transfer  increases  as  integral 
length  scale  decreases.  A  correlation  was  developed  which  modified  those  of  previous 
cylinder  studies  to  include  the  integral  length  scale  in  the  form  of  the  ratio  of  the  length 
scale  to  cylinder  leading  edge  diameter  (Ax/d).  This  type  of  correlation  may  have  direct 
application  to  turbine  blade  leading  edge  heat  transfer  predictions,  using  an  effective 
diameter  based  on  potential  flow  approaching  the  leading  edge.  Since  combustor 
turbulence,  as  will  be  discussed  later,  has  been  shown  to  be  of  fairly  large  scale,  the  lower 
bound  of  the  correlation  presented  by  Van  Fossen  may  not  be  of  particular  interest  in 
direct  application  to  turbine  design. 

Dullenkopf  and  Mayle  (1995)  presented  a  new  correlation  that  introduced  the  idea 
of  an  “effective”  turbulence  level.  The  correlation  was  developed  based  on  stability 
theory  for  a  laminar  boundary  layer.  Fluctuations  at  high  frequencies  (small  scale)  are 
viscously  damped,  while  disturbances  at  low  frequency  (large  scale)  are  seen  as  quasi¬ 
steady  similar  to  bulk  flow  pulsations,  which  in  the  mean  do  not  affect  heat  transfer. 
Only  a  small  frequency  band  around  a  dominant  frequency  is  effective  at  increasing  heat 
transfer.  A  new  turbulence  parameter  (Tux)  was  introduced  which  accounts  for  effective 
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turbulence  level  and  length  scale.  They  introduced  a  Stanton  number  that  was  linearly 
proportional  to  the  effective  turbulence  level  (Uefr),  similar  to  the  results  proposed  by 
Maciejewski  and  Moffat  (1992). 

Van  Fossen  and  Bunker  (2002)  measured  stagnation  region  augmentation 
downstream  of  a  section  of  a  dual-annular  combustor.  The  turbulence  intensity  and 
length  scale  downstream  of  the  combustor  section  were  mapped  out  using  a  two 
component  hot-wire.  The  turbulence  field  was  shown  to  have  a  maximum  intensity  of 
approximately  *32%  with  an  integral  length  scale  of  1.29  cm.  Heat  transfer  augmentation 
was  predicted  well  by  correlations.  The  reported  turbulence  levels  exiting  the  combustor 
were  somewhat  higher  than  measured  by  other  researchers  (Goldstein,  et  al.,  1983, 
Kuotmos  and  McGuirk,  1989  and  Ames  1994,  1995  and  1997)  where  levels  of  20-25% 
were  measured. 

1.2.1. 3  Turbine  Vane  and  Blade  Experiments 

Several  research  groups  have  performed  experiments  on  turbine  vane  geometries 
in  low  speed  cascade  facilities.  Cascade  experiments  present  a  more  realistic  simulation 
of  the  effects  of  turbulence  on  turbine  heat  transfer  and  results  are  generally  more 
applicable  to  turbine  conditions  than  flat  plate  data. 

Ames  (1995,  1997)  performed  measurements  of  the  effects  of  large-scale 
freestream  turbulence  in  a  low  speed  turbine  vane  cascade.  In  these  studies,  three 
different  turbulence  fields  were  generated  with  varying  intensity  and  scale.  Freestream 
turbulence  was  shown  to  significantly  augment  heat  transfer  in  the  laminar  regions  of  the 
blade  (stagnation  point,  pressure  surface  and  favorable  pressure  gradient  region  of  the 
suction  surface).  Heat  transfer  augmentation  was  reported  to  be  highest  in  the  stagnation 
region  and  on  the  pressure  surface.  Stagnation  region  increases  agreed  well  with  cylinder 
stagnation  point  correlations.  The  average  augmentation  of  heat  transfer  on  the  pressure 
surface  was  found  to  scale  well  with  the  inlet  turbulence  intensity  and  the  length  scale  to 
the  -1/3  power  (TuLu'^^).  Comparisons  of  Stanton  number  increases  on  both  the 
pressure  and  suction  surfaces  of  the  vane  demonstrated  that  turbulence  has  a  much 
greater  effect  on  the  pressure  surface  of  the  vane  as  compared  to  the  suction  surface. 
Increased  turbulence  slightly  shifted  suetion  surface  transition  towards  the  leading  edge. 
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In  the  turbulent  region  of  the  suction  surface,  increased  turbulence  had  a  smaller  effect  on 
heat  transfer. 

Radomsky  and  Thole  (1998,  2000,  2002)  documented  surface  heat  transfer  and 
the  turbulent  flowfield  around  a  scaled  up  stator  vane.  Several  different  turbulence  levels 
as  high  as  25%  were  generated.  Observations  of  heat  transfer  augmentation  were 
consistent  with  results  reported  by  Ames  (1995),  with  heat  transfer  augmentation  being 
greatest  on  the  pressure  surface  and  in  the  stagnation  region  and  suction  surface  transition 
shifting  toward  the  leading  edge  at  elevated  freestream  turbulence  levels.  In  Radomsky 
and  Thole  (1998),  detailed  measurements  of  unsteady  velocity  upstream  of  the  vane 
passage  and  along  a  streamline  in  the  center  of  the  vane  passage  were  used  to  document 
the  levels  of  turbulence  kinetic  energy  and  all  three  unsteady  components  of  velocity 
through  the  vane  passage.  The  turbulent  kinetic  energy  was  shown  to  increase  dirough 
most  of  the  blade  passage.  The  streamwise  velocity  fluctuations  (u*)  were  shovm  to 
decrease  through  the  blade  passage,  while  cross-stream  (v')  and  lateral  (w^  fluctuations 
increased.  Detailed  measurements  using  a  3  component  LDV  system  were  reported  in 
Radomsky  and  Thole  (2000)  and  showed  that  turbulence  remains  quite  high  throughout 
the  passage.  Length  scale  was  shown  to  decrease  through  the  blade  passage,  although  not 
considerably.  The  effect  of  turbulent  length  scale  was  reported  to  have  the  largest  effect 
on  the  pressure  surface  of  the  vane,  with  lower  Stanton  numbers  at  the  larger  scales. 

Mehendale,  et  al.  (1994)  and  Wang,  et  al.  (1998)  measured  the  effects  of  high 
intensity,  large  scale  turbulence  on  turbine  blade  heat  and  mass  transfer  in  a  low  speed 
cascade  at  varying  turbulence  levels.  In  general,  results  of  their  tests  agreed  with  those 
from  turbine  vane  testing  in  terms  of  locations  of  highest  heat  and  mass  transfer  and  early 
boundary  layer  transition  on  the  suction  surface.  Leading  edge  results  agreed  well  with 
published  correlations  from  Kestin  and  Wood  (1971)  and  Dullenkopf  and  Mayle  (1995). 

Holmberg  (1996)  studied  the  effects  of  small-scale,  relatively  low  intensity 
turbulence  on  heat  transfer  in  high-turning  turbine  blades  in  transonic  flow  and  measured 
time-resolved  heat  flux  and  velocity.  Coherence  between  freestream  velocity  and  surface 
heat  flux  was  found  to  be  useful  in  determining  the  scale  and  frequency  range  of 
turbulent  structures  interacting  with  surface  heat  flux.  This  study  is  similar  to  the  current 
work  except  with  lower  intensity,  smaller  scale  turbulence. 
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1.2.2  Turbulence  Generation  Methods 

Experimental  measurements  of  turbulent  heat  transfer  cannot,  of  course,  be 
performed  without  the  means  of  generating  the  desired  turbulent  flowfield.  One 
consequence  of  the  many  turbulence  studies  discussed  in  the  previous  section  is  a  vast 
array  of  turbulence  generation  techniques.  A  fairly  comprehensive  discussion  is 
presented  in  Nix,  et  al.  (2002)  in  Appendix  A.  The  paper  breaks  down  the  generation 
techniques  into  two  categories,  mock  combustor/jets  in  crossflow  and  grid  or  bar 
generators.  A  summaiy  of  the  grid  generation  techniques  is  presented  in  Table  1.1. 

From  the  summaiy,  it  can  be  seen  that  active  blowing  (either  through  jets  in  cross 
flow  or  through  active  grids  or  bar  arrays)  and  passive  bar  arrays  with  large  flow 
blockage  are  effective  at  producing  high  levels  of  turbulence.  The  current  work 
attempted  to  use  active  blowing  in  the  transonic  facility,  as  discussed  in  detail  later  in 
section  2.4,  but  the  momentum  flux  ratio  (I=pju//pooU«>^)  required  was  not  achievable. 
The  grid  used  in  this  work  was  a  large  blockage  grid  based  on  the  results  reported  by 
Polanka  and  Bogard  (1999). 

1.2.3  Time-Resolved  Heat  Transfer 

The  analyses  presented  later  in  this  work  rely  on  the  measurement  of  time- 
resolved  heat  transfer  and  velocity.  Limited  time-resolved  heat  flux  experiments  exist  in 
the  open  literature.  Heat  flux  microsensors  have  been  used  to  measure  time-resolved  heat 
flux  at  Virginia  Tech  by  Johnson  (1995),  Holmberg  (1996),  Nix,  et  al.  (1997a,  b),  and 
Popp,  et  al.  (2000).  Time-resolved  measurements  were  performed  to  measure  unsteady 
fluctuations  and  unsteady  shock  passing  events. 

Thin  film  heat  flux  gages  have  been  used  to  measure  unsteady  heat  transfer  in 
rotating  rigs  at  CALSPAN  by  Dunn  (1986,  1989),  at  MIT  by  Guenette  (1989)  and  at 
Oxford  by  Ashworth  (1989),  Garside  et  al.  (1994)  ,  by  Moss  and  Oldfield  (1992,  1996) 
and  by  Anthony,  et  al  (1999)  at  Oxford  and  Syracuse. 

Of  the  time-resolved  heat  flux  data  sets  published  in  the  open  literature,  only  three 
included  simultaneous  measurements  of  velocity.  Moss  and  Oldfield  (1992, 1996)  and 
the  previous  work  by  Holmberg  (1996)  and  Holmberg  and  Pestian  (1996)  are  the 
exceptions. 
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Table  1.1:  Review  of  Turbulence  Generation  Techniques 


Investigator 


Turbulence  Length  Scale 
Level  cm  (A^c) 


Method 


Flow 

Condition 


Mehendale  and  Han 
(1992) 

13% 

3 

Blown  Grid 

Low  speed 

Maciejewski  and 
Moffat  (1992) 

20-60% 

- 

Turbulent  free  jet 

- 

Sahm  and  Moffot 
(1992) 

4-9% 

13-21% 

- 

■ 

Grid,  blown  grid  and 
jets  in  crossflow 

- 

Thole,  et  al.  (1994) 

20% 

2-3 

- 

Jets  in  crossflow 

Low  speed 

Holmberg  (1996) 

5% 

moil 

<0.1 

Grid 

Transonic 

Schauer  and 

Pestian  (1996) 

17% 

8 

- 

Jets  in  crossflow 

Low  speed 

Ames  (1997) 

13% 

1.6 

0.11 

Mock  combustor 

Low  speed 

Boyle,  et  al.  (1998) 

10% 

1.5 

■ 

Blown  and  unblown 
grid 

Low  speed 

Wang,  et  al.  (1998) 

18% 

8 

0.43 

Mock  combustor 

Low  speed 

Radomsky  and 

Thole  (2000) 

20% 

5-6 

- 

Blown  bar  array 

Low  speed 

Polanka  and 

Bogard  (1 999) 

20% 

5 

- 

Passive  bar  array 
(large  blockage) 

Low  speed 

Giel,etal.(2000) 

10% 

2.2 

0.17 

Active  grid 

Transonic 

Barringer,  et  al. 
(2002) 

18% 

■ 

0.11 

Combustor 
simulator  with  large 
jets 

Low  speed 

Moss  and  Oldfield  (1992, 1996)  traversed  a  hot-wire  probe  above  a  flat  plate  with 
thin  film  gages  on  the  surface.  They  performed  time-domain  analyses  of  the  time- 
resolved  signals,  comparing  raw  signals  and  cross  correlations  between  the  velocity  and 
heat  flux  signal  and  between  heat  flux  signals  of  different  gages.  The  raw  heat  flux  time 
signal  showed  large  and  small  fluctuations.  The  large  fluctuations  were  attributed  to 
freestream  turbulent  eddies  acting  on  the  surface,  while  the  small  fluctuations  were  due  to 
turbulence  within  the  boundary  layer.  Through  correlations  of  velocity  and  heat  flux  and 
spatial  correlations  between  heat  flux  gages  they  concluded  that  the  turbulent  eddy 
structure  of  the  boundary  layer  is  dominated  by  the  freestream  turbulence  and  eddies 
within  the  boundary  layer  convected  at  the  freestream  velocity,  U,  rather  than  at  ~0.8U 
characteristic  of  boundaiy  layers. 
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Holmberg  (1996),  discussed  earlier,  performed  measurements  of  unsteady  surface 
heat  flux  with  a  hot-wire  collocated  with  the  heat  flux  gage.  Holmberg  and  Pestian 
(1996)  measured  u'  and  v'  velocity  components  and  fluctuating  surface  temperature  (f) 
and  heat  flux  (q*)  in  a  low  speed  wall-jet  facility.  Surface  heat  flux  was  shown  to  be 
dominated  by  the  u'  fluctuating  velocity- field  primarily  by  large  freestream  eddies  which 
reach  deep  into  the  boundary  layer.  These  results  are  consistent  with  the  observations  of 
Moss  and  Oldfield  (1992, 1996). 

1.3  Overview  OF  Current  Work 

The  current  experimental  work  at  Virginia  Tech  studying  the  effects  of  freestream 
turbulence  extends  work  performed  by  Holmberg  (1996).  The  blade  geometry  is  the 
same  as  that  studied  by  Holmberg,  but  has  been  scaled  up  three  times  larger  than  the 
previous  work.  The  advantage  of  the  increase  in  blade  size  is  that  it  decreases  the  relative 
size  of  the  heat  flux  sensor  footprint  compared  to  the  size  of  the  blade  surface.  This 
decreases  spatial  averaging  due  to  better  spatial  resolution  compared  to  previous 
experiments.  Another  advantage  to  the  increase  in  scale  is  that  the  pitch  to  chord  ratio 
(p/c)  has  not  changed,  so  the  ratio  of  the  diameter  of  the  hot-wire  probes  used  to  make 
passage  measurements  to  the  passage  width  is  much  smaller.  This  minimizes  any  flow 
disturbance  which  may  be  caused  by  the  hot-wire  probe  support.  Probe  disturbances 
would  be  seen  in  the  frequency  spectrum  of  velocity  and  heat  flux. 

The  turbulence  field  in  the  current  work  is  also  considerably  different  from  that 
investigated  in  the  previous  study.  The  turbulence  field  of  Holmberg  was  not  what  is 
believed  to  be  representative  of  engine  conditions.  The  turbulence  intensity  approaching 
the  turbine  blade  cascade  in  the  current  experiments  is  almost  three  times  greater,  and  the 
integral  scale  to  blade  chord  ratio  (Ax/c)  is  between  0.13  and  0.16,  as  compared  to  a 
maximum  value  of  only  0.06  in  the  work  of  Holmberg.  The  desired  Ax/c  value  was 
estimated  based  on  combustor  turbulence  measurements  by  Van  Fossen  and  Bunker 
(2002).  The  value  used  in  this  work  closely  matches  that  reported  by  Van  Fossen. 

The  work  of  Holmberg  (1996)  focused  on  a  frequency  domain  analysis  of  surface 
heat  transfer  and  freestream  turbulence  interactions.  This  work  employs  many  of  the 
tools  used  by  the  previous  work,  but  focuses  more  on  quantifying  the  effects  of 
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turbulence  on  heat  transfer  and  developing  an  analytical  model  for  explaining  and 
predicting  turbulent  heat  transfer. 

This  work  reports  the  measurements  of  baseline  (low  turbulence)  time-averaged 
heat  transfer  and  the  increases  in  heat  transfer  due  to  the  generated  freestream  turbulence 
field  on  both  the  suction  and  pressure  .surfaces  of  the  turbine  blades  at  three  locations. 
The  results  of  time-averaged  heat  transfer  measurements  are  presented  in  Section  3.2.3. 

Time-resolved  measurements  of  heat  flux  and  velocity  are  also  reported  and 
analyzed.  Analytical  modeling  of  a  mechanism  by  which  freestream  turbulence  affects 
surface  heat  transfer  is  documented  in  Section  3.4.  The  analytical  model  development 
and  application  is  detailed  in  both  a  “single-scale”  and  “multiple-scale”  framework.  The 
“single-scale”  model  (Section  3.4.1)  utilizes  time-resolved  velocity  data  in  the  time- 
domain  to  determine  model  input  parameters.  The  model  is  directly  applied  to  the  time- 
averaged  heat  transfer  results  measured  in  the  current  work,  as  well  as  to  results  reported 
by  other  research  groups.  The  “multiple-scale”  model  looks  at  both  time-resolved 
velocity  and  heat  flux  data  in  the  frequency  domain  (Section  3.4.4).  The  purpose  of  this 
application  of  the  analytical  model  was  to  attempt  to  determine  which  frequencies  (eddy 
scales)  had  the  largest  impact  on  surface  heat  transfer. 
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Chapter  2 


Experimental  Setup 

2.1  Transonic  Cascade  Wind  Tunnel 

The  experimental  facility  in  which  the  majority  of  the  research  for  this  work  was 
performed  was  the  transonic  cascade  wind  tunnel  at  Virginia  Tech.  The  wind  tunnel  is  a 
blowdown  facility,  which  controls  tunnel  inlet  total  pressure  to  a  desired  value  through  a 
feedback  control  loop.  The  tunnel  is  capable  of  run  times  of  up  to  30  seconds  with 
constant  inlet  pressure  and  heated  flow.  A  schematic  of  the  facility  is  show  in  Figure  2.1 . 

High-pressure  air  at  approximately  1240  kPa  (180  psig)  is  supplied  by  a  four 
stage,  reciprocating  compressor  to  two  large  storage  tanks.  Moisture  and  contaminants 
are  removed  from  the  compressed  air  by  an  activated  alumina-drying  unit.  Between  the 
facility  and  the  supply  pipe  from  the  storage  tanks  are  two  automatic  valves  in  series. 
The  first  valve  is  an  open-close  valve,  which  isolates  the  high-pressure  air  from  the 
tunnel  until  the  valve  is  opened.  Once  this  valve  is  opened,  air  flows  through  the 
downstream  valve,  which  is  initially  opened  slightly  awaiting  a  feedback  pressure.  The 
downstream  valve  will  actuate  based  on  feedback  pressure  measured  by  a  total  pressure 
probe  at  the  test  section  inlet  to  maintain  a  constant  inlet  pressure.  A  PC  based  control 
program  using  both  a  feedback  and  feed  forward  algorithm  allows  the  desired  objective 
total  pressure  to  be  set  and  controls  the  actuation  of  the  valve.  Figure  2.2  shows  a  sample 
time  history  of  tunnel  inlet  total  pressure.  To  achieve  the  design  Mach  number  and 
Reynolds  number  for  the  cascade,  an  inlet  objective  total  pressure  was  set  to  a  value  of  69 
kPa  (10  psig). 

Prior  to  entering  the  test  section,  the  tunnel  flow  passes  through  two  bundles  of 
copper  tubing,  which  serve  as  a  passive  heat  exchanger.  The  heat  exchanger  tube  bundles 
are  charged  by  a  heat  exchange  loop  built  into  the  wind  tunnel.  Two  valves,  shown  in 
Figure  2.1,  are  closed  to  isolate  the  tube  bundles  from  the  mainstream  flow  path.  Air  is 
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heated  by  two  36  kW  electric  heaters  and  circulated  through  the  heating  loop  by  an  axial 
fan  at  approximately  2  m/s,  heating  the  copper  tubes.  The  valves  are  closed  prior  to 
wind  tunnel  startup,  restoring  the  path  for  mainstream  airflow  through  the  tunnel.  The 
tunnel  air  passes  over  the  copper  tube  bundles,  which  transfer  heat  to  the  mainstream 
flow.  The  air  leaving  the  heat  exchanger  and  entering  the  blade  cascade  test  section  can 
be  heated  as  high  as  120°C  (248°F).  Thermocouples  and  a  total  temperature  probe  are 
located  within  the  heating  loop  and  tube  bundle  to  monitor  the  tube  temperature  during 
heating  and  the  total  temperature  of  the  flow  during  a  tunnel  run,  respectively.  Due  to  the 
passive  nature  of  the  heat  exchanger,  the  total  temperature  of  the  air  decreases  during  a 
tunnel  run  as  the  tubes  lose  energy  and  transfer  heat  to  the  flow.  As  will  be  discussed 
later,  heat  transfer  data  reduction  techniques  take  advantage  of  the  transient  nature  of  the 
flow  temperature. 

For  experiments  with  high  freestream  turbulence,  a  passive  turbulence  grid  is 
placed  between  the  heat  exchanger  and  the  test  section  inlet  (shown  in  Figure  2.1).  The 
turbulence  grid  is  briefly  discussed  in  Section  2.4  and  details  of  the  design  and  testing  of 
the  grid  are  presented  in  Appendix  A. 


Figure  2.1:  Virginia  Tech  transonic  blowdown  wind  tunnel  facility 
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Time  (s) 


Figure  2.2:  Typical  tunnel  inlet  total  pressure  plot 

2.2  Cascade  Test  Section  and  Blade  Instrumentation 

The  turbine  blade  cascade  test  section  used  in  the  facility  consists  of  4  full  and 
two  half  blades,  creating  a  total  of  five  flow  passages.  Figure  2.3  shows  an  illustration  of 
the  turbine  cascade.  The  turbine  blades  are  held  in  place  by  two  Plexiglas  end  walls, 
allowing  for  shadowgraph  or  Schleiren  flow  visualization.  The  blade  design  is  a  high 
turning,  generic  GE  first  stage  rotor  geometry.  The  predicted  Mach  number  profile  for 
the  blade  is  shown  in  Figure  2.4.  The  flow  exits  the  test  section  and  is  exhausted  to 
atmospheric  conditions  through  a  muffler.  The  non-dimensional  conditions  of  the  flow  in 
the  cascade  are  representative  of  operating  conditions  of  modem  gas  turbine  engines. 
The  non-dimensional  conditions  and  the  cascade  conditions  and  geometry  are  shown  in 
Table  2.1  below. 
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Table  2.1:  Non-dimensional  conditions  and  cascade  conditions  and  geometry 


Non-dimensional  Conditions 

Engine 

Experiment 

Mach  No.  Exit 

Transonic 

1.2 

Reynolds  No.  Inlet 

5x10* 

1x10* 

Temperature  Ratio  (Tg/Tj,) 

1.4 

1.4 

1  Cascade  Conditions  and  Geometry  | 

Inlet  Mach 

0.4 

Exit  Mach 

1.2 

Inlet  total  pressure 

164kPa(10psig) 

Inlet  total  temperature 

lOOT 

Reynolds  No.  Exit  (ROex) 

5x10* 

Blade  aerodynamics  chord  (cj 

13.6  cm 

Suction  surface  chord 

24.8  cm 

Pressure  surface  chord  (Cpj) 

14.9  cm 

Blade  span  (sp) 

15.3  cm 

Blade  pitch  (p) 

12.2  cm 

sp/c 

1.12 

p/c 

0.9 

2-0 

Contraction 


Figure  2  Cascade  test  section  illustration  showing  instrumented  blade  location 
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Figure  2.4:  Predicted  Mach  number  distribution  for  the  turbine  blades 

One  blade  near  the  center  of  the  cascade  is  instrumented  on  both  the  pressure  and 
suction  surfaces  at  roughly  mid-span  to  avoid  endwall  boundary  layer  effects.  The 
sensors  were  staggered  in  the  spanwise  direction  to  prevent  small  flow  disturbances  in  the 
boundaiy  layer  caused  by  the  gage  mounting  from  effecting  downstream  gages. 
Measurements  were  made  at  locations  on  both  the  pressure  and  suction  surfaces  of  the 
blade  to  compare  the  effects  of  turbulence  on  heat  transfer  under  greatly  different  flow 
conditions  (differences  in  mean  velocity,  pressure  gradient  and  curvature  and  velocity 
fluctuations).  For  the  results  presented  in  this  work,  one  location  on  the  suction  surface 
(SSI)  was  used  and  two  locations  on  the  pressure  surface  (PS  1  and  PS2).  Each 
measurement  location  consists  of  an  array  of  three  sensors,  as  shown  in  Figure  2.5.  Each 
triplet  of  sensors  includes  a  Vatell  HFM-7  heat  flux  microsensor,  a  Kulite  XCQ-062-50a 
pressure  transducer  and  a  K-type  thermocouple.  The  heat  flux  sensors  and  pressure 
transducers  are  discussed  in  Section  2.3. 
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Pressure  Surface 


Location  SSI . 


,  KuRte  Pressure  Transducer 
~  Heat  Flux  Mtaosensor 
Surface  Thermocouple 


Suction  Surface 


Figure  2.5:  Turbine  blade  showing  measurement  locations  and  instrumentation 

details 


The  location  of  the  gages  on  the  two  surfaces  was  selected  such  that  the  gages  were 
in  regions  of  high  heat  transfer  and  for  ease  of  mounting.  Since  the  gages  are  insert  type, 
only  regions  of  the  blade  where  there  is  sufficient  material  for  machining  and  mounting 
of  the  sensors  could  be  instrumented.  The  instrumented  blade  was  machined  in  two 
pieces  such  the  blade  could  be  split  in  half  to  mount  the  gages  so  that  they  are  flush  with 
the  blade  surface  (Bubb,  1999).  The  suction  surface  gage  (SSI)  is  located  at  a  distance  of 
x/Css=0.20  from  the  stagnation  point  of  the  turbine  blade.  The  two  gages  on  the  pressure 
surface  of  the  blade  (PSl  and  PS2)  are  located  at  x/Cps=0.16  and  x/cps=0.40,  respectively. 
The  chord  length  of  the  blade  on  each  surface  is  denoted  as  Css  for  the  suction  surface  and 
Cps  for  the  pressure  surface.  Based  on  aerodynamic  chord  length  (ca),  SSI  is  located  at 
x/ca=0.36,  PSl  is  at  x/ca=0.20,  and  PS2  is  at  x/Ca=0.43.  The  location  of  the  gages  relative 
to  the  Mach  number  distribution  is  shown  in  Figure  2.4.  It  should  be  noted  that  the  two 
locations  on  the  pressure  surface  are  located  in  regions  of  relatively  different  curvature 
and  pressure  gradient.  Table  2.2  summarizes  the  conditions  at  the  gage  locations. 
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Table  2.2:  Blade  Measurement  Locations 


Location 

x/Css  or  ps 

X/Ca 

Approx. 

Mach  Number 

Pressure 

Gradient  (+/-) 

SSI 

0.20 

0.36 

0.55 

+  Favorable 

PSl 

0.18 

0.20 

0.22 

-  Adverse 

PS2 

0.40 

0.43 

0.25 

+  Favorable 

2.3  Instrumentation  and  Data  Acquisition  Systems 

In  order  to  perform  measurements  for  the  analysis  of  turbulent  heat  transfer, 
instrumentation  measuring  both  steady  and  unsteady  mainstream  flow  and  blade  surface 
conditions  was  employed.  Steady  measurements  were  used  to  determine  flow  conditions 
within  the  turbine  cascade  such  as  inlet  Mach  number  and  Reynolds  number,  as  well  as 
steady  heat  transfer  coefficients  and  static  pressures  on  the  blade  surface.  Unsteady 
measurements  were  performed  to  measure  time-resolved  velocity  and  heat  flux  data. 
Following  is  a  detailed  description  of  the  instrumentation  used  in  all  experiments. 

2.3.1  Vatell  HFM-7  Heat  Flux  Sensors 

Measurements  of  heat  flux  on  the  blade  surface  were  performed  using  Vatell 
HFM-7  heat  flux  sensors.  The  heat  flux  microsensor  (HFM)  technology  is 
documented  in  published  papers  by  Holmberg  and  Diller  (1995),  Hager,  et.al.  (1991),  and 
Baker  and  Diller  (1 993).  Calibration  techniques  for  the  heat  flux  sensor  are  discussed  in 
detail  in  Smith,  et  al.  (1999).  The  heat  flux  microsensor  consists  of  two  sensors,  one  for 
measurement  of  heat  flux  (HFS)  and  a  resistance  temperature  sensor  (RTS).  The  sensor 
is  an  insert  type  gage  that  is  designed  to  be  inserted  into  a  machined  slot  such  that  the 
sensor  face  is  flush  with  the  measurement  surface.  The  sensor  face  is  approximately  6.35 
mm  (0.25  in.)  in  diameter,  which  includes  both  the  heat  flux  sensor  footprint  (4.5  nun) 
and  the  resistance  temperature  sensor.  The  gages  are  sufficiently  small  compared  to  the 
blade  chord  and  span  to  provide  reasonably  good  spatial  resolution  and  the  ratio  of  the 
turbulence  integral  length  scale  to  the  sensor  diameter  is  approximately  Ax/Ds=5.5.  Using 
the  sensor  diameter  and  a  convective  velocity  of  approximately  100  m/s,  the  frequency 
resolution  of  the  gages  is  estimated  as  approximately  20  kHz.  The  heat  flux  sensors 
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attenuate  energy  above  this  frequency.  As  illustrated  in  the  results  presented  in  the  next 
chapter,  this  is  well  above  the  frequency  range  of  interest,  since  there  is  no  coherence 
between  velocity  and  heat  flux  above  5  kHz.  The  scales  of  interest  discussed  later  are  all 
on  the  order  of  the  gage  diameter  or  larger.  An  illustration  of  the  HFM-7  is  shown  in 
Figure  2.6a.  The  heat  flux  sensors  are  made  on  a  substrate  of  aluminum  nitride,  which 
has  thermal  properties  similar  to  those  of  aluminum,  minimizing  error  in  measurement  of 
heat  flux  from  a  mismatch  of  materials. 


Figure  2.6a.  Vatell  HFM-7  Heat  Flux  Microsensor 
Source:  Vatell  Corporation  website 

The  heat  flux  sensor  (HFS)  gives  an  output  voltage  that  is  linearly  proportional  to 
heat  flux.  The  sensitivity  of  the  HFS  is  approximately  150-200  pVAVW,  depending 
on  the  gage.  The  HFS  sensitivity  includes  a  small  correction  based  on  gage  temperature. 
The  resistance  temperature  sensor  (RTS)  generates  a  voltage  output  with  temperature 
changes  by  applying  a  0.1  mA  excitation  to  the  resistance  element  on  the  sensor.  The 
sensor  resistance  (Ro)  is  measured  when  the  gage  is  at  thermal  equilibrium  with  its 
surroundings  and  the  voltage  ouq>ut  is  zeroed  at  this  temperature.  Measured  voltages  are 
converted  to  a  resistance  using  Ohm’s  Law  (V=IR)  and  the  known  excitation  current  and 
amplifier  gain  and  are  then  added  to  the  zeroed  resistance  (Ro).  The  sensor  temperature  is 
then  determined  through  a  calibrated  linear  resistance/temperature  curve  for  each  gage,  as 
shovm  in  Equations  2.1  and  2.2. 

R=-r^ — (2.1) 

^  RTS  '^RTS 
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T  =  c-R  +  d 


(2.2) 


where  Vrts  is  the  measured  voltage  output  of  the  RTS  sensor,  Irts  is  the  applied 
excitation  current  (0.1  mA),  Grts  is  the  RTS  amplifier  gain  and  c  and  d  are  the 
coefficients  of  the  linear  resistance/temperature  curve  for  the  gage.  Heat  flux  from  the 
measured  HFS  voltage  signal  was  determined  by: 


"  _YhfsJ_£hfs_ 
^  g-T  +  h 


(2.3) 


where  Vhfs  is  the  measured  voltage  from  the  HFS  sensor,  Ghfs  is  the  HFS  amplifier  gain 
and  g  and  h  are  coefficients  of  the  sensitivity  of  the  sensor,  accounting  for  gage 
temperature  (T).  The  sensors  are  each  connected  to  a  Vatell  AMP-6  amplifier,  which 
provide  the  HFS  and  RTS  gains  and  RTS  excitation  current.  The  gain  for  both  the  HFS 
and  RTS  was  typically  set  to  100,  providing  an  easily  measured  signal  with  low  signal-to- 
noise  ratio.  The  output  from  the  amplifier  is  connected  to  the  data  acquisition  system 
through  a  BNC  connection. 

The  heat  flux  sensors  were  used  to  perform  measurements  of  both  steady  heat 
transfer  coefficient  and  for  time-resolved  measurements  of  heat  flux.  For  steady  heat 
transfer  measurements,  output  voltages  from  the  sensor  were  recorded  at  100  Hz 
sampling  frequency  and  were  sampled  over  the  course  of  a  30  second  tunnel  run.  For 
measurements  of  time-resolved  heat  transfer,  data  were  recorded  at  100  kHz  for  a  short 
duration  of  1.3  s  and  filtered  at  40  kHz.  The  sampling  period  gives  a  sufficient  sample 
size  for  statistical  analysis  while  limiting  the  output  file  sizes.  Popp  (1999)  investigated 
the  dynamic  behavior  of  the  HFM-7  gage  in  detail.  It  was  demonstrated  that  the  gage 
responds  like  a  first  order  system  with  a  time  constant  of  17  ps,  or  a  frequency  response 
of  approximately  60  kHz.  This  provides  a  wide  enough  frequency  band  for  sampling 
time-resolved  heat  flux  data  for  correlation  with  high  frequency  velocity  data. 

2.3.2  Kulite  Pressure  Transducers 

The  pressure  transducers  used  in  this  experiment  were  Kulite  XCQ-062-50  high- 
response  miniature  transducers  (Figure  2.6b).  The  transducers  were  used  to  measure 
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surface  pressure  to  determine  the  Mach  number  near  the  blade  surface.  The  transducers 
have  a  diameter  of  1.7  mm  (0.067  in)  and  are  approximately  1  cm  (0.4  in.)  in  length. 
Each  transducer  is  equipped  with  a  B-screen,  which  protects  the  internal  diaphragm  from 
being  damaged  by  debris  moving  through  the  test  section.  The  sensor  is  also  an  insert 
type  gage  that  is  flush  mounted  with  the^jneasurement  surface. 

The  Kulite  sensors  are  located  at  the  same  chordwise  locations  as  the  heat  flux 
sensors,  but  are  slightly  offset  from  the  HFM  gages  in  the  spanwise  direction,  as  shown 
in  Figure  2.5. '  The  lead  wires  for  all  of  the  Kulites  are  connected  to  a  Measurements 
Group  2310  signal  conditioner/amplifier  via  a  bridge  board.  The  2310  amplifler  provides 
the  excitation  voltage  to  the  transducer  and  applies  a  gain  to  the  output  signal.  The 
voltage  output  from  the  transducer  is  zeroed  at  ambient  conditions  such  the  transducer 
gives  an  output  proportional  to  gage  pressure  (psig).  The  amplifler  also  has  built  in  filter 
settings.  For  the  experiments  here,  the  gages  were  not  used  to  perform  time-resolved 
pressure  measurements;  therefore  the  frequency  bandwidth  of  the  transducers  was  not  of 
importance  and  the  amplifiers  were  set  to  “wide-band”.  The  transducers  are  capable  of 
high-frequency  measurements  of  static  pressure  and  have  a  frequency  range  of  up  to  25 
kHz. 

The  output  voltages  from  all  Kulites  were  reduced  to  gage  pressure  (psig)  using  a 
known  sensitivity  of  each  transducer.  The  sensitivity  of  the  Kulites  used  for  these 
experiments  is  typically  between  2.0  and  4.0  mV/psi.  The  general  formula  for  reducing 
the  Kulite  data  is  as  follows: 

P(psig)=-fe^  (2,4) 

p 

where  Vp  is  the  amplified  output  voltage  measured  from  the  transducer,  Sp  is  the 
sensitivity  of  the  transducer,  and  Gp  is  the  gain  setting  on  the  amplifier.  The  gain  for  all 
pressure  transducers  used  during  experimentation  was  set  to  100. 
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Figure  2.6b:  Unmounted  Kulite  pressure  transducer 


2.3.3  Hot-wire  Probes  and  Anemometer 
2.3.3.1  Probe  Design 

Measurements  of  time-resolved  velocity  were  performed  with  single  sensor  hot¬ 
wire  probes  (straight  wire).  The  probes  were  used  to  make  measurements  in  the 
mainstream  flow  upstream  of  the  cascade  and  in  the  turbine  blade  passages  near  the  blade 
surface.  Mainstream  measurements  were  performed  to  characterize  the  turbulence  field 
downstream  of  the  grid  (turbulence  intensity  and  length  scale)  and  to  measure  the  decay 
of  turbulence  intensity  and  dispersion  of  length  scale  as  the  flow  progresses  through  the 
test  section,  towards  the  cascade  of  blades.  Passage  measurements  were  performed  to 
measure  time-resolved  velocity  near  the  surface  of  the  instrumented  blade. 

For  mainstream  flow  measurements,  a  Dantec  Type  55P71  hot-wire  probe  and 
several  custom  made  hot-wire  probes  were  used.  The  wire  sensing  area  on  all  probes  is 
an  annealed  tungsten  wire  with  a  diameter  of  5  pm  and  length  between  1-1.5  mm.  The 
probes  are  shown  in  Figure  2.7a  and  Figure  2.7b.  Figure  2.7a  details  the  dimensions  of 
the  custom  made  hotwires,  while  Figure  2.7b  shows  a  photograph  of  all  wires  used  in  the 
experiments. 
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Passaae  Hot-Wire  (PHW)  Probe 

Single  U-wire  probe  for  passage  measurements 


Upstream  Hot-Wire  (UHW)  Probe 


Dimensions  are  in  mm 


Figure  2.7a:  Drawings  of  custom  hot-wire  probes  used  in  experiments 
(Edited  from  Holmberg,  1996) 
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Dantec  Hotwire 


Passage  Hot-wire  (PHW)  Upstream  Hot-wire  (UHW) 


Figure  2.7b:  Pictures  of  hot-wire  probes  used  in  experiments 
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The  Dantec  hot-wire  is  mounted  on  a  90“  probe  holder  mounting  so  that  the  wire 
is  oriented  in  the  proper  direction  to  measure  streamwise  velocity  (U)  and  velocity 
fluctuations  (u').  The  probe  holder  is  mounted  on  the  sidewall  of  the  wind  tunnel.  The 
Dantec  probes  seemed  to  be  more  robust  than  the  custom  made  probes,  withstanding 
many  runs  at  high  velocity  without  breaking  the  wire.  The  disadvantage  to  this  probe  is 
that  in  order  to  move  the  probe  to  perform  streamwise  surveys  one  sidewall  of  the  tunnel 
test  section  needed  to  be  removed  to  move  the  probe  mounting  since  the  probe  holder 
extends  far  upstream  of  where  the  probe  is  mounted.  This  wire  was  primarily  used  in  low 
speed  testing  of  the  turbulence  grid  (see  Appendix  A).  The  probe  was  used  for  upstream 
measurements  in  the  transonic  wind  tunnel  prior  to  ordering  the  custom  probes. 

The  custom  made  passage  probes  (PHW)  were  primarily  designed  for  time- 
resolved  velocity  measurements  above  the  surface  of  the  turbine  blades,  but  the  probes 
were  used  for  many  of  the  mainstream  flow  measurements  as  well.  The  probes  were 
custom  made  by  Auspex  Corporation  (1416  Union  Blvd.,  Allentown,  PA  18109).  The 
probe  body  was  designed  such  that  it  would  slip  fit  into  a  compression  fitting  for  easy 
mounting  on  the  sidewalls  of  the  tunnel.  As  can  be  seen  in  Figure  2.7a,  the  prongs 
supporting  the  wire  are  at  a  90°  angle  to  the  probe  body  so  that  the  wire  could  be  oriented 
parallel  to  the  surface  of  the  turbine  blade.  The  short  length  of  the  prongs  allowed  the 
wire  to  be  moved  in  and  out  of  the  tunnel  test  section  to  different  measurement  locations 
without  removing  the  sidewalls  of  the  test  section.  The  probe  body  extending  from  the 
mounting  reduces  down  to  a  diameter  of  approximately  2.4  mm  (0.093  in.)  so  that  the 
probe  body  provides  as  little  disturbance  of  the  passage  flow  as  possible.  Figure  2.8 
shows  a  picture  of  the  UHW  mounted  such  that  the  prongs  are  above  the  heat  flux  gage. 
In  this  configuration,  the  wire  was  approximately  5  mm  above  the  surface  of  the  turbine 
blade. 

One  final  probe  (UHW)  that  was  used  in  the  experiments  was  an  Auspex  probe 
(also  shown  in  Figure  2.7a)  with  the  prongs  oriented  90°  from  the  probe  body  designed 
for  vertical  traverses  behind  the  turbulence  grid.  The  body  of  the  probe  was  sized  such 
that  it  would  fit  into  a  traverse  mechanism  mounted  at  the  top  of  the  turbine  cascade 
upstream  of  the  cascade  passages. 
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Figure  2.8:  Picture  of  PHW  probe  mounted  above  heat  flux  gage  on  blade  surface 

Calibration  of  the  hot-wire  probes  for  use  in  the  compressible  flow  field  are 
described  in  Nix,  et  al.  (2002)  in  Appendix  A  and  in  Holmberg  (1996).  Briefly,  the  wire 
is  calibrated  via  a  blowdown  calibration  method  where  the  tunnel  is  run  and  allowed  to 
blow  down  to  low  pressures.  The  cascade  will  unchoke  and  thus  the  Reynolds  number 
will  decrease  with  time  as  the  tunnel  inlet  pressure  drops.  The  wire  is  placed  near  a  pitot- 
static  probe  or  surface  pressure  transducer  to  determine  the  flow  conditions  (static 
pressure  and  Mach  number).  The  wire  voltage  output  is  measured;  along  with  wire 
temperature  and  flow  conditions,  and  a  calibration  curve  is  determined  of  the  form: 


Nu'  =  C-Re^  (2.5) 

where  Rew  is  the  wire  Reynolds  number,  Nu'  is  the  wire  Nusselt  number  and  C  and  x  are 
determined  through  a  power  law  curve  fit  of  the  blowdown  calibration  data.  Blowdown 
calibration  was  performed  with  the  wire  upstream  of  the  cascade  collocated  with  a  pitot- 
static  probe  and  also  at  each  of  the  three  measurement  locations  near  the  blade  surface. 
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The  wire  calibrations  for  each  of  these  four  locations  were  combined  to  yield  a  single 
calibration  curve.  To  minimize  error  due  to  calibration  drift,  the  wire  calibration  was 
redone  if  significant  time  passed  between  calibration  and  sampling. 

2.3.3.2  Anemometer 

The  hot-wire  probes  were  Controlled  using  a  Dantec  55M01  hot-wire 
anemometer.  The  anemometer  is  used  to  set  and  maintain  the  wire  resistance  of  the  hot¬ 
wire  probe  during  testing  and  to  tune  the  frequency  response  of  the  hot-wire.  The 
overheat  ratio  of  the  hot-wire  probes  was  typically  set  to  between  1 .5  and  1 .8.  Overheat 
ratios  in  this  range  provided  frequency  response  up  to  approximately  50  kHz.  The 
fi^quency  response  of  the  probe  was  tuned  by  applying  a  square  wave  to  the  probe  and 
adjusting  settings  on  the  anemometer.  The  response  of  the  wire  to  the  square  wave  input 
was  observed  on  an  oscilloscope.  The  overheat  ratio  was  held  constant  between 
calibration  runs  and  experiments.  The  resistance  of  the  probe  assembly  (cable,  mounting 
support,  probe  body  and  wire  prongs)  was  compensated  for  with  the  anemometer  either 
through  measurement  of  the  resistance  using  a  shorting  probe,  or,  in  the  case  of  the 
custom  probes,  through  manufacturer  supplied  data.  The  operating  resistance  of  the  hot¬ 
wire  was  set  based  on  reported  probe  resistance  at  a  specified  temperature  (typically 
22®C). 

2.3.4  Low  Frequency  Data  Acquisition-  Labview 

Low  frequency  data  acquisition  was  performed  via  a  National  Instruments  AT- 
MIO-16XE-50  digital  8  channel  A/D  board.  The  sampling  capability  of  the  data 
acquisition  card  is  increased  through  the  use  of  two  National  Instruments  AMUX-64T 
multiplexer  boards.  Each  multiplexer  board  has  32  channels,  16  channels  for  input  via 
BNC  cables,  and  16  thermocouple  input  channels.  Data  is  recorded  through  a  LabView 
code  designed  to  store  data  from  64  channels  with  the  option  of  displaying  up  to  8 
channels  immediately  after  acquiring  data.  Due  to  the  use  of  the  multiplexer,  the 
maximum  sampling  frequency  is  limited.  Since  the  data  acquisition  system  is  for  low 
frequency  or  “steady”  sampling,  the  limitations  of  sampling  frequency  were  not  an  issue. 
Data  was  typically  sampled  at  1 00  Hz  for  a  30  second  tunnel  run.  The  data  acquisition 
system  was  manually  triggered  at  the  start  of  every  run.  Data  sampled  with  the  low  speed 
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data  acquisition  system  included  all  pressures  (tunnel  total  pressure,  wall  static  pressure, 
data  from  test  section  pitot-static  probes  and  blade  surface  pressures),  tunnel  total 
temperature  and  blade  surface  temperature,  heat  flux  (HFS)  and  surface  temperature 
(RTS)  data  from  the  HFM-7  sensors  and  surface  thermocouples  during  steady  heat 
transfer  tests  and  hot-wire  data  during  blpwdown  calibration  of  the  hot-wires. 

2.3.5  High  Frequency  Data  Acquisition  -  LeCroy 

High  frequency  data  acquisition  was  performed  using  a  LeCroy  6810  waveform 
recorder.  The  LeCroy  module  used  in  this  work  consists  of  four  single  channel  modules, 
each  with  the  capability  of  sampling  data  up  to  1  MHz.  For  most  high  frequency  data 
sets,  the  sampling  frequency  was  set  at  100  kHz  (At  =  10  |Xs)  and  data  were  sampled  for 
1.31  seconds  (131,070  data  points  per  sample).  The  frequency  resolution  of  the  sampled 
data  was  then  Af  =  0.763  Hz.  Data  were  filtered  at  40  kHz.  The  voltage  limits  and 
offsets  were  set  on  the  LeCroy  to  maximize  dynamic  resolution  of  the  measured  signals. 
The  data  acquisition  was  triggered  manually  once  the  tunnel  had  reached  a  steady  inlet 
pressure.  The  high  speed  data  acquisition  system  was  used  to  sample  all  hot-wire  data 
(except  during  blowdown  calibration)  and  high  frequency  surface  heat  flux  from  the 
HFM  sensors. 

2.4  Turbulence  Grid 

Turbulence  was  generated  in  the  facility  using  a  passive  turbulence  grid  placed 
upstream  of  the  cascade  test  section  inlet  (see  Figure  2.1  and  Figure  2.3).  The  grid 
(shown  in  Figure  2.9)  consists  of  three  5.08  cm  (2  in.)  wide  bars  spaced  7.62  cm  (3  in.) 
apart.  Each  bar  is  2.54  cm  (1  in.)  thick. 
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Modified  Design:  Passive  Grid 


Original  Design:  Active  Grid 

Figure  2.9:  Turbulence  generating  grid  used  in  the  facility 

The  grid  was  developed  and  evaluated  through  extensive  testing  in  both  the 
transonic  facility  and  a  low  speed  flow  facility.  Details  of  the  development  of  the  grid 
arc  discussed  in  Nix,  et  al.  (2002)  and  the  attached  corresponding  presentation  (which 
extends  the  data  presented  in  the  paper)  in  Appendix  A.  The  grid  was  originally  designed 
to  generate  turbulence  through  active  blowing.  Testing  of  the  active  grid  in  the  low  speed 
facility  produced  the  desired  intensity  (approximately  20%)  and  scale  of  turbulence  (3-4 
cm).  The  objective  turbulence  intensity  and  length  scale  were  based  on  data  reported  by 
Van  Fossen  and  Bunker  (2002)  and  Radomsky  and  Thole  (1998).  Van  Fossen  and 
Bunker  (2002)  used  a  60°  section  of  a  GE  90  dual-annular  combustor  in  a  wind  tunnel  to 
model  the  exit  flow  of  a  combustor,  including  inlet  swirl  vanes,  film  cooling  holes  and 
dilution  holes.  All  three  of  these  combustor  components  contribute  strongly  to  the 
combustor  exit  turbulence  conditions.  Hot  wire  surveys  were  performed  downstream  of 
the  combustor  section  to  determine  turbulence  intensity  and  length  scale.  Turbulence 
intensities  as  high  as  20%  were  measured  with  length-scales  between  1  and  1.5  cm. 
Using  these  results,  the  desired  length  scale  was  increased  to  match  the  length  scale  to 
chord  length  ratio  (Ax/c)  of  the  blades  of  a  GE  90  engine  to  the  blades  used  in  the 
transonic  turbine  cascade  facility.  The  desired  turbulence  intensity  was  based  on  these 
rcsults  and  reported  decay  data  through  a  turbine  vane  passage  by  Radomsky  and  Thole 
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(1998).  The  turbulence  intensity  of  the  core  flow  believed  to  be  consistent  with  engine 
conditions  based  on  these  two  works  is  in  the  range  of  6-1 0%. 

Results  of  grid  evaluation  in  the  low  speed  rig  and  in  the  transonic  facility  are 
presented  in  detail  in  Nix,  et  al.  (2002).  When  the  turbulence  grid  was  moved  to  the 
transonic  facility,  testing  showed  that  the  active  grid  was  generating  inadequate  levels  of 
turbulence  (only  about  5%).  It  was  desired  to  match  the  momentum  flux  ratio 
(I=Pju//pooUoo^)  of  the  blowing  from  the  bars  between  the  two  facilities  in  order  to  obtain 
similar  turbulence  fields.  The  mainstream  flow  velocity  (u«))  of  the  low  speed  facility 
was  only  20  m/s  and  the  test  section  was  at  ambient  conditions.  Due  to  the  high  density 
(p«,)  and  velocity  (u<„)  of  the  mainstream  flow  in  the  transonic  facility,  the  pressure 
required  within  the  internal  plenum  of  the  blowing  bars  to  match  momentum  flux  ratio 
with  the  low  speed  testing  was  much  higher  than  achievable.  The  momentum  flux  ratio 
at  the  maximum  plenum  supply  pressure  to  the  bars  in  the  transonic  facility  was  10-12 
times  lower  than  the  low  speed  testing.  At  this  momentum  flux  ratio,  low  speed  test 
demonstrated  that  the  blowing  from  the  bars  essentially  only  acted  to  fill  in  the  wake 
from  the  bar,  reducing  the  effectiveness  of  the  grid  in  producing  elevated  levels  of 
turbulence. 

The  grid  was.  modified  to  be  a  passive  grid  with  relatively  large  blockage  area 
(approximately  40%).  The  modification  was  made  based  on  results  obtained  by  Polanka 
and  Bogard  (1999)  using  large  bars  to  generate  high  turbulence.  The  modified  grid 
operated  as  a  passive  grid  generated  an  acceptable  level  of  turbulence  with  intensity  of 
15%  and  2  cm  length  scale  near  the  cascade  test  section  inlet.  The  turbulence  intensity 
decays  to  approximately  10%  near  the  entrance  to  the  cascade  passages  with  a  length 
scale  of  2  cm.  This  level  is  representative  (if  not  slightly  higher)  than  the  levels  believed 
to  exist  in  the  core  flow  following  the  first  stage  vane  row. 

One  major  concern  with  the  modified  grid  was  uniformity  of  the  flow  downstream 
of  the  grid.  The  modified  grid  has  an  open  area  of  60%,  with  40%  flow  blockage  area, 
and  thus  there  was  concern  with  inlet  flow  non-uniformity  due  to  large  bar  wakes.  Cross¬ 
stream  traverses  of  both  a  total  pressure  probe  and  the  upstream  hot-wire  probe  (UHW) 
were  performed  behind  the  turbulence  grid,  approximately  10.2  cm  (4  in.)  upstream  of 
the  leading  edge  of  the  center  cascade  passage  in  order  to  measure  flow  and  turbulence 
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field  uniformity  behind  the  modified  grid.  Each  traverse  covered  about  half  the  passage 
in  the  pitchwise  direction.  The  results  of  the  surveys  are  shown  in  Appendix  A.  Total 
pressure  surveys  demonstrate  that  although  there  is  a  significant  drop  in  total  pressure 
across  the  grid  (approximately  7%),  the  total  pressure  of  the  flow  approaching  the  leading 
edge  of  the  cascade  is  relatively  uniform,  with  total  pressure  uniformity  of  ±1%.  The 
main  tunnel  flow  passes  through  a  2:1  area  contraction  downstream  of  the  grid  (as  shown 
in  Figure  2.3)  before  entering  the  cascade  test  section,  which  is  believed  to  help  in 
merging  and  mixing  of  the  relatively  large  bat  wakes,  thus  increasing  flow  uniformity. 

Also  of  concern  was  the  possibility  of  vortex  shedding  off  of  the  bars  generating 
coherent  structures  in  the  flow.  Analysis  of  the  power  spectral  density  (PSD)  of  the  hot¬ 
wire  signals  at  various  locations  in  the  flow  (varying  in  both  streamwise  and  pitchwise 
location)  demonstrated  that  the  spectra  were  free  of  peaks  which  would  indicate  vortex 
shedding.  The  vortex  shedding  fi^quency  of  the  bars  as  estimated  based  on  a  Strouhal 
number  of  0.13  (versus  0.21  for  circular  cylinders).  Estimating  the  expected  frequency 
provided  a  guideline  of  where  to  look  for  evidence  of  shedding  in  the  spectra.  Plots  of 
streamwise  and  spanwise  spectral  data  are  shown  in  the  presentation  in  Appendix  A. 

2.5  Measurement  Uncertainty 

Measurement  uncertainty  in  the  current  work  is  detailed  in  Appendix  B.  Table  2.3 
below  summarizes  the  total  uncertainty  determined  for  all  relevant  parameters  in  this 
work. 


Table  23  Summary  of  total  uncertainty  (A) 


Parameter 

Total  Uncertainty 

h 

(average  of  low  and  high  turbulence) 

6.3% 

U^rms 

5.0% 

Tu 

5.5% 

A, 

6.1% 

Aht 

4.5% 
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Chapters 


Data  Analysis  and  Results 


3.1  Turbulence 

3.1.1  Definition  of  Turbulence 

Turbulence,  as  defined  by  Hinze  (1975),  is  an  “irregular  condition  of  flow  in 
which  the  various  quantities  show  a  random  variation  with  time  and  space  coordinates,  so 
that  statistically  distinct  average  values  can  be  discerned.”  Further,  it  is  stated  that,  when 
viewing  an  oscillogram  (or  time  history  of  velocity  for  instance),  “given  a  point  in  a 
turbulent  domain  a  distinct  pattern  is  repeated  more  or  less  regularly  in  time.”  If  different 
turbulent  fields  are  compared,  there  may  exist  a  difference  in  the  size  (or  amplitude)  of 
the  patterns,  as  well  as  the  “violence”  of  the  fluctuations.  The  author  states: 

“...to  describe  a  turbulent  motion  quantitatively,  it  is  necessary  to  introduce  the 
notion  of  scale  of  turbulence. . .”  and 

“It  is  apparent  that  it  is  insufficient  to  characterize  a  turbulent  motion  by  its  scale 
alone,  since  to  do  so  does  not  tell  us  anything  about  the  violence  of  the  motion.  One 
cannot  take  the  average  value  of  the  velocity  as  a  measure  of  this  violence,  because  the 
violence  of  the  fluctuations  with  respect  to  this  average  velocity  is  just  what  one  wants  to 
know.” 

These  two  quantities,  the  turbulence  intensity  (violence  of  motion)  and  turbulence 
scale  (scale  of  motions)  can  be  determined  through  statistical  and  spectral  analysis  of  a 
time-resolved  measurement  of  velocity.  Experimentally,  this  is  typically  measured 
through  intrusive  measurements  such  as  hot-wire  anemometry  or  high  frequency  pressure 
measurement  devices,  or  by  non-intrusive  measurements  such  as  LDV  (Laser  Doppler 
Velocimetry).  The  following  sections  will  discuss  how  turbulence  is  characterized 
experimentally  and  will  detail  the  analysis  tools  employed  to  do  so. 
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3.1.2.  Characterizing  Turbulence 

3. 1.2. 1  Time-Resolved  Velocity  Signal  Decomposition 

In  the  current  work,  time-resolved  velocity  is  measured  using  a  single  sensor  hot¬ 
wire  probe.  The  hot-wire  probe  is  limited  in  its  measurement  of  turbulence  by  the 
frequency  response  of  the  probe  and  spatial  averaging  due  to  a  finite  length  sensing  wire. 
The  hot-wire  probe  produces  a  time-varying  voltage  signal  which  can  be  converted  to  a 
streamwise  velocity  signal  through  calibration  of  the  wire  in  a  known  flowfield  (see  Nix, 
et  al.,  2002  in'  Appendix  A).  From  this  velocity  signal,  one  can  determine  the  mean 
velocity  (C/)  and  the  fluctuating  component  of  velocity  {u).  The  measured  value  of 
velocity  at  any  given  sample  time  (m  )  can  be  decomposed  as  follows  (Sauer,  1996): 

u=U+u  (3.1) 

An  example  of  a  time  varying  signal  and  illustration  of  the  components  of  the 
velocity  signal  is  shovm  in  Figure  3.1. 


Time 


Figure  3.1:  Illustration  of  velocity  time  signal  and  decomposition 


Data  Analysis  and  Results 


33 


The  mean  velocity  is  determined  by: 


(3.2) 


where  N  is  the  number  of  discretized  data  points  in  the  sampling  period  of  length  T.  The 
fluctuating  component  of  veilocity  is  then  determined  by: 

u=u-U  (3.3) 

3.1. 2.2  Turbulence  Intensity  and  Velocity  Spectrum 

From  the  decomposition  of  the  turbulent  velocity  time  signal,  the  turbulence 
intensity  can  be  determined.  Using  the  mean  and  fluctuating  velocity  components  of  the 
signal,  the  turbulence  intensity  (Tu)  of  the  flow  at  a  given  point  is  calculated  by: 


Tu  = 


U  U 


(3.4) 


The  value  of  ,  the  root  mean  square  (RMS)  for  the  discretized  fluctuating  velocity 

signal,  can  be  determined  in  either  the  time  domain,  or  the  frequency  domain. 
Mathematically,  this  quantity  is  determined  in  the  time  domain  as  follows: 

(3-5) 

Another  term,  the  heat  flux  intensity  (Tuq),  can  be  defined  in  a  similar  way,  using  the 
time-resolved  heat  flux  signal  instead  of  the  velocity  signal. 

In  the  frequency  domain,  the  RMS  can  be  determined  through  the  power  spectral 
density  (PSD)  of  the  fluctuating  velocity  signal  («'(/))•  The  velocity  spectrum  (PSD)  is 
essentially  the  energy  at  each  frequency  of  turbulent  fluctuation  (or  eddy  scale).  The 
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integral  of  this  spectrum  gives  the  fluctuating  energy.  Based  on  the  discretized  signal, 
this  is  the  summation  of  the  energy  over  the  entire  spectrum,  which  is  equal  to  the  square 
of  the  RMS  of  the  velocity  signal.  In  the  frequency  domain,  the  RMS  is  then  determined 
as  follows: 

"1/2  -1/2 

(ecm  =  Z^(/)rAf  (3.6) 

.0  J  L/=i  J 


where  E(f)  is  the  energy  at  each  frequency  and  Af  is  the  frequency  resolution,  which  is 
equal  to  the  sampling  frequency  divided  by  the  number  of  sample  points  (N).  The  PSD  is 
determined  by: 


£(/)  = 


i^xjfyxjfS 


where  Sxx  is  the  auto-spectrum,  Xf  is  the  FFT  of  the  fluctuating  velocity  time-signal  and 
the  *  denotes  the  complex  conjugate  of  the  FFT.  The  factor  of  is  a  product  of  the  FFT 
function  in  the  Matlab  software  used  to  process  the  data. 

The  PSD  of  a  measured  turbulent  velocity  signal  reveals  much  about  the  flow. 
Through  observation  of  the  PSD  of  one  velocity  signal  compared  to  the  PSD  of  a  signal 
from  a  flow  with  different  turbulence  conditions,  assuming  they  have  similar  mean 
velocities,  one  can  easily  determine  which  has  the  higher  intensity  and  in  some  cases, 
larger  or  smaller  scale  (scale  will  be  discussed  in  the  next  section).  This  is  illustrated  in 
Figure  3.2.  The  data  in  Figure  3.2  (as  well  as  Figure  3.3  and  Figure  3.4  later)  is  from  the 
low  speed  testing  of  the  turbulence  grid.  These  plots  are  merely  for  illustration  of  the 
data  analysis  methods.  In  the  figure,  the  two  spectra  are  measured  at  different  streamwise 
locations  in  a  turbulent  flow  using  identical  sampling  frequencies  where  the  mean 
velocity  is  relatively  constant.  It  is  apparent  from  the  figure  that  the  turbulence  level  has 
decayed  from  signal  I  to  signal  2  as  integrating  under  the  PSD  curve  of  signal  1  will 
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yield  a  higher  value  of  than  integrating  under  the  PSD  of  signal  2.  With  similar 
mean  velocity,  this  then  indicates  that  signal  one  has  higher  turbulence  intensity. 
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Figure  3.2:  Example  of  velocity  power  spectral  density  (PSD) 


3.1. 2.3  Length  Scale  and  Autocorrelation 

The  length  scale  of  the  turbulent  flovvfield  used  in  this  work  is  based  on  the 
integral  length  scale  (Ax).  The  integral  length  scale  is  representative  of  the  size  of  the 
most  energetic  eddies  in  the  flow  and  is  determined  through  analysis  of  the 
autocorrelation  of  the  fluctuating  velocity  time  signal.  The  autocorrelation  (Ruu)  of  the 
fluctuating  component  of  velocity  (normalized  such  that  Ruu  is  equal  to  1  at  zero  lag, 
T=0),  was  calculated  by: 


r  N 


uXt)-uXt  +  r) 


u 


ms 


/=! 


u 


ms 


where  r  =  j-bJ. 


(3.8) 


In  the  limit  as  T  goes  to  zero  (zero  lag  time),  the  value  of  the  unnormalized 
autocorrelation  is  equal  to  the  square  of  the  RMS  of  the  velocity  signal.  The 
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autocorrelation  was  used  to  determine  the  integral  time  scale  (r  )  by  integrating  under  the 
curve  to  the  first  zero  crossing: 


No 


T=  {lHT)dT  =  ^R^r^^ 

1=0  <=' 


(3.9) 


where  No  is  the  point  of  the  first  zero  crossing. 

Theoretically,  the  curve  should  be  integrated  to  infinity  (all  data  in  the  discretized 
time  signal),  however,  in  practice  it  is  typically  only  integrated  out  to  the  first  zero 
crossing,  since  noise  may  be  present  in  the  auto-correlation.  This  noise  is  an  effect  of 
dominant  peaks  in  the  frequency  spectrum,  in  the  case  of  this  work  from  small  amplitude, 
high  frequency  vibrations  of  the  hot-wire  probe  support.  Integrating  the  autocorrelation 
out  beyond  the  first  zero  crossing  may  inflate  the  value  of  the  time  scale,  even  if  the 
oscillations  are  small.  From  analysis  of  the  autocorrelation  of  most  measurements  in  this 
work,  the  oscillations  in  the  “tail”  of  the  autocorrelation  (after  the  first  zero  crossing)  are 
of  very  small  amplitude.  If  there  were  a  problem  with  the  measurements  (excessive 
probe  vibrations  or  60  Hz  electrical  noise  for  example),  the  autocorrelation  would  show 
significant  oscillations  in  the  decaying  region  and  possibly  large  amplitude  oscillations  in 
the  tail. 

The  streamwise  integral  length  scale,  which  is  representative  of  the  largest  eddies 
in  the  turbulent  flowfield,  is  determined  by  invoking  Taylor’s  hypothesis  of  frozen 
turbulence  which  gives: 

A,  =U-7  (3.10) 

the  time  scale  is  multiplied  by  the  mean  velocity  to  yield  the  scale  of  the  most  energetic 
eddies  in  the  flow. 

Figure  3.3  shows  two  typical  velocity  autocorrelation  curves.  Assuming  that  the 
mean  velocity  is  constant  between  the  two  curves,  it  can  often  be  determined  from 
comparison  of  two  autocorrelations  which  signal  has  the  larger  scale.  A  slower  decay 
indicates  that  correlation  exists  within  the  flow  for  larger  time  delays,  and  thus  flow 
structures  (e.g.  turbulent  eddies)  are  relatively  larger.  This  function  helps  describe  time 
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scales  (and  length  scales)  present  in  the  turbulent  flowfield.  In  Figure  3.3,  integrating 
under  the  autocorrelation  curve  for  both  signals  indicates  that  signal  1  has  a  larger  time 
scale  than  signal  2.  Since  the  mean  velocity  for  both  signals  is  the  same,  the  figure 
indicates  that  signal  1  has  a  larger  length  scale  than  signal  2. 


The  function  used  to  determine  the  autocorrelation  (Equation  3.8)  and  time  scale 
(Equation  3.9)  can  also  be  used  to  determine  the  time  scale  of  the  measured  time-resolved 
heat  flux  signal  by  substituting  q’  for  u’.  Usefulness  of  the  heat  flux  time  scale  will  be 
discussed  later.  A  similar  function  can  also  be  used  to  obtain  time  domain  correlation 
(cross  correlation)  between  two  time-resolved  signals  (velocity  and  heat  flux  for 
instance). 

3.1. 2.4  Energy  Spectra  Relations  and  Isotropic  Turbulence 

As  was  indicated  earlier,  the  PSD  can  also  be  used  to  get  an  estimate  of  the 
integral  length  scale.  Referring  to  Figure  3.2,  the  turbulent  spectrum  is  shown  to  be 
relatively  flat  in  the  low  frequency  range.  At  a  certain  frequency,  the  spectrum  begins  to 
slope  downwards  and  tail  off  at  higher  frequencies.  The  region  where  the  curve  begins  to 
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slope  is  the  approximate  frequency  corresponding  to  the  most  energetic  eddies  in  the 
flowfield.  As  will  be  shown  below,  this  frequency  can  be  related  to  the  integral  scale. 

It  is  often  convenient  to  look  at  the  turbulent  velocity  spectra  in  a  different  way. 
Grid  generated  turbulence,  in  the  absence  of  outside  forces,  will  tend  towards  isotropy. 
Isotropic  turbulence  is  characterized  by  a  Cascade  of  energy  from  the  larger  eddies  to  the 
smaller  eddies  until  the  energy  is  viscously  dissipated  from  the  smallest  scales  (highest 
frequencies).  The  spectrum  of  isotropic  turbulence  is  typically  represented  by  the  Von 
Karman  relation  for  one-dimensional  spectra  (Hinze,  1975): 


^RMS  '^x 
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=  4- 

r 

1  + 

2> 
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(3.11) 


where  the  wave  number,  K,  is  equal  to  2nf  tU  and  the  most  energetic  wave  number,  Ke  is 
equal  to  3/4A, .  Figure  3.4  shows  an  example  of  a  normalized  turbulent  spectrum 

compared  to  the  theoretical  Von  Karman  spectrum.  From  the  raw,  un-normalized 
spectrum  (Figure  3.2),  the  estimated  frequency  (f)  where  the  spectrum  tails  off  (most 

3'U 

energetic  frequency)  can  be  converted  to  the  integral  scale  by  A  = - .  In  Figure  3.4, 

Snf 

this  frequency  can  be  seen  in  the  normalized  spectrum  where  k  •  A^^  =  3/4 . 

All  experimental  turbulence  data  was  normalized  and  compared  to  the  theoretical 
Von  Kdrmdn  spectrum  during  data  analysis  to  verify  that  the  flow  was  behaving  as 
expected.  The  tools  presented  in  the  preceding  sections  will  be  applied  frequently 
throughout  the  course  of  this  work. 
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Figure  3.4:  Normalized  PSD  and  Von  Karman  relation 
3.1. 2.5  Coherence 

The  correlation  of  two  simultaneously  sampled  signals  can  be  determined  in  the 
frequency  domain  through  the  coherence  function,  which  is  a  measure  of  what 
frequencies  are  correlated  between  the  two  signals.  Significant  coherence  at  a  given 
frequency  indicates  that  the  two  signals  have  energy  at  the  same  frequency,  which  is 
correlated  in  time  at  that  frequency.  The  coherence  function  is  defined  as  the  square  of 
the  cross-spectrum  between  the  two  signals  divided  by  the  product  of  the  auto-spectrum 
of  the  two  signals: 

(3-12) 

^xx'^yy 

By  definition,  the  coherence  function  is  equal  to  unity  at  all  frequencies,  however  with  a 
significant  number  of  averages,  the  coherence  function  decreases  from  a  value  of  unity  at 
each  frequency  and  asymptotically  approaches  a  value  representative  of  the  correlation  at 
each  frequency.  In  essence,  averaging  of  the  coherence  function  eliminates  uncorrelated 
content  at  each  frequency  such  that  only  correlated  content  remains.  For  the  current  data. 
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the  time-resolved  signals  of  velocity  and  heat  flux  (each  with  131,070  data  points)  were 
broken  up  into  smaller  blocks  (1024  data  points)  with  50%  overlap,  such  that 
approximately  256  averages  were  used  to  determine  the  coherence  function. 

3.2  Steady  Surface  Heat  Transfer 

This  section  presents  the  methodology  for  measuring  steady  surface  heat  transfer 
coefficients  in  the  transonic  wind  tunnel,  as  well  as  results  of  steady  measurements  with 
and  without  high  freestream  turbulence.  Surface  heat  transfer  measurements  were 
performed  at  three  different  locations  on  the  blade  surface,  one  on  the  suction  surface 
(SSI)  and  two  on  the  pressure  surface  (PSl  and  PS2).  The  locations  and  flow  conditions 
at  these  three  locations  are  described  in  section  2.2. 

3.2.1  High  Speed  Flow  Heat  Transfer  (edited  from  Nix,  1996) 

For  low  speed  flows,  the  boundary  layer  energy  equation  in  terms  of  enthalpy  is: 


P- 


u 


ok  cfy  ^ 


(3.13) 


For  compressible,  high-speed  flows,  the  kinetic  energy  of  the  flow  becomes  important 
and  the  total  enthalpy  is  used.  Also,  due  to  the  speed  of  the  flow,  frictional  heating  can 
be  significant  and  the  viscous  dissipation  term  must  be  considered.  The  energy  equation 
in  terms  of  total  enthalpy  for  high  speed  flows  then  becomes: 


<3*7,  <^7, 
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(3.14) 


When  Pr  (Prandtl  number)=l,  the  last  term  goes  to  zero  and  the  energy  equation  is  of  the 
same  form  as  the  low  speed  energy  equation  (equation  3.13)  and  the  heat  transfer 
coefficient  between  the  air  and  the  blade  is  defined  using  the  flow  total  temperature 
instead  of  the  freestream  static  temperature: 
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q”=h(T,-TJ 


(3.15) 


where  h  is  the  convection  heat  transfer  coefficient,  Tt  is  the  total  temperature  of  the  flow, 
and  Tw  is  the  surface  (wall)  temperature  of  the  blades  (measured  at  each  gage  location). 
When  the  Prandtl  number  is  not  equal  .to  one  (Pr^l),  which  is  the  general  case,  then  a 
proper  temperature  must  be  defined  for  the  flow.  The  temperature  typically  used  is  the 
adiabatic  wall  or  recovery  temperature  (Tr).  The  recovery  temperature  is  the  temperature 
the  wall  would*  attain  if  q"=0.  The  heat  transfer  is  then  given  by: 

q”=h-{T,-TJ  (3.16) 


The  recovery  temperature  (T,)  is  calculated  using: 

U  ^ 

7;=r+r-^  (3.17) 

2-C^ 

where  U-  is  the  local  freestream  velocity  of  the  flow,  and  r  is  the  recovery  factor.  The 

recovery  factor  for  laminar  flows  is  approximated  by  r  =  ,  while  r  =  Pr*^^  for 

turbulent  flow  (Kays  and  Crawford,  1993).  The  recovery  temperature  can  be  calculated 
using  the  flow  total  temperature  and  the  Mach  number  to  find  the  freestream  velocity, 
U».  Using  the  total  temperature,  the  static  temperature  can  be  calculated  using: 


T=- 


1  +  ^M' 


(3.18) 


3.2.2  Steady  Heat  Transfer  Coefficient  Methodology 

The  methodology  for  determining  heat  transfer  coefficients  in  the  transient  flows 
of  the  transonic  wind  tunnel  has  been  described  by  Popp,  et  al.  (1999)  and  in  Smith,  et  al. 
(2000).  The  following  section  outlines  the  methodology. 
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It  has  been  established  that  the  driving  force  in  convective  heat  transfer  is  the 
recovery  temperature  (Tr)  and  Equation  3.16  governs  heat  transfer.  Similar  to  Equation 
3.17  for  the  recovery  temperature,  the  total  temperature  of  the  flow  is  defined  as: 


U 

r  =r,+— =- 

'  ‘  2-C„ 


(3.19) 


The  difference  between  the  total  temperature  and  the  recovery  temperature  can  be 
presented  as: 

T-T,={\-ry^^T,  (3.20) 


Rewriting  equation  3.16: 


g"=  /j .  (r,  -  rj  =  /j  •  (t;  -  rj  -  A  •  (3.21) 

Of  the  variables  in  this  equation,  the  surface  heat  flux  (q"),  total  temperature  of  the  flow 
(Tt)  and  the  surface  temperature  (Tw)  are  directly  measured  in  the  experiments.  Due  to 
the  nature  of  the  passive  heat  exchanger  in  the  transonic  wind  tunnel,  the  total 
temperature  of  the  flow  drops  throughout  the  course  of  an  experiment  and  the  wall 
temperature  increases  as  heat  is  transferred  from  the  mainstream  flow.  Figure  3.5a  shows 
a  sample  time  history  of  the  relevant  temperatures. 
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Figure  3.5a:  Time  history  of  relevant  temperatures  (SSI) 

The  dependent  variable  in  Equation  3.21  (q")  can  be  plotted  against  the 
independent  variable  (Tt-Tw)  and  the  slope  of  the  resultant  curve  yields  the  heat  transfer 
coefficient  (h).  The  intercept  of  the  equation,  where  q=0,  results  in  Td=TrTw.  Popp,  et 
al.  (1999)  demonstrated  that  the  heat  transfer  coefficient  and  Ta  did  not  significantly 
depend  on  the  temperatures  in  Equation  3.21.  These  values  are  relatively  constant  during 
a  tunnel  run.  Figure  3.5b  shows  a  plot  of  the  experimental  data.  Linear  regression 
analysis  is  performed  on  the  data  to  determine  the  slope  and  x-axis  intercept.  With  the 
value  of  Td  known,  the  time  history  of  the  recovery  temperature  can  be  determined  by 
Equation  3.20. 
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Figure  3.5b:  Heat  transfer  data  analysis  (SSI) 

For  the  purposes  of  the  current  analyses,  the  value  of  Td  has  little  utility  other  than 
to  allow  the  time  history  of  the  heat  transfer  coefficient  to  be  plotted  from  the  time 
histories  of  q",  Tt,  and  Tw  Figure  3.6  plots  the  time  history  of  the  heat  flux  and  the 
calculated  heat  transfer  coefficient.  As  the  recovery  temperature  and  wall  temperature 
converge  (as  shown  in  Figure  3.5a)  the  heat  flux  into  the  blade  approaches  zero.  At  this 
point  the  heat  transfer  coefficient  calculations  go  unstable  due  to  the  driving  temperature 
difference  tending  towards  zero. 
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Figure  3.6:  Time  history  of  heat  Hux  and  heat  transfer  coefficient  (SSI) 

3.2.3  Steady  Heat  Transfer  Results 

Using  the  steady  heat  transfer  analysis  developed  above,  the  heat  transfer 
coefficients  were  determined  from  experimental  measurements  of  the  surface  heat  flux 
and  relevant  temperatures.  As  noted,  the  temperatures  in  the  wind  tunnel  are  transient,  so 
the  heat  transfer  is  in  essence  a  quasi-steady  measurement.  The  flow  and  heat  transfer  m 
the  wind  tunnel  never  reach  a  steady  state,  however,  the  heat  transfer  coefficient  is 
relatively  steady  throughout  the  tunnel  run.  The  results  for  measurements  with  and 
without  high  freestream  turbulence  are  discussed  below. 

3.2.3.i  Low  Turbulence  Heat  Transfer  Coefficients 

As  outlined  in  section  2.2,  three  locations  on  the  instrumented  turbine  blade  near 
the  center  of  the  cascade  were  instrumented  with  surface  pressure,  temperature  and  heat 
flux  sensors.  Two  locations  on  the  pressure  surface  and  one  location  on  the  suction 
surface  were  instrumented.  During  a  tunnel  run,  the  necessary  data  were  collected  at 
each  location  in  order  to  determine  the  heat  transfer  coefficients  at  each  location.  Several 
runs  were  performed  over  the  course  of  several  days  to  get  a  number  of  sample  points  to 
demonstrate  repeatability  and  for  uncertainty  analyses.  Table  3.1  summarizes  the  results 
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of  the  heat  transfer  analyses  with  low  turbulence  along  with  measurement  uncertainty. 
Background  or  “low  turbulence”  levels  in  the  facility  are  less  than  1%. 

The  measured  heat  transfer  coefficient  data  on  the  suction  surface  are  in  fairly 
good  agreement  with  data  reported  by  Smith,  et  al.  (2002).  The  data  on  the  pressure 
surface  were  compared  with  data  collected  on  the  pressure  surface  by  Holmberg  (1996). 
The  data  were  considerably  lower  than  reported  in  his  work,  where  heat  transfer 
coefficients  in  the  range  of  800-850  W/m^C  were  reported.  The  data  of  Holmberg 
(1996)  were  reported  to  have  large  uncertainty  in  measurements  and  flow  conditions  were 
slightly  different  with  the  smaller  turbine  blades.  The  data  collected  in  the  current  work 
are  believed  to  be  more  accurate  and  are  in  better  agreement  with  code  predictions.  In 
prior  work  in  the  transonic  facility,  calculations  were  performed  by  GE  on  the  turbine 
blades  using  the  KEP  K-eps  code  modified  from  STAN-5.  Details  of  the  KEP  code  can 
be  found  in  2^rkle  (1988).  Regarding  the  KEP  model  results,  Holmberg  (1996)  mentions 
“it  is  important  to  note  that  heat  transfer  predictions  with  this  code  are  tailored  to  match 
experience  from  full  rig  tests  while  the  measurements  here  are  made  in  a  2-D  cascade.” 
For  comparison  purposes,  the  data  collected  in  this  work  are  plotted  versus  the  GE  KEP 
code  results  in  the  next  section  (Figure  3.7),  including  results  with  high  turbulence. 


Table  3.1:  Summary  of  heat  transfer  coefficients  with  low  turbulence 


Run 

Heat  Transfer  Coefficient  < 

Number 

SSI 

PSl 

PS2 

1 

717 

644 

568 

2 

703 

623 

551 

3 

734 

655 

567 

4 

725 

651 

567 

5 

699 

649 

555 

6 

736 

691 

588 

Mean 

719 

652 

566 

STDEV 

16 

22 

13 

Uncertainty 

6.3% 

6.6% 

6.2% 

3.2.3.2  High  Turbulence  Heat  Transfer  Coefficients 

Measurements  of  steady  heat  transfer  with  high  turbulence  were  performed  and 
compared  with  the  low  turbulence  measurements  and  code  predictions.  The  turbulence 


DATA  Analysis  and  Results 


47 


grid  was  installed  and  identical  procedures  were  used  to  acquire  the  data.  The  upstream 
turbulence  field  approaching  the  cascade  blade  passages  was  characterized  with  a 
streamwise  traverse  of  a  hot-wire  probe  (see  Appendix  A)  in  order  to  determine  the 
turbulence  conditions  near  the  entrance  to  the  cascade  passages.  The  results  indicated  a 
turbulence  intensity  of  approximately  10%  with  a  length  scale  of  2.0  cm  near  the  entrance 
to  the  cascade  passages  (see  plot  on  pg.  130  of  Appendix  A).  These  conditions  are 
heretofore  referred  to  as  “high  turbulence”.  The  results  for  high  turbulence  experiments 
sampled  at  each  of  the  three  blade  measurement  locations  are  presented  in  Table  3.2 
along  with  measurement  uncertainty  and  compared  with  the  low  turbulence  results. 
Several  runs  were  performed  to  demonstrate  repeatability  and  for  uncertainty  analysis. 
The  uncertainty  in  the  increase  in  heat  transfer  coefficient  due  to  turbulence  is  only  a 
function  of  precision  uncertainty  (from  repeatability)  since  the  bias  uncertainty  of  the 
heat  flux  gages  is  present  in  both  measurements.  The  uncertainty  in  the  measured 
increase  is  shown  in  Table  3.2  and  is  an  average  of  2.9%  (much  lower  than  the  measure 
percentage  increase). 


Table  3.2:  Summary  of  heat  transfer  coefficients  with  high  turbulence 


Run 

Heat  Transfer  Coefflcienf 

SRSZSSSHI 

Number 

SSI 

PSl 

PS2 

1 

795 

776 

671 

2 

788 

768 

668 

3 

771 

761 

657 

4 

762 

750 

650 

Mean 

779 

764 

662 

STDEV 

15 

11 

10 

Uncertainty 

6.4% 

6.2% 

6.2% 

Low  Turbulence  HTC 

719 

652 

566 

%  Increase  over  Low  Tu 

8.3% 

17.1% 

16.9% 

Uncertainty  in  Increase  in  h 

2.9% 

3.2% 

2.5% 

The  results  shown  in  Table  3.2  indicate  that  there  is  a  more  significant  increase  in 
heat  transfer  coefficient  on  the  pressure  surface  due  to  high  turbulence  as  compared  to  the 
suction  surface.  The  increase  in  heat  transfer  on  the  pressure  surface  (17%)  is 
approximately  two  times  greater  than  that  measured  on  the  suction  surface  (8%).  This 
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corresponds  to  increases  in  blade  surface  temperature  of  4-10%  (between  50  and  100°C). 
Increases  in  temperature  of  25“C  or  more  can  significantly  decrease  the  life  of  turbine 
materials.  This  result  is  in  good  agreement  with  work  by  previous  researchers  on  both 
turbine  blades  and  vanes  (see  section  1.2. 1. 3).  The  increase  in  heat  transfer  on  the 
pressure  surface  gages  is  similar  in  term?  of  percentage  increase,  but  the  absolute  increase 
in  heat  transfer  coefficient  is  larger  at  location  PSl  (closer  to  the  leading  edge)  than  at 
location  PS2.  As  noted  earlier,  the  two  gages  are  located  in  regions  of  relatively  similar 
velocity;  however,  the  flow  at  PSl  experiences  a  slightly  adverse  pressure  gradient  while 
location  PS2  sees  a  favorable  pressure  gradient. 

The  results  with  low  and  high  turbulence  were  compared  with  GE  KEP  code 
predictions  for  the  turbine  blade  and  with  data  from  Holmberg  (1996)  and  Smith,  et  al. 
(2000).  The  results  are  presented  in  Figure  3.7  with  measurement  uncertainty  error  bars. 
As  shown  in  the  figure,  the  code  predicts  a  much  higher  increase  in  heat  transfer  due  to 
turbulence  on  the  pressure  surface  as  compared  with  the  suction  surface.  This  result  is 
consistent  with  the  current  data,  although  the  current  data  with  high  turbulence  agrees 
better  with  the  code  predictions  at  5%  inlet  turbulence.  It  should  be  noted  that  the  heat 
transfer  predictions  with  the  KEP  code  were  tailored  to  match  experience  fi'om  full-rig 
tests,  while  measurements  here  are  in  a  linear  2-D  cascade  (Holmberg,  1996).  At  present, 
the  KEP  code  data  are  the  only  numerical  results  available  for  this  blade  geometry  for 
comparison  with  experimental  measurements.  The  data  for  the  suction  surface  location 
(SSI)  without  turbulence  agree  well  with  the  data  from  Smith,  et  al.  (2000).  The  data 
from  Holmberg  (1996)  on  the  pressure  surface  are  much  higher  than  the  current  data. 
Holmberg  reported  very  high  uncertainty  in  his  measured  data,  as  high  as  20% 
uncertainty  (shown  in  figure). 
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Figure  3.7:  Heat  transfer  coefficient  data  compared  with  KEP  code  predictions  and 

previous  data 

3.3  Time-Resolved  Velocity  and  Heat  Transfer 

Time-averaged  (steady)  heat  transfer  results  for  the  current  work  with  and  without 
high  turbulence  have  been  presented  in  the  previous  sections.  In  addition  to  measuring 
the  effects  of  turbulence  on  the  mean  heat  transfer,  the  goal  of  this  work  was  to 
investigate  time-resolved  surface  heat  flux  and  the  interaction  of  unsteady  velocity  of  the 
mainstream  flow  near  the  blade  surface  with  surface  heat  flux.  Analysis  of  the  time- 
resolved  signals  may  provide  insight  into  the  mechanism  by  which  freestream  turbulence 
affects  surface  heat  transfer. 

As  discussed  earlier,  measurements  of  time-resolved  velocity  were  performed 
with  a  hot-wire  probe  (PHW)  located  approximately  5  mm  above  the  surface  of  the  three 
instrumented  locations  of  the  turbine  blade  (see  Figure  2.8).  This  distance  is  outside  of 
the  boundary  layer  and  far  enough  from  the  blade  surface  to  insure  the  probe  will  not 
impact  the  blade  surface  (and  destroy  the  wire)  if  severe  vibration  occurred  or  during 
startup  and  shutdown  of  the  wind  tunnel.  The  boundary  layer  thickness  was  estimated 
based  on  measurements  by  Wesner  (1996)  and  scaling  of  analyses  by  Holmberg  (1996)  to 
be  approximately  0.3  mm  over  most  of  the  suction  surface  and  between  0.3  and  1.5  mm 
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on  the  pressure  surface.  Data  from  Radomsky  (1999)  and  Ames  (1994)  demonstrated 
that  levels  of  velocity  fluctuations  outside  the  boundary  layer  but  near  the  vane  surface 
are  relatively  constant.  The  5  mm  distance  is  outside  of  the  boundary  layer,  but  close 
enough  to  measure  mean  and  fluctuating  velocities  near  the  edge  of  the  boundary  layer. 


Figure  3.8  shows  an  illustration  of  the  hot-wire  setup.  The  heat  flux  gages  at 
these  instrumented  locations  are  capable  of  high  frequency  measurements  to  obtain  a 
time-resolved  signal  of  surface  heat  transfer.  The  following  sections  discuss  the 
comparisons  of  the  signals  through  analysis  in  both  the  time  domain  and  the  frequency 
domain.  These  analyses  are  limited,  but  some  key  observations  are  made  about  how  the 
freestream  velocity  fluctuations  interact  with  surface  heat  transfer. 


Figure  3.8:  Illustration  of  simultaneous  velocity  and  heat  flux  measurements 


3.3.1  Time  Domain  Analysis 

3.3.1. 1  Direct  Time  Signal  Comparison 

Comparisons  of  the  time-resolved  velocity  and  heat  flux  signals  were  performed 
in  the  time-domain  to  investigate  the  correlation  between  the  two  signals.  Figure  3.9a  - 
Figure  3.9c  show  a  section  of  the  time-resolved  velocity  and  heat  flux  signals  normalized 
by  their  mean  values  for  all  three  measurements  locations.  The  traces  of  velocity  and 
heat  flux  demonstrate  that  there  is  low  frequency  correlation  between  the  two  signals  on 
the  pressure  surface  (PSl  and  PS2),  as  peaks  in  velocity  show  corresponding  peaks  in 
heat  flux.  The  low  frequency  energy  of  the  velocity  signal  appears  to  be  in  phase  with 
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the  corresponding  energy  in  the  heat  flux  signal.  It  is  apparent  from  the  figures  that 
higher  frequency  fluctuations  in  the  signals  do  not  appear  to  be  correlated.  These  results 
are  similar  to  data  reported  in  work  by  Moss  and  Oldfield  (1994).  From  the  plot  for  the 
suction  surface  gage  (SSI)  it  can  be  seen  that  the  fluctuations  in  velocity  are  much 
smaller  when  normalized  by  the  mean  velocity  (low  turbulence  intensity)  and  that  there 
does  not  appear  to  be  a  strong  correlation  between  the  fluctuations  in  heat  flux  with 
fluctuations  in  velocity.  The  signals  for  all  three  locations  in  Figure  3.9  are  plotted  on  the 
same  scale  for  direct  comparison  between  locations. 

Comparisons  of  the  velocity  and  heat  flux  signals  in  this  fashion  are  a  weak 
analysis  technique  except  when  there  is  a  significant  level  of  fluctuation  in  both  signals 
(as  is  the  case  with  the  pressure  surface  locations).  Frequency  domain  correlation  is  a 
much  stronger  tool  for  direct  comparison  of  two  simultaneously  sampled  signals,  which 
will  be  presented  later.  Frequency  domain  analysis  allows  for  correlation  of  the  energy 
across  the  spectrum  of  frequencies  as  well  as  phase  information. 

Using  the  analysis  tools  presented  in  Equations  3.3  -  3.5,  the  turbulence  intensity 
and  heat  flux  intensity  of  the  signals  at  each  gage  location  were  determined  and  are 
shown  in  Table  3.3.  The  turbulence  intensity  is  normalized  by  the  local  mean  velocity. 
Interestingly,  the  intensity  of  turbulent  fluctuations  (independent  of  normalization  by  the 
local  mean  velocity)  is  shown  to  be  very  low  on  the  suction  surface  and  much  higher  on 
the  pressure  surface.  The  turbulent  fluctuations  on  the  pressure  surface  are  shown  to 
decrease  slightly  as  the  flow  progresses  from  PSl  to  PS2. 

Measurements  of  cross-stream  (or  pitch-wise)  variations  of  turbulence  intensity 
downstream  of  the  grid  (presented  in  Appendix  A)  showed  little  variation  in  intensity  so 
differences  in  suction  surface  versus  pressure  surface  turbulence  is  not  due  to  non¬ 
uniformity  of  inlet  turbulence.  The  higher  intensity  of  fluctuations  on  the  pressure  side 
versus  the  suction  surface  qualitatively  agrees  with  the  data  from  Holmberg  (1996), 
although  the  magnitude  of  fluctuations  and  relative  difference  between  the  two  surfaces  is 
much  higher  in  this  case.  However,  as  pointed  out  earlier,  the  approaching  turbulence 
intensity  and  scale  are  much  higher  in  this  work.  The  low  speed  turbine  vane  turbulence 
data  from  Radomsky  and  Thole  (1998)  also  demonstrated  that  the  turbulence  levels  on 
the  pressure  surface  were  higher  than  those  measured  on  the  suction  surface. 
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Table  3.3:  Turbulence  and  heat  flux  intensity 


Gage  Location 

U'rm,  (m/s) 

Tu  (%) 

Tu,  (%) 

SSI 

931 

4.5 

0.39 

15.9 

PSl 

20.97 

18.8 

0.28 

18.6 

PS2 

17.47 

12.5 

0.38 

18.1 

From  the  table,  there  does  not  appear  to  be  any  correlation  between  turbulence 
intensity  (Tu)  and  heat  flux  intensity  (Tuq).  This  is  consistent  with  results  reported  by 
Holmberg  (1996).  The  two  pressure  surface  locations  have  very  similar  values  of  Tuq, 
however  there  is  a  fairly  signiflcant  difference  in  turbulence  intensity  between  the  two 
locations.  The  suction  surface  location  shows  a  large  value  of  Tuq,  however  the 
turbulence  intensity  at  this  location  is  very  low  as  is  apparent  by  observation  of  Figure 
3.9a.  Through  comparison  of  the  absolute  values  of  u'rms  and  q'nns  (rather  than  the  values 
of  Tu  and  Tuq  which  are  normalized  by  local  measured  mean  velocity  and  heat  flux)  there 
still  does  not  appear  to  be  good  correlation. 

If  the  value  of  u’nns  is  compared  to  the  value  of  Tuq,  there  seems  to  be  a  fairly 
linear  relationship.  In  these  tests,  the  mainstream  flow  was  not  heated,  rather  the  driving 
temperature  difference  for  heat  transfer  was  achieved  by  cooling  the  blade  surface  prior 
to  running  the  tests.  Using  unheated  mainstream  air  reduced  uncertainties  that  may  exist 
in  whether  fluctuations  in  surface  heat  flux  were  due  to  fluctuations  in  heat  transfer 
coefficient,  or  due  to  fluctuations  in  temperature.  In  addition,  unheated  mainstream  flow 
reduces  uncertainty  in  hotwire  calibration  and  conversion  of  hotwire  voltage  into  velocity 
data.  Since  there  is  little  variation  in  the  mainstream  temperature  and  the  blade  surface 
temperature  does  not  change  significantly  in  the  short  sampling  time,  the  quantity  Tuq  is 
essentially  equal  to  h^/h .  Since  there  appears  to  be  a  correlation  between  u’rms  and  Tuq, 

it  would  stand  to  reason  that  u'rms  would  be  valuable  in  predicting  increases  in  heat 
transfer.  This  has  been  demonstrated  in  many  of  the  correlations  on  flat  plates  and 
cylinders  and  will  be  used  later  in  this  work  in  the  development  of  an  analytical  model  to 
predict  increases  in  heat  transfer  due  to  turbulence. 
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PS1  Normalized  Velocity  (U/Um.n)  SSI  Normalized  Velocity  (U/U 
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b.  Location  PSl 
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Normalized  Heat  Flux  (q"/q”m««n)  SSI  Normalized  Heat  Flux  (q^/q' 


1.8 


c.  Location  PS2 

Figure  3.9:  Time-resolved  velocity  and  heat  flux  signals 

3.3.1. 2  Time  Correlation  and  Scales 

The  autocorrelation  function  for  both  the  velocity  and  heat  flux  signals  was 
calculated  to  determine  the  time  scale  of  the  fluctuations  in  each.  Figure  3.10a  -  Figure 
3.10c  present  the  autocorrelation  results  for  each  gage  location.  As  shown  in  the  figures, 
the  measured  velocity  and  surface  heat  flux  have  very  similar  time  scales  at  the  two 
pressure  surface  locations  (PSl  and  PS2).  The  suction  surface  location  (SSI) 
demonstrates  a  smaller  time  scale  for  the  heat  flux  signal  as  compared  to  the  velocity 
signal.  The  results  of  calculated  velocity  and  heat  flux  time  scales  and  streamwise 
integral  length  scale  of  velocity  are  presented  in  Table  3.4. 


Table  3.4:  Time  scales  and  integral  length  scales 


Gage  Location 

Tu  (ms) 

Ax  (cm) 

Tq  (ms) 

SSI 

0.18 

3.46 

0.14 

PSl 

0.21 

2.39 

0.20 

PS2 

0.18 

2.51 

0.18 
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As  the  table  shows,  the  integral  length  scale  calculated  on  the  suction  surface  is 
larger  than  the  pressure  surface  locations.  The  approaching  length  scale  is  approximately 
2  cm  and  increases  slightly  as  the  flow  proceeds  along  the  two  pressure  surface  locations. 
There  is  a  fairly  large  uncertainty  in  determination  of  the  integral  length  scale  (as  will  be 
presented  in  Appendix  B),  but  the  scale  determined  for  the  suction  surface  location  is 
considerably  larger  than  the  measured  upstream  and  pressure  surface  length  scales.  The 
time  scales  are  relatively  similar,  but  due  to  the  high  velocity  in  the  suction  surface,  the 
calculated  length  scale  is  higher.  Measurements  on  the  suction  surface  have 
demonstrated  smaller  turbulence  fluctuations  and  larger  scale  than  the  measurements  on 
the  pressure  surface. 
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PS2  Autocorrelation  •  R^q  PS1  Autocorrelation  -  Ryui  Rq 
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Figure  3.10:  Autocorrelation  of  velocity  and  heat  flux  signals 
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3.3.2  Frequency  Domain  Analysis 
3.3.2.1  Velocity  and  Heat  Flux  Spectra 

The  time  resolved  measurements  of  velocity  and  heat  flux  were  converted  into  the 
frequency  domain  to  determine  the  velocity  frequency  spectra  above  the  heat  flux  gages 
and  the  heat  flux  spectra  at  each  location.  Figure  3.1  la  —3.1  Ic  show  the  velocity  and  heat 
flux  spectra  for  each  gage  location.  Figure  3.12  shows  a  comparison  of  the  un¬ 
normalized  velocity  spectra  for  all  three  measurement  locations  for  comparison  of  the 
energy  at  each,  location,  while  Figure  3.13  shows  a  comparison  of  the  heat  flux  spectra  at 
all  three  measurement  locations.  Figure  3.14  shows  the  velocity  spectra  from  the  three 
locations  normalized  and  plotted  against  the  theoretical  Von  Karmdn  relation  for  one¬ 
dimensional  spectra. 


a.  Location  SSI 
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PS2  Velocity  and  Heat  Flux  PSD  PSI  Velocity  and  Heat  Flux  PSD 


Frequency  (kHz) 

b.  Location  PSI 


Frequency  (kHz) 

c.  Location  PS2 

Figure  3.11:  Velocity  and  heat  flux  spectra 
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PSD  of  Heat  Flux  PSD  of  Velocity 


Figure  3.14:  Normalized  velocity  spectra  and  Von  Karmen  relation 


From  Figure  3.11a  -  Figure  3.11c  it  can  be  seen  that  the  velocity  and  heat  flux 
spectra  have  very  similar  shapes  and  have  significant  energy  at  low  frequencies  which 
attenuates  at  frequencies  of  approximately  5  kHz  and  above.  The  spikes  in  the  velocity 
spectra  between  10  and  15  kHz  are  due  to  vibration  of  the  hot-wire  probe  support. 
Vibration  of  the  hot-wire  probe  at  this  frequency  would  be  driven  by  vortex  shedding 
from  the  probe  support.  The  shedding  frequency  of  a  body  can  be  described  by  the 
dimensionless  parameter  of  Strouhal  (S), 


S  = 


fs-D 

U 


(3.22) 


where  f*  is  the  vortex  shedding  frequency,  D  is  the  diameter  of  the  body  (probe  support), 
and  U  is  the  mean  velocity  approaching  the  probe  support.  The  Strouhal  number  is  a 
ftinction  of  the  geometry  of  the  body  and  Reynolds  number.  For  a  circular  cylinder  in 
cross-flow,  the  Strouhal  number  has  been  demonstrated  to  be  0.21  (Blevins,  1994). 
Using  a  value  of  S=0.21,  the  mean  velocity  (U)  near  the  surface  of  each  gage  location  and 
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the  probe  support  diameter  (D)  of  3.2mm,  the  calculated  shedding  frequency  ranges  from 
9  kHz  at  PSl  to  14  kHz  at  SSI .  These  calculated  values  agree  with  results  in  Figure  3.1 1 . 
This  peak  in  the  spectra  is  much  higher  at  the  suction  surface  location  (SSI)  due  to  higher 
velocity  at  this  location.  The  probe  vibration  is  not  an  issue  since  it  is  in  a  region  of  the 
spectrum  where  most  of  the  energy  has  been  attenuated  and  does  not  affect  calculations 
of  the  RMS  value  or  coherence  between  the  velocity  and  heat  flux  signals  (discussed  in 
the  following  section). 

The  date  in  Figure  3.12  indicate  the  relative  levels  of  RMS  velocity  fluctuations  at 
the  three  locations  (similar  to  Figure  3.9  in  the  time  domain).  As  shown,  the  fluctuating 
velocity  spectrum  for  the  suction  surface  location  (SSI)  has  much  less  energy  than  the 
two  pressure  surface  locations  (PSl  and  PS2).  This  agrees  with  the  date  presented  in 
Table  3.3.  The  spectra  were  noimalized  using  the  local  mean  velocity,  RMS  of 
fluctuating  velocity  and  length  scale  as  shown  in  Figure  3.14  and  compared  to  the 
theoretical  Von  Kdrman  relation  for  one-dimensional  spectra  (Equation  3.11).  The  date 
fit  well  with  the  Von  Kdrmdn  equation  except  for  a  deviation  from  the  characteristic  -5/3 
slope  region  between  10  and  20  kHz  (approximately  A^f/U  of  2-3).  As  mentioned,  this  is 
likely  due  to  vibration  of  the  hot-wire  probe  from  vortex  shedding. 

3.3.2.2  Coherence 

As  discussed  in  section  3. 1.2.5,  correlation  of  two  simultaneously  sampled  signals 
in  the  frequency  domain  is  accomplished  through  the  coherence  function.  The  coherence 
function  was  applied  to  the  time-resolved  velocity  and  heat  flux  date  in  the  frequency 
domain  to  determine  which  frequencies  were  correlated  between  the  two  signals.  The 
coherence  function  for  each  gage  location  is  shown  in  Figure  3.15. 
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Frequency  (kHz) 


Figure  3.15:  Coherence  between  velocity  and  heat  flux  signals 


The  figure  demonstrates  that  significant  coherence  at  all  three  gage  locations  is 
limited  to  a  relatively  narrow  band  of  frequencies  below  2  kHz  with  zero  coherence  at 
frequencies  starting  at  about  5  kHz.  Shown  on  the  figure  is  a  band  of  frequencies 
corresponding  to  the  integral  length  scale  range  at  the  three  gage  locations  ±10%.  The 
frequency  corresponding  to  the  integral  length  scale  at  each  location  was  determined 
using  the  analysis  presented  in  section  1.3.2 .4  where  frequency  was  converted  to  the 


integral  scale  by  A,  = 


-j-U 

Snf  ‘ 


At  each  location,  the  frequency  corresponding  to  the 


integral  length  scale  was  then  determined  by: 


3-U 


SitA, 


(3.23) 


For  location  SSI,  this  frequency  is  approximately  650  Hz,  and  570Hz  and  700  Hz  at 
locations  PSl  and  PS2,  respectively. 
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The  frequency  range  and  magnitude  of  significant  coherence  agrees  very  well 
with  data  reported  by  Holmberg  (1996).  The  coherence  measured  on  the  pressure  surface 
at  PSl  and  PSl  shows  no  “preferential  frequencies”,  rather  a  broad  spectrum  of 
frequencies  centered  on  the  integral  scale,  out  to  approximately  2  kHz,  show  significant 
coherence.  On  the  suction  surface,  the  coherence  shows  a  peak  at  approximately  550  Hz. 
This  peak  is  well  correlated  with  the  frequency  of  the  integral  scale.  The  comparison  of 
the  coherence  on  the  two  surfaces  may  indicate  that  a  much  smaller  range  of  coherent 
structures  in  the  freestream  penetrate  the  boundary  layer  on  the  suction  surface  to  affect 
surface  heat  transfer.  On  the  pressure  surface,  the  range  of  turbulent  vortices  in  the 
freestream  which  show  significant  coherence  with  surface  heat  transfer  is  broader, 
possibly  indicating  the  ability  of  lower  energy  scales  to  penetrate  into  the  boundary  layer 
to  affect  surface  heat  transfer. 

3.4  Analytical  Model 

Using  the  time-resolved  velocity  and  heat  flux  data  and  results  (both  steady  heat 
transfer  and  frequency  domain  results),  an  analytical  model,  which  may  be  used  to 
predict  increases  in  steady  heat  transfer  due  to  freestream  turbulence,  was  developed. 
The  following  sections  document  the  development  of  the  model  and  present  results  of 
application  of  the  model  to  the  current  data,  as  well  as  to  data  from  the  work  of  other 
researchers  on  varying  geometries  and  under  different  flow  conditions. 

3.4.1  Analytical  Model  Development 

The  analytical  model  developed  is  based  on  the  concept  of  large-scale  eddies 
penetrating  the  boundary  layer  and  directly  impacting  surface  heat  transfer.  Results  from 
analyses  of  the  time-resolved  velocity  and  heat  flux  data  in  the  time  and  frequency 
domains  have  demonstrated  that  the  large-scale  (low  frequency)  eddies  present  in  the 
freestream  are  highly  correlated  with  surface  heat  transfer.  Work  by  Moss  and  Oldfield 
(1994)  demonstrated  through  time-resolved  measurements  of  velocity  and  heat  flux  on  a 
flat  plate  that  eddies  in  the  boundary  layer  convected  at  the  freestream  velocity,  rather 
than  at  the  ~0.8U  characteristic  of  boundaiy  layers.  They  concluded  through  time 
domain  and  spectral  analysis  of  hot-wire  and  heat  flux  gage  data  and  through  correlation 
analysis  between  heat  flux  gages  that  the  turbulent  eddy  structure  of  the  boundary  layer  is 
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dominated  by  freestream  turbulence  structures  and  that  heat  transfer  enhancement  is 
caused  by  these  large  eddies  being  carried  deep  down  into  the  boundary  layer. 

The  model  developed  in  the  current  work  assumes  that  freestream  eddies 
penetrate  to  and  directly  transfer  heat  to  the  surface.  The  resultant  heat  transfer  is 
superimposed  onto  the  normal  boundary  layer  heat  transfer.  The  model  assumes  that  the 
heat  transfer  process  is  linear.  Figure  3.16  illustrates  the  idea  behind  the  model. 


Figure  3.16:  Illustration  of  turbulent  heat  transfer  model 

Assuming  the  eddies  have  indeed  penetrated  the  thermal  boundary  layer  and 
reached  the  surface  of  the  turbine  blade,  a  model  of  how  the  heat  “packaged”  in  these 
eddies  is  transferred  to  the  surface  was  needed.  The  model  employs  the  concept  of  the 
semi-infrnite  medium,  whereby  a  “slug”  of  hot  gas  with  a  temperature  representative  of 
the  freestream  (Too)  impacts  the  surface  and  resides  for  a  short  time,  transferring  heat  to 
the  constant  temperature  wall  (Tw)  as  a  semi-infinite  medium  of  gas.  The  equation  for 
heat  transfer  for  a  semi-infinite  medium  with  constant  surface  temperature  is  given  by 
(Incropera  and  DeWitt,  1990): 
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where  k  is  the  thermal  conductivity  and  a  is  the  thermal  diffusivity  of  the  hot  gas  slug. 
Assuming  that  the  increase  in  heat  transfer  due  to  freestream  turbulence  over  measured 
heat  transfer  without  turbulence  was  due  to  the  combined  effect  of  these  eddies 
enhancing  surface  heat  transfer,  equation  3.24  is  then  equal  to: 


where  the  thermal  diffusivity  (a)  and  thermal  conductivity  (k)  are  based  on  the 
representative  eddy  temperature,  Te  (static  temperature  of  flow).  The  model  results  yield 


infinite  heat  flux  at  time  t=0,  followed  by  heat  flux  which  decays  as  -7= .  The  increase  in 

yt 


heat  transfer  (AqO  over  low  turbulence  measurements  is  determined  by  averaging  q  over 
the  time  of  the  event  which  is  estimated  as  t=A,/u'n„s,  where  A,,  and  u’^^  are  measured 
experimentally  with  a  hot-wire.  This  time  is  essentially  the  residence  time  of  a 
freestream  eddy  that  has  penetrated  to  the  surface.  Using  the  values  of  A^  and  u'n„s  to 
determine  this  time  assumes  that  all  of  the  energy  contained  within  the  freestream 
turbulence  field  (u’n„s)  is  packed  into  a  representative  sized  eddy,  which  for  the  model  the 
most  energetic  eddies  (A^)  were  used.  The  increase  in  heat  transfer  coefficient  over  low 
turbulence  heat  transfer  is  then  determined  by  the  following  series  of  equations: 


Aq,=k- 


(T.-TJ 


ni/2 


mx- 


(3.26) 


Aq,=Ahr(T^-TJ 


(3.27) 
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ni/2 


Tta- 


(3.28) 


so  the  increases  in  heat  transfer  (Aht)  can  then  be  determined  through  measurements  of 
steady  heat  transfer  at  low  turbulence  conditions  (ho)  and  measurement  of  the  length  scale 
(A^)  and  RMS  of  turbulent  fluctuations  (u'n„s)  along  with  knowledge  of  a  representative 
temperature  of  the  flow.  The  representative  eddy  temperature  could  be  the  static,  total  or 
recovery  temperature  of  the  flow,  but  since  the  eddy  temperature  is  only  important  in 
determining  the  values  of  the  thermal  conductivity  and  diffusivity,  variations  between 
these  temperatures  would  only  have  a  small  effect  on  the  results  of  the  model 
calculations.  The  turbulent  heat  transfer  coefficient  could  then  be  determined  by  adding 
the  increases  predicted  by  the  model  to  the  measured  heat  transfer  at  low  turbulence: 

+  (3.29) 

Several  key  points  can  be  made  about  the  model: 

•  The  model  is  independent  of  the  geometry  of  the  surface,  although  both  the 
turbulent  fluctuations  (u'nns)  and  the  length  scale  (A^)  may  be  affected  by  surface 
geometry. 

•  The  model  is  independent  of  the  state  of  the  boundary  layer.  The  model  does  not 
rely  on  any  knowledge  of  whether  the  boundary  layer  is  laminar,  transitional  or 
turbulent,  since  the  primary  assumption  is  that  the  large-scale  ffeestream  turbulent 
structures  from  the  ffeestream  superimpose  on  the  boundary  layer  heat  transfer 
without  high  ffeestream  turbulence. 

•  The  model  is  independent  of  the  boundary  layer  thickness,  but  assumes  that  the 
eddy  scales  are  much  larger  than  the  boundary  layer  thickness. 

•  The  model  assumes  that  the  increase  in  heat  transfer  due  to  turbulence  over  the 
existing  heat  transfer  with  low  freestream  turbulence  is  a  linear  process  (linear 
superposition). 
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These  points  simply  state  how  the  model  is  applied  to  predict  surface  heat 
transfer.  The  ability  of  the  model  to  accurately  predict  increases  in  surface  heat  transfer 
due  to  freestream  turbulence  may  be  affected  by  geometry  and  boundary  layer 
characteristics,  however,  the  model  is  not  a  function  of  these  parameters.  In  the 
following  sections,  the  analytical  model  will  be  applied  to  the  current  data  set  and  results 
presented  in  sections  3.2.3. 1  and  3.2.3.2'and  also  to  the  data  of  other  researchers. 

3.4.2  Single-scale  Model  Results 

The  analytical  model  in  the  form  developed  in  the  previous  section  is  referred  to 
in  this  text  as  the  “single-scale”  model.  The  model  uses  integral  parameters  of  the 
freestream  turbulence  energy  and  length  scale  to  determine  the  increase  in  steady  heat 
transfer.  Later  a  variation  of  the  model  will  be  developed  which  uses  the  frequency 
domain  data  and  will  be  referred  to  as  the  “multiple-scale”  model.  The  single-scale 
model  was  applied  to  the  current  data  set  using  measured  steady  heat  transfer  coefficients 
with  and  without  turbulence  and  the  measured  turbulence  at  each  gage  location.  The 
results  are  presented  below. 

3.4.2.1  Application  to  Transonic  Facility  Data 

Using  the  measured  turbulence  data  presented  in  Table  3.3  and  Table  3.4  and 
determining  the  freestream  static  temperature  at  each  gage  location  from  the  measured 
test  section  total  temperature  and  the  blade  Mach  number  distribution,  the  single-scale 
model  was  applied  to  several  data  runs  for  each  gage  location  on  both  the  pressure  and 
suction  surfaces  of  the  transonic  turbine  blades.  Table  3.5  summarizes  the  data  from 
Table  3.1  through  Table  3.4  and  applies  the  model  equation  (Equation  3.28)  to  the  data  to 
predict  the  turbulent  heat  transfer  coefficient  at  each  location.  The  length  scales  listed  in 
Table  3.5  compared  to  the  boundary  layer  thickness  estimates  detailed  in  section  3.3 
correspond  to  a  ratio  of  length  scale  to  boundary  layer  thickness  (Ax/8)  of  approximately 
100  on  the  suction  surface  and  20-30  on  the  pressure  surface. 

Due  to  the  nature  of  the  measurement  techniques  in  the  transonic  facility,  the 
steady  heat  transfer  coefficients  with  and  without  high  freestream  turbulence  and  the 
time-resolved  velocity  and  heat  flux  data  were  measured  in  separate  runs.  The  steady 
heat  transfer  runs  were  performed  with  heated  flow  and  were  typically  of  20-30  seconds 
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in  length.  The  runs  with  the  hot-wire  in  place  were  performed  without  heated  flow,  but 
with  a  pre-cooled  turbine  blade  for  the  driving  temperature  difference  for  heat  transfer,  as 
discussed  earlier.  Separating  the  runs  was  necessary  for  several  reasons.  Besides  the 
desire  for  unheated  flow  when  performing  hot-wire  measurements,  the  run  times  of  the 
hot-wire  tests  were  of  short  length  to  minimize  the  possibility  of  breaking  the  hot-wire 
which  is  extremely  fragile,  especially  in  such  an  aggressive  flowfield.  In  addition,  in 
order  to  measure  steady  and  time-resolved  heat  transfer  simultaneously,  both  low  and 
high  frequency  data  acquisition  were  required.  Recording  data  with  the  two  data 
acquisition  systems  simultaneously  increases  the  risk  of  cross  talk  and  noise  in  the  time- 
resolved  data,  increasing  uncertainty  in  the  data.  The  data  in  Table  3.5  then  consists  of 
averaged  heat  transfer  coefficient  data  without  turbulence  from  Table  3.1,  averaged 
turbulent  heat  transfer  coefficient  data  from  Table  3.2  and  u'^^  and  data  recorded  from 
time-resolved  data  runs  with  the  hot-wire  collocated  with  the  heat  flux  gage  (Table  3.3 
and  Table  3.4)  along  with  uncertainty  in  the  analytical  model  predicted  increases. 


Table  3.5:  Results  of  single-scale  model  applied  to  transonic  tunnel  data 


Meas. 

Location 

Measured 
tt’rms 
m/s  - 

Measured 

Ax 

cm 

Ah, 

W/m^K 

LowFST 

b 

W/m^K 

Predicted 

Turbulent 

h 

W/m*K 

Uncertainty 
In  Modeled 
h 

W/m^K 

Measured 

Turbulent 

h 

W/m*K 

9.31 

3.81 

57 

719 

776 

3 

779 

SSI 

9.09 

3.51 

58 

719 

777 

3 

779 

8.56 

3.08 

60 

719 

779 

3 

779 

20.97 

2.34 

116 

652 

768 

5 

764 

PSl 

21.79 

2.41 

116 

652 

768 

5 

764 

21.30 

2.34 

117 

652 

769 

5 

764 

21.29 

2.48 

114 

652 

766 

5 

764 

PS2 

17.47 

2.54 

101 

667 

5 

16.90 

2.48 

100 

5 

The  results  demonstrate  very  good  agreement  between  the  model  predicted  and 
average  measured  turbulent  heat  transfer  coefficients  on  both  blade  surfaces.  In  addition, 
the  analytical  model  closely  predicts  the  observed  differences  in  heat  transfer 
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augmentation  between  the  suction  surface  and  pressure  surface.  These  results  have 
demonstrated  that  a  relatively  simple  analytical  model  with  a  physical  basis  (penetration 
of  large  scale  eddies  deep  into  the  boundary  layer)  coupled  with  a  theoretical  heat  transfer 
mechanism  (semi-infinite  medium  assumption)  is  very  successful  at  predicting  increases 
in  heat  transfer  due  to  freestream  turbulence.  The  model  requires  only  the  knowledge  of 
the  local  turbulence  characteristics  (velocity  fluctuations  and  scale)  near  the  surface  of 
the  turbine  blade.  In  order  to  further  test  the  validity  of  the  model,  its  application  to 
additional  data  sets  needs  to  be  explored.  The  following  sections  continue  the  application 
of  the  model  to  the  data  sets  of  other  researchers. 

3.4.2.2  Application  to  Other  Data  Sets 

With  the  model  data  applied  to  the  current  data  set  with  extremely  good  results, 
the  model  was  then  applied  to  the  work  of  other  researchers  to  test  its  validity  against 
other  data  sets.  These  data  sets  include  several  different  geometries  and  both  laminar  and 
turbulent  boundary  layers.  According  to  the  primaiy  assumptions  of  the  model,  the 
geometry  and  state  of  the  boundary  layer  should  not  affect  the  ability  of  the  model  to 
predict  increases  in  heat  transfer  due  to  turbulence.  Requirements  for  application  of  the 
model  were  that  the  experiments  measured  surface  heat  transfer  coefficients  with  and 
without  turbulence  and  included  local  measurements  (i.e.  collocated  with  the  surface  heat 
transfer  measurements)  of  time-resolved  velocity  to  determine  u'^^^  and  Knowledge 
of  the  freestream  temperature  reduces  uncertainty  in  the  modeling,  but  estimated  values 
based  on  reported  data  should  be  adequate,  since  this  quantity  only  affects  the  values  of 
the  thermal  conductivity  and  thermal  diffusivity  according  to  Equation  3.28,  which  have 
a  second  order  effect  on  the  model  results.  The  local  velocity  measurements  in  these 
works  were  performed  by  either  hot-wire  anemometiy  or  LDV  measurements.  The  data 
sets  that  were  used  for  comparison  with  the  single-scale  model  were  as  follows: 

•  Stagnation  point  data  -  Diller  and  Van  Fossen,  1 996 

•  Low  speed  turbine  vane  data  -  Radomsky  and  Thole,  1 998 

•  Low  speed  turbine  vane  pressure  surface  data  -  Ames,  1 994 

•  Turbulent  flat  plate  data  -  Moss  and  Oldfield,  1 992 

•  Transonic  turbine  blade  data  -  Holmberg,  1 996 
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The  model  application  in  these  cases  spans  cylinder  stagnation  point  data,  flat  plat 
data,  turbine  vane  and  turbine  blade  data,  and  includes  both  laminar  and  turbulent 
boundary  layers  on  different  geometries.  The  following  sections  document  the 
application  of  the  model  to  these  cases  and  present  the  results,  in  some  cases  discussing 
limitations  of  the  model  based  on  these  results. 

3.4.2.3  Stagnation  Point  Data  -  Diiier  and  Van  Fossen,  1996 

The  effects  of  freestream  turbulence  on  stagnation  point  heat  transfer,  measuring 
time-resolved  velocity  and  heat  flux,  were  studied  in  unpublished  work  by  Diller  and  Van 
Fossen  (1996).  The  work  was  performed  in  the  same  facility  as  published  work  by  Van 
Fossen,  et  al.  (1995),  using  the  same  facility,  turbulence  grids  and  heat  transfer  models. 
The  heat  transfer  models  had  an  elliptical  leading  edge  with  a  wedge  shaped  afterbody 
extending  downstream  to  eliminate  vortex  shedding.  The  heat  transfer  model  was  placed 
in  a  low  speed  wind  tunnel  with  a  freestream  velocity  of  approximately  45  m/s.  The 
cylinder  model  had  an  HFM  high  frequency  heat  flux  microsensor  mounted  in  the 
stagnation  region  and  a  hot-wire  anemometer  probe  mounted  such  that  it  could  be 
traversed  in  the  streamwise  direction  to  measure  time-resolved  velocity  approaching  the 
stagnation  region.  The  turbulence  grid  was  placed  upstream  of  the  heat  transfer  models 
and  consisted  of  a  biplanar  grid  of  1/16  inch  (0.159  cm)  bars  spaced  0.72  cm  apart.  The 
results  of  the  heat  transfer  tests  for  this  grid  (G4  in  the  work  of  Van  Fossen,  et  al.)  are 
presented  in  Table  3.6  below. 


Table  3.6:  Diller  (1996)  heat  transfer  data 


Case 

Distance 

(in) 

Mean 

Velocity 

(m/s) 

Mean 

Heat 

Flux 

(W/m^) 

RMS 

Velocity 

(m/s) 

Measured 

HTC 

(W/m^K) 

Tu 

% 

Lx 

(cm) 

G4201 

1.89 

32.18 

4861 

2.4295 

262 

7.55 

0.86 

G4202 

0.27 

15.15 

4925 

2.3937 

267 

15.80 

0.65 

G4203 

0.22 

13.31 

4971 

2.5867 

272 

19.43 

0.98 

G4204 

0.17 

11.09 

5149 

2.7481 

271 

24.78 

0.77 

G4205 

0.12 

8.53 

5403 

2.6242 

274 

30.78 

0.42 

G4206 

0.07 

6.06 

5504 

2.2664 

273 

37.41 

0.24 
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The  data  in  Table  3.6  represents  six  cases  where  the  hot-wire  probe  was  located  at 
different  streamwise  locations  (distance  in  second  column)  upstream  of  the  stagnation 
region  of  the  heat  transfer  model.  The  mean  velocity  decreases  as  the  wire  is  moved 
closer  to  the  stagnation  region.  The  streamwise  turbulence  intensity  is  normalized  by  the 
mean  velocity,  so  the  turbulence  intensity  increases  at  locations  nearer  the  stagnation 
region.  The  streamwise  velocity  fluctuations  are  relatively  constant  approaching  the 
stagnation  region  of  the  blade,  while  the  length  scale  decreases.  The  analytical  model 
was  applied  to  these  data  and  the  predicted  heat  transfer  coefficients  with  turbulence  were 
compared  to  measured  data.  The  results  are  shown  in  Table  3.7.  The  mean  heat  transfer 
coefficient  at  low  turbulence  (less  than  0.5%  for  this  facility)  was  approximately 
180  W/m^K.  Uncertainty  in  the  model  predictions  is  likely  of  the  same  order  as  the 
model  predictions  from  the  current  research  work  (±4.5%  on  the  value  of  Aht). 

As  the  results  show,  the  model  predictions  of  the  turbulent  heat  transfer 
coefficient  are  quite  good.  The  location  where  the  turbulence  characteristics  are 
measured  closest  to  the  stagnation  region  (0.07  inches)  provides  the  best  prediction  of 
heat  transfer  increases  due  to  turbulence,  indicating  the  dependence  of  the  model  on  local 
turbulence  measurements,  rather  than  measurements  far  upstream.  In  the  case  of  a  flat 
plate  or  turbine  vane  or  blade,  this  would  translate  to  measurements  near  the  surface, 
close  to  the  boundary  layer  edge,  rather  than  far  out  into  the  freestream. 


Table  3.7:  Results  of  single-scale  model  applied  to  Diller  (1996)  data 


Case 

Measured 

u’rms 

(m/s) 

Measured 

Ax 

(cm) 

Aht 

(W/m^K) 

Predicted 

Turbulent 

h 

(W/m^K) 

Measured 

Turbulent 

h 

(W/m^K) 

G4201 

2.4295 

0.86 

52.56 

180 

233 

262 

G4202 

2.3937 

0.65 

60.01 

180 

240 

267 

G4203 

2.5867 

0.98 

50.80 

180 

231 

272 

G4204 

2.7481 

0.77 

59.07 

180 

239 

271 

G4205 

0.42 

180 

258 

274 

G4206 

0.24 

180 

276 

273 
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3.4.2.4  Low  speed  turbine  vane  data  -  Radomsky  and  Thole,  1998 

The  effect  of  freestream  turbulence  intensity  and  length  scale  on  turbine  vane  heat 
transfer  in  a  low  speed  facility  was  studied  in  detail  by  Radomsky  and  Thole  (1998).  The 
study  focused  on  documenting  the  surface  heat  transfer  and  highly  turbulent  flowfield 
around  a  stator  vane  placed  in  a  large-scale  wind  tunnel.  The  test  section  consists  of  a 
stator  vane  and  the  adjacent  blade  leading  edges  (forming  two  vane  passages)  to  correctly 
simulate  the  pressure  distribution  and  secondary  flows  about  the  vane.  An  active 
turbulence-generating  grid  was  placed  upstream  of  the  test  section  to  simulate  the 
turbulence.  Three-component  laser  Doppler  velocimetry  (LDV)  measurements  were 
performed,  mapping  the  turbulence  intensity  in  the  two  passages.  The  turbulent  length 
scale  was  measured  using  a  single-sensor  hot-wire.  In  order  to  accurately  apply  the 
analytical  model,  the  measurements  nearest  to  the  vane  surface  (near  the  edge  of  the 
boundary  layer)  were  used  as  input  data.  Two  inlet  turbulence  conditions  (10%  and  20%) 
were  compared.  The  integral  length  scale  of  the  inlet  turbulence  was  4.6  cm  for  the  10% 
turbulence  case  and  5.0  cm  for  the  20%  turbulence  case.  Local  length  scale 
measurements  (matching  the  LDV  measurement  locations)  were  not  performed,  but 
measurements  along  a  streamline  in  the  center  of  the  vane  passage  demonstrated  that  the 
length  scale  did  not  change  significantly  throughout  the  vane  passages.  A  constant  length 
scale  was  then  assumed  for  the  model  application.  The  model  was  used  to  predict  the 
increase  in  heat  transfer  due  to  turbulence  by  Equation  3.28,  then  the  value  of  Aht  was 
added  to  the  measured  Stanton  number  for  each  chord  wise  location  on  the  vane.  The 
results  of  the  model  application  are  shown  graphically  in  Figure  3.17  and  3.18  with 
uncertainty  bars  of  ±4.5%  added  to  selected  data  points  for  the  values  of  the  Stanton 
number  and  ±2.4%  applied  to  the  St/Sto  data. 
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Figure  3.17:  Comparison  of  model  predicted  and  measured  Stanton  numbers  (St) 

from  Radomsky  and  Thole  (1998) 
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Figure  3.18:  Comparison  of  model  predicted  and  measured  St/Sto  from  Radomsky 

and  Thole  (1998) 


The  results  shown  in  Figure  3.17  demonstrate  that  the  model  is  fairly  good  at 
predicting  the  increases  in  heat  transfer  due  to  turbulence  for  both  the  low  (10%)  and 
high  (20%)  turbulence  cases  in  the  region  where  the  boundary  layer  is  laminar.  The 
model  slightly  over  predicts  the  heat  transfer  for  the  10%  turbulence  case  and  slightly 
under  predicts  the  heat  transfer  for  the  20%  turbulence  case  over  most  of  the  blade 
surface.  In  both  cases,  there  is  a  large  discrepancy  between  the  model  predicted  and 
measured  heat  transfer  on  the  suction  surface  at  approximately  s/c=0.6  and  greater.  The 
boundary  layer  is  transitional  and  then  fully  turbulent  above  this  value.  The  large  scale 
eddies  in  the  flow  may  interact  with  the  larger  scale  turbulent  eddies  in  the  turbulent 
boundary  layer,  reducing  their  effectiveness  at  increasing  heat  transfer.  The  model  does 
not  account  for  any  interaction  of  scales  and  would  thus  over  predict  the  measured  heat 
transfer.  Similarly,  on  the  pressure  surface  for  the  20%  turbulence  case,  the  measured 
heat  transfer  indicates  possible  boundary  layer  transition  at  s/c=-0.5  and  above  and  the 
model  predictions  are  lower  than  the  measured  heat  transfer  and  does  not  predict  the 
turbulence  induced  transition.  In  both  of  these  regions,  the  results  demonstrate  the 
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inability  of  the  model  to  predict  early  transition  due  to  turbulence  and  over  prediction  of 
heat  transfer  in  the  fully  turbulent  boundary  layer  region. 

The  heat  transfer  predictions  with  turbulence  are  compared  to  the  measured 
values  normalized  by  the  low  turbulence  heat  transfer  (St/Sto)  in  Figure  3.18.  The  results 
in  the  figure  clearly  demonstrate  that  the  model  is  unable  to  predict  early  transition  due  to 
turbulence,  since  the  measured  St/Sto  results  show  a  large  peak  where  heat  transfer 
increases  are  quite  large.  The  figure  also  shows  that  the  results  predicted  by  the  model 
for  10%  and  20%  turbulence  are  not  greatly  different  from  one  another  and  are  both 
inside  the  “envelope”  of  the  measured  heat  transfer  for  the  two  turbulence  conditions 
over  most  of  the  blade  surface.  This  further  illustrates  the  under  prediction  of  the  20% 
turbulence  measurements  and  the  over  prediction  of  the  10%  turbulence  measurements. 
The  relatively  small  difference  in  the  modeled  heat  transfer  increases  is  a  result  of  the 
change  in  the  value  of  u’rms  being  offset  by  the  slightly  larger  length  scale  for  the  20% 
turbulence  case,  since  Equation  3.28  which  leads  to: 


(  / 


^l/2 


;rQr- 


(3.29) 


The  model  predictions  are  sensitive  to  both  u'rms  and  Ax  as  evidenced  by  Equation 
3.29.  Overall,  the  model  is  shown  to  adequately  predict  increases  in  heat  transfer,  with 
some  limitations. 

3.4.2.5  Low  speed  turbine  vane  pressure  surface  data  -  Ames,  1994 

Work  by  Ames  (1994)  also  measured  increases  in  heat  transfer  due  to  turbulence 
on  a  turbine  vane  geometry.  The  work  was  performed  in  a  low  speed  wind  tunnel  with 
several  different  turbulence  generators.  Turbulence  is  generated  using  a  combustor 
simulator  and  using  a  passive  biplanar  grid.  The  combustor  simulator  was  located  at  two 
different  distances  upstream  of  the  vane  cascade  to  give  different  turbulence  levels.  The 
turbulence  conditions  are  referred  to  by  their  generator  configuration,  where  the  grid¬ 
generated  turbulence  (Grid)  has  an  intensity  of  7.5%  and  length  scale  of  1.1  cm  at  the 
cascade  inlet,  the  combustor  simulator  closest  to  the  vane  passages  (Comb  1)  generates 
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inlet  turbulence  with  an  intensity  of  approximately  13%  and  1.6  cm  length  scale  and  the 
combustor  simulator  located  further  from  the  cascade  (Comb  2)  generates  inlet 
turbulence  with  8.5%  intensity  and  2.1  cm  length  scale. 

Measurements  of  local  turbulent  fluctuations  and  length  scale  were  performed 
using  a  hot-wire  probe  traversed  perpendicular  to  the  turbine  vane  at  four  locations  on  the 
pressure  surface.  Only  three  of  these  locations  were  used  in  the  model  application,  since 
the  fourth  location  is  at  the  trailing  edge  of  the  vane  and  heat  transfer  measurements  were 
reported  by  Ames  to  have  large  uncertainties.  The  data  presented  in  Ames  (1994) 
reported  measured  turbulent  fluctuations  and  length  scale  at  varying  distances  from  the 
surface  of  the  blade.  The  data  for  the  point  corresponding  to  the  edge  of  the  boundary 
layer  was  used  for  the  model  input.  The  model  calculated  values  of  Aht  were  determined 
and  added  to  the  measured  heat  transfer  under  low  turbulence  conditions.  Figure  3.19 
and  Figure  3.20  present  the  results  of  the  model  predictions  compared  to  the  measured 
values  with  ±4.5%  uncertainty  on  the  modeled  data  for  Stanton  number  and  ±2.4% 
uncertainty  on  the  St/Sto  data. 


Figure  3.19:  Comparison  of  model  predicted  and  measured  Stanton  numbers  (St) 

from  Ames  (1994) 
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Figure  3.20:  Comparison  of  model  predicted  and  measured  St/Sto  from  Ames 

(1994) 


Figure  3.19  indicates  that  the  model  is  fairly  good  at  predicting  the  heat  transfer 
due  to  freestream  turbulence.  With  so  few  data  points  to  compare  with  (as  opposed  to  the 
comparisons  performed  with  the  work  of  Radomsky  and  Thole  (1998)),  the  results  do  not 
give  a  strong  indication  as  to  the  ability  of  the  model  to  accurately  predict  heat  transfer 
increases.  There  does  appear  to  be  a  fairly  large  discrepancy  in  the  model  results  and  the 
measured  data  at  the  location  furthest  along  the  pressure  surface.  As  indicated  by  Ames 
(1994),  locations  far  along  the  blade  surface  (near  the  trailing  edge)  have  larger 
measurement  uncertainties. 

Figure  3.20  plots  the  model  results  against  measured  data  for  the  heat  transfer 
with  freestream  turbulence  normalized  by  the  low  turbulence  value  (St/Sto).  The  large 
difference  in  the  furthest  streamwise  measurement  location  is  well  illustrated  in  this  case. 
For  the  other  two  measurement  locations,  the  model  is  shown  to  over  predict  the 
turbulent  heat  transfer.  In  general,  the  model  does  not  appear  to  well  predict  the  Comb  2 
case.  Based  on  reported  uncertainty  in  the  measured  data,  error  bars  were  generated  for 
the  measured  heat  transfer  data  for  comparison  with  the  model  data  using  the  perturbation 
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technique  described  in  Moffat  (1988).  Figure  3.21  plots  the  data  from  Comb  2  with 
several  error  bars  on  the  curve  along  with  the  model  predicted  data.  As  shown,  the  data 
fall  within  the  uncertainty  in  the  measurement  of  the  value  of  the  normalized  heat  transfer 
data  (St/Sto). 


Figure  3.21:  Comparison  of  model  predicted  and  measured  St/Sto  from  Ames 
(1994)  for  Comb  2  with  error  estimates 

3.4.2.6  Turbulent  flat  plate  data  -  Moss  and  Oldfield,  1992 

Measurements  of  the  effect  of  freestream  turbulence  and  length  scale  on  heat 
transfer  to  a  flat  plate  with  a  turbulent  boundary  layer  were  reported  by  Moss  and 
Oldfield  (1992).  Turbulence  was  generated  by  a  series  of  upstream  turbulence  grids  of 
varying  bar  diameter  and  pitch.  Heat  transfer  was  measured  at  three  streamwise  locations 
on  the  surface  of  the  plate  with  and  without  freestream  turbulence.  Fluctuating  velocity 
was  measured  with  a  hot-wire  probe  near  the  edge  of  the  boundary  layer  at  the  surface  of 
the  plate.  The  flat  plate  boundary  layer  was  tripped  to  create  a  turbulent  boundary  layer 
at  downstream  locations.  The  heat  transfer  due  to  freestream  turbulence  was  modeled, 
similar  to  previous  cases,  by  adding  the  value  of  Aht  to  the  low  freestream  turbulence  heat 
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transfer  measurements  reported  in  the  work.  For  this  comparison,  the  grids  with  the 
highest  intensity  were  used  (Grids  D  and  E  in  Moss  and  Oldfield  (1992))  for  the 
comparisons.  The  results  of  the  modeling  are  presented  in  Table  3.8.  Uncertainty  in  the 
modeled  values  of  Nu/Nuo  was  estimated  through  the  perturbation  method  (Appendix  B) 
to  be  ±2.4%  (±0.03  Nu/Nuo). 

Table  3.8:  Model  results  compared  to  data  from  Moss  and  Oldfield  (1992) 


GridD 


Dist.  X 

Tu 

U^rins 

n 

Aht 

ho 

Calc. 

Reported 

(mm) 

(%) 

(m/s) 

(W/m^K) 

(W/m^K) 

Nu/Nuo 

Nu/Nuo 

100 

16.0 

11.2 

1.46 

86.64 

210.4 

1.41 

1.30 

180 

15.0 

10.5 

1.53 

81.95 

182.6 

1.45 

1.39 

300 

11.7 

8.2 

1.66 

69.48 

157.8 

1.44 

1.41 

GridE 


Dist.  X 

Tu 

U'rms 

A, 

Aht 

ho 

Calc. 

Reported 

(mm) 

(%) 

(m/s) 

(cm) 

(W/m^K) 

(W/m^K) 

Nu/Nuo 

Nu/Nuo 

100 

12.1 

8.5 

0.85 

98.75 

210.4 

1.47 

1.30 

180 

9.4 

6.6 

0.97 

81.47 

182.6 

1.45 

1.30 

300 

6.9 

.  4.8 

1.21 

62.50 

157.8 

1.40 

1.30 

For  grid  D,  the  model  predicts  the  measured  increases  in  heat  transfer  fairly  well, 
especially  at  the  locations  farthest  downstream.  The  first  column  is  the  streamwise 
distance  along  the  plate  measured  from  the  leading  edge.  The  model  slightly  over 
predicts  the  measured  heat  transfer.  The  model  does  not  predict  the  heat  transfer  for  grid 
E  as  well,  over  predicting  the  measured  heat  transfer  at  ail  three  streamwise  locations. 
This  is  similar  to  the  results  reported  earlier  in  the  comparison  with  the  work  of 
Radomsky  and  Thole  (1998)  where  heat  transfer  in  the  regions  of  the  turbine  vane  where 
the  boundary  layer  was  turbulent  was  over  predicted  by  the  model.  As  mentioned  earlier, 
this  may  be  due  to  the  large  scale  eddies  in  the  boundary  layer  interacting  with  the  eddies 
fi-om  the  freestream,  reducing  the  effectiveness  of  the  large  scale  eddies  in  increasing 
surface  heat  transfer. 
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Of  interest  is  the  apparent  correlation  of  the  error  in  the  predicted  measurements 
with  length  scale.  The  over  prediction  of  measured  heat  transfer  increases  as  the  length 
scale  decreases,  even  though  the  turbulence  intensity  is  decreasing.  This  suggests  that 
large-scale  turbulence  is  well  predicted  in  its  effects  on  turbulent  heat  transfer.  This  may 
indicate  that  the  primary  assumption  of  the  model  that  large-scale  eddies  penetrate  to  the 
surface  to  affect  heat  transfer  does  not  hold  as  well  for  smaller  scale  turbulence.  In  later 
work  by  the  authors.  Moss  and  Oldfield  (1994),  the  convective  velocity  of  eddy 
structures  in  the  boundary  layer  was  measured  through  spatial  correlation  of  high 
frequency  heat  flux  gages  at  different  locations.  It  was  shown  that  the  convective 
velocity  normalized  by  the  freestream  velocity  was  near  1 ,0  for  grid  D,  whereas  the  value 
was  slightly  lower  (around  0.9)  for  Grid  E,  which  has  smaller  length  scale.  These  results 
lend  credence  to  the  possibility  that  smaller  scale  freestream  turbulence  does  not 
completely  penetrate  the  boundary  layer,  and  thus  does  not  affect  surface  heat  transfer  to 
the  extent  that  the  model  would  predict,  especially  for  a  turbulent  boundary  layer. 

The  authors  (Moss  and  Oldfield  (1994))  state  that,  “The  main  heat  transfer 
enhancing  mechanism  is  due  to  the  penetration  of  freestream  turbulent  eddies  deep  into 
the  boundary  layer,  rather  than  enhancement  of  existing  boundary  layer  turbulence.” 
From  this  statement  and  the  results  of  the  analytical  model  application,  it  appears  that  this 
holds  true  when  the  eddies  are  large  enough  and  have  sufficient  energy  (intensity)  to 
penetrate  the  boundary  layer.  The  current  results  (including  the  comparison  with  the 
work  of  Radomsky  and  Thole  (1998))  suggest  that  for  turbulent  boundary  layers,  the 
eddies  may  not  completely  penetrate  to  the  surface  unless  the  freestream  eddies  are  of 
sufficient  intensity  and  scale,  leading  to  the  over  prediction  of  results  for  turbulent 
boundary  layers. 

3.4.2.7  Transonic  turbine  biade  data  -  Hoimberg,  1996 

The  final  data  set  analyzed  was  data  from  Hoimberg  (1996)  in  the  same  transonic 
turbine  cascade  facility  as  this  work.  The  data  was  taken  with  a  hot-wire  co-located  with 
a  heat  flux  gage  on  both  the  pressure  and  suction  surfaces  of  the  transonic  turbine  blade. 
Two  different  turbulence  grids  were  used.  The  hot-wire  probe  was  located  approximately 
1  mm  above  die  blade  surface.  As  discussed  in  section  1 .3,  the  work  of  Hoimberg  was  in 
a  turbulence  flowfield  with  much  lower  turbulence  intensity  (less  than  5%)  and  smaller 
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length  scale  (1-5  mm).  Table  3.9  presents  the  results  of  the  model  applied  to  the  data  for 
both  grids,  with  measurements  at  four  locations  on  the  blade  surface.  The  locations  on 
the  pressure  surface  were  performed  with  insert  heat  flux  gages  (similar  to  the  ones  used 
in  the  current  work),  while  the  suction  surface  gage  was  directly  sputtered  on  the  blade 
surface.  For  the  results  reported  for  gage  DSl,  the  grids  were  moved  to  different 
locations  to  give  different  turbulence  levels  (G#u  and  G#b). 


Table  3.9.  Model  results  compared  data  from  Holmberg  (1996) 


Location 

Grid 

Tu 

(%) 

u'rms 

(m/s) 

Aht 

(W/m^K) 

h„ 

(W/m^K) 

Predicted 

Turbulent 

h 

Measured 

Turbulent 

h 

DSl 

G2b 

2.41 

4.09 

2.37 

179 

1000 

1179 

1580 

G2u 

3.15 

5.33 

1.79 

235 

1000 

1235 

1785 

G4b 

1.84 

3.36 

119 

1000 

1119 

1360 

G4u 

2.72 

4.91 

159 

1000 

1159 

1507 

PSl 

G2 

4.57 

4.74 

1.52 

247 

800 

1047 

1085 

G4 

7.15 

Mi 

3.92 

186 

800 

986 

970 

PS2 

G2 

4.66 

4.29 

2.64 

183 

820 

1003 

930 

G4 

5.61 

5.10 

3.71 

167 

820 

987 

875 

PS3 

G2 

2.35 

2.96 

1.62 

188 

845 

1033 

1005 

G4 

3.07 

3.91 

2.50 

173 

845 

1018 

925 

The  results  demonstrate  that  the  single-scale  model  predicts  the  increase  in  heat 
transfer  reasonably  well  for  the  pressure  surface  locations,  while  the  predicted  results  for 
the  suction  surface  are  considerably  lower  than  the  measured  values.  This  is  possibly  due 
to  errors  in  the  measurement  of  heat  transfer  with  the  direct  sputtered  gage.  Reported 
uncertainties  in  the  data  were  quite  high,  with  error  as  high  as  10-15%  on  both  surfaces. 
The  luicertainty  in  the  modeled  values  of  Aht  was  estimated  to  be  4.5%. 

3.4.2.8  Film-cooled  Transonic  turbine  blade  data  -  Smith,  et  al.,  2000 

In  addition  to  applying  the  analytical  model  to  data  sets  that  measured  heat 
transfer  coefficients  with  and  without  turbulence,  the  model  was  applied  to  film-cooled 
heat  transfer  data  from  the  work  of  several  researchers  at  Virginia  Tech  on  the  same 
geometry  as  the  blades  in  the  current  work.  The  work  of  Smith,  et  al.  (2000),  Smith 
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(1999),  Popp  (1999)  and  Bubb  (1999)  investigated  the  effects  of  film  cooling  on  heat 
transfer  coefficients  and  film  cooling  effectiveness  on  the  suction  surface  of  a  transonic 
turbine  blade.  Measurements  were  performed  at  6  locations  on  the  suction  surface  of  the 
blade.  Measurements  of  time-resolved  velocity  were  only  performed  at  one  gage  location 
on  the  suction  surface  (location  SSI  in  this  work,  which  is  location  SS2  in  the  previous 
work).  With  high  ffeestream  turbulence,  the  increase  in  heat  transfer  coefficient  (Aht)  for 
this  gage  location  based  on  local  measurements  of  turbulence  intensity  and  length  scale 
was  determined  to  be  an  average  of  59  W/m^K  (see  Table  3.5).  In  order  to  apply  the 
model  to  estimate  the  heat  transfer  coefficient  with  film-cooling  and  high  ffeestream 
turbulence,  this  value  was  applied  to  each  of  the  gage  locations  from  the  previous  work. 
This  assumes  the  turbulence  intensity  and  scale  is  constant  in  the  streamwise  direction  in 
the  range  of  these  gage  locations.  This  may  not  be  the  case,  but  it  is  the  only  available 
measurement  from  the  current  work.  One  of  the  primary  model  assumptions  is  that  the 
state  of  the  boundary  layer  is  not  a  factor  in  the  analytical  model,  thus  the  model  would 
be  applicable  to  a  film-cooling  layer  as  well  as  normal  boundary  layer  flow. 
Measurements  with  film  cooling  were  not  performed  in  this  work,  but  are  recommended 
for  future  study  in  this  facility.  This  analysis  provides  an  estimate  of  what  may  be 
expected. 

Using  the  measured  film-cooled  heat  transfer  coefficients  (ho,c)  (Bubb^  1999),  the 
predicted  heat  transfer  coefficient  with  film-cooling  and  high  freestream  turbulence  was 
determined  by  adding  the  value  of  Aht  to  the  measured  film-cooled  heat  transfer 
coefficient.  Table  3.10  below  reports  the  predicted  results.  As  a  comparison  with  this 
data,  if  the  low,  turbulence  no  film-cooling  data  on  the  suction  surface  (Table  3.1  and 
Smith,  et  al,,  2000)  is  compared  to  the  high  turbulence,  suction  surface  film  cooling  data 
in  Table  3.10,  the  result  would  be  a  ratio  (ht,c/ho)  of  approximately  1.16-1.34.  Johnston, 
et  al.  (1999)  reported  data  of  this  type  for  heat  transfer  data  with  film  cooling  and  high 
freestream  turbulence  on  a  quarter  cylinder  model.  They  reported  ratios  as  high  as  1 .6  on 
areas  with  curvature,  which  is  significantly  higher  than  the  predicted  data  in  the  current 
work.  Factors  such  as  the  film  cooling  blowing  rate,  freestream  turbulence  level  and 
curvature  will  have  a  large  effect  on  the  results.  This  is  evidenced  by  the  ratios  of  film 
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cooled  heat  transfer  coefficient  and  baseline  heat  transfer  coefficient  (ho,c/ho),  which  in 
the  work  of  Johnston,  et  al.  (1999)  are  between  1.3  and  1.5  in  areas  of  curvature,  where  in 
Smith,  et  al.  (2000)  this  ratio  is  between  1 .06  and  1 .26. 


Table  3.10.  Analytical  model  applied  to  film-cooled  data  from  Smith,  et  al.  (2000) 


Gage 

X/Cp, 

Ah| 

(WWK) 

ho,c 

(W/m^K) 

Predicted 

Turbulent 

ht,c 

(W/m^K) 

ht,c/ho,C 

ht,c/ho 

SSI 

0.16 

59 

864 

923 

1.07 

1.34 

SS2 

0.20 

59 

785 

844 

1.08 

1.20 

SS3 

0.24 

59 

622 

681 

1.10 

1.16 

SS4 

0.28 

59 

783 

842 

1.08 

1.16 

SS5 

0.32 

59 

806 

1.08 

1.16 

SS6 

0.36 

59 

mm 

934 

1.07 

1.27 

3.4.3  Comparison  with  Other  Correlations 

The  existing  correlations  in  the  open  literature  typically  correlate  normalized  heat 
transfer  coefficients  (Nu,  St)  with  integral  parameters  (Tu,  Ax).  Several  of  the  existing 
correlations,  discussed  in  section  1.2.1,  were  applied  to  the  data  in  the  current  work. 
Many  of  these  correlations  require  information  on  boundary  layer  and  momentum 
thickness,  values  that  have  not  been  determined  for  the  current  work.  This  limits  the 
correlations  that  could  be  applied  for  comparison.  The  following  sections  present  the 
results  of  the  application  of  several  correlations  to  the  current  data. 

3.4.3.1  Dullenkopf  and  Mayle  (1995)  Effective  Turbulence  Correlation 

Dullenkopf  and  Mayle  (1995)  proposed  a  correlation  for  the  effects  of  freestream 
turbulence  on  laminar  heat  transfer  downstream  of  a  cylinder  stagnation  point  and  on  the 
pressure  surface  of  some  airfoils  (assuming  regions  of  constant  freestream  strain  rate). 
The  correlation  uses  an  “effective”  turbulence  level  and  takes  into  account  the  strain  rate. 
The  Nusselt  number  was  presented  as  a  function  of  this  effective  turbulence  level  by: 

Nu„  •  Pr^”  =  0.571  +  0.01  •  Tu;^  (3.30) 
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where  the  value  of  Tux  is  determine  by: 


a=strainrate  =  ■ 


dx 


..  h-^lvla 
Nu„  = - 


Tu„  = 


-Ja-v 


L  = 


A, 


■JWa 


Tu.K 


u. 


~  T  ^  1  |5/12  “  / 

1  +  0.004-47  vat' 


(3.31) 

(3.32) 

(3.33) 

(3.34) 

(3.35) 


The  strain  rate  is  used  to  form  both  a  non-dimensional  length  scale,  L®,  and  to 
form  Tux.  The  u'eir  represents  an  “effective”  turbulence  level  considering  strain  rate  and 
turbulent  length  scale.  The  correlation  was  applied  to  the  data  for  both  the  pressure  and 
suction  surfaces  of  the  transonic  turbine  blades.  The  results  are  presented  below  in 
Figure  3.22.  The  data  reflects  uncertainty  in  the  measured  heat  transfer  coefficient  of 
±6.3%  and  ±5.5%  on  turbulence  intensity.  The  data  from  the  pressure  surface  (PSl  and 
PS2)  are  in  fairly  good  agreement  with  the  Dullenkopf  and  Mayle  (1995)  correlation, 
while  the  data  for  the  suction  surface  deviates  from  the  correlation.  This  is  most  likely 
due  to  the  fact  that  the  suction  surface  location  is  in  a  region  of  highly  accelerating  flow, 
where  the  freestream  strain  rate  is  changing  dramatically,  while  the  pressure  surface 
locations  are  in  a  region  of  relatively  constant  strain  rate  (as  evidenced  by  the  Mach 
number  distribution  shown  earlier  in  Figure  2.4) 
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2.5 


Figure  3.22  Comparison  with  Dullenkopf  and  Mayle  (1995)  effective  turbulence 

correlation 

3.4.3.2  Maciejewski  and  Moffat  (1992)  St'  Correlation 

As  discussed  in  section  1.2. 1.1,  Maciejewski  and  Moffat  (1992)  presented  a 
correlation  for  the  effects  of  high  freestream  turbulence  on  flat  plate  turbulent  boundary 
layer  heat  transfer  data.  They  proposed  a  simple  method  of  analyzing  high  turbulence 
heat  transfer  data  and  observed  a  linear  relationship  between  measured  heat  transfer 
coefficient  (h)  and  turbulent  fluctuations  in  the  freestream  (u').  Their  correlation  was 
independent  of  geometry  or  length  scale,  and  was  simply  a  function  of  near  wall 
turbulence  level: 

- ; —  =  constant  (3.36) 

O-C'U 

r'  ^  *^rms 

where  ht  is  the  measured  turbulent  heat  transfer  coefficient.  In  their  paper,  the  authors 
present  a  plot  that  shows  the  quantity  /i,//?-c-M^plotted  against  turbulence  intensity  in 

a  range  up  to  approximately  60%  turbulence.  The  value  of  A,  //?  •  c  •  is  equivalent  to 
the  parameter  St',  where: 
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St'= - ^L_  =  [o.24  +  0.12  -exp- [(rw  -0.1  l)/0.055]']-/(Pr)  (3.37) 

P'C'‘*irms 

The  correlation  was  applied  to  the  current  data  and  the  results  are  plotted  in 
Figure  3.23.  The  figure  plots  the  current  data  along  with  the  correlation  with  ±  15% 
uncertainty  curves  and  uncertainty  of  ±6.3%  on  the  measured  heat  transfer  coefficient 
and  ±5.5%  on  turbulence  intensity 


Figure  3.23  Comparison  with  Maciejewski  and  Moffat  St'  correlation 


From  the  figure,  it  is  demonstrated  that  the  correlation  matches  the  measured  data 
well  for  the  pressure  surface  locations,  but  the  suction  surface  data  St’  results  are  much 
higher  than  the  correlation  predicts.  This  is  similar  to  the  result  observed  for  the 
Dullenkopf  and  Mayle  (1995)  correlation. 
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3.4.3. 3  Van  Fossen,  Simoneau  and  Ching  (1995)  Frossling  Number 
Correlation 

Van  Fossen,  et.  al  (1995)  studied  the  influence  of  turbulence  on  heat  transfer  on 
cylinders  and  cylinders  with  wedge  shaped  afterbodies.  A  correlation  was  developed 
which  included  the  integral  length  scale  in  the  form  of  the  ratio  of  the  length  scale  to 
cylinder  leading  edge  diameter  (Ax/d).  -  The  correlation  is  independent  of  body  shape 
(other  than  cylinder  diameter)  and  presented  a  correlation  in  terms  of  the  Frossling 
number  augmentation  as: 

MO),^  =  0.0085  lJraRe/‘f%l  +1  (3.38) 

where  d  is  the  leading  edge  diameter  and  Re  and  Nu  (Fr=Nu/VRe)  are  based  on  leading 
edge  diameter  and  inlet  velocity. 

Hie  correlation  was  applied  to  the  current  data  using  the  leading  edge  diameter 
(approximately  0.8  cm)  and  inlet  velocity  (140  m/s),  turbulence  intensity  (10%)  and 
length  scale  (2  cm).  Unlike  the  other  correlations  that  used  local  turbulence 
measurements,  the  correlation  of  Van  Fossen,  et  al.  (1995)  uses  turbulence  parameters 
upstream  of  the  body  (turbine  blade  in  this  case).  The  results  are  shown  in  Figure  3.24 
below,  including  error  curves  of  ±  10%  on  the  correlation  and  ±2.4%  on  the  ratio  of  the 
Frossling  number  for  the  measured  data.  The  correlation  fits  the  data  for  the  pressure 
surface  locations  extremely  well.  The  data  for  both  PSl  and  PS2  fall  on  top  of  each  other 
on  the  curve.  Similar  to  the  results  from  other  correlations  presented  previously,  the 
suction  surface  data  (SSI)  does  not  match  the  correlation  as  well,  however,  the  data  do 
fall  within  the  ±  10%  uncertainty. 
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Figure  3.24  Comparison  with  Van  Fossen,  et  al.  (1995)  Frossling  number 

correlation 


3.4.4  Multiple-scale  Model 

The  single-scale  model  results  presented  earlier  demonstrated  that  the  analytical 
model  developed  shows  promise  in  predicting  increases  in  heat  transfer  due  to 
turbulence.  The  single-scale  model  assumes  that  all  of  the  turbulent  energy  in  the  flow  is 
packed  into  one  representative  sized  eddy  (the  integral  scale).  This  is  an  idealization  of 
the  actual  flow  physics,  since  in  actuality  each  eddy  size  has  energy  associated  with  it  as 
demonstrated  by  the  power  spectrum  of  the  fluctuating  velocity.  Using  the  time-resolved 
signal  of  fluctuating  velocity,  it  is  possible  to  determine  the  energy  content  in  the 
turbulent  flowfield  for  each  individual  eddy  size  (i.e.  at  each  frequency).  Using  this 
information,  it  may  be  possible  to  predict  the  increase  in  heat  transfer  that  is  associated 
with  each  frequency  or  eddy  size  of  freestream  velocity  fluctuations.  The  model  analysis 
in  this  form  is  done  in  the  frequency  domain  and  is  referred  to  in  this  text  as  the 
“multiple-scale”  model.  The  following  sections  develop  the  multiple-scale  model  and 
analyze  its  applicability  and  present  some  key  results.  The  results  of  the  multiple-scale 
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model  analysis  demonstrate  the  frequencies  of  turbulent  eddies  which  have  the  largest 
effect  on  heat  transfer.  The  application  of  the  multiple-scale  model  to  the  current  data  set 
is  analyzed  in  various  forms  and  the  results  demonstrate  that  the  model,  coupled  with  the 
coherence  function  between  unsteady  velocity  and  heat  flux  gives  reasonable  predictions 
of  the  measured  heat  transfer  increase,  although  the  model  results  are  3  times  higher  than 
the  measured  values.  The  analysis  presents  some  interesting  points  and  attempts  to 
determine  the  range  of  frequencies,  centered  on  the  integral  length  scale,  which  will  give 
accurate  predictions  of  the  measured  increases. 

3.4.4.1  Development  of  Multiple-scale  Model 

The  multiple-scale  model  basis  is  similar  to  the  single-scale  model  with  the 
exception  that  individual  values  of  the  fluctuating  energy  and  length  scale  at  each 
frequency  are  used  as  input  into  the  model  to  come  up  with  a  “spectrum”  of  increase  in 
heat  transfer  versus  frequency  (Aqt(f)).  The  model  equation  in  the  multiple-scale  form 
presented  as: 


k^AT 


|I/2 


na 


M£) 

«'(/) 


(3.39) 


In  order  to  apply  this  equation  in  the  frequency  domain,  expressions  for  intensity  and 
length  scale  of  turbulent  fluctuations  as  a  function  of  frequency  are  needed.  From  section 
3. 1.2.4,  the  expression  for  the  integral  length  scale  was  used: 


A.(/)  =  2±  (3.40) 

For  the  turbulent  energy  at  each  frequency,  the  magnitude  of  the  fluctuating  component 
of  velocity  was  determined  in  the  frequency  domain.  The  time-resolved  fluctuating 
velocity  signal  (u'(t))  can  be  transformed  into  the  frequency  domain  by  performing  the 
Fourier  Transform  (FFT)  of  the  time  signal: 
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Uf  =  FFT{u) 


(3.41) 


The  FFT  was  performed  using  the  FFT  function  in  Matlab  software  code.  The  results  of 
the  FFT  give  the  real  and  imaginary  components  of  the  fluctuating  component  of  velocity 
at  each  frequency.  The  magnitude  of  Uf  was  then  determined  by: 


=  Vl™g(W/))“  +  Real(W))’  (3-42) 

This  magnitude  is  related  to  the  PSD  of  fluctuating  velocity  presented  in  section  3.1. 2.2 
by: 

PSD  (u)  =  (3.43) 

¥ 

since  Imag(C^)^  +  Real(t^)^  =  FFT(u)  •  FFT\u)  .  Thus,  if  we  want  to  determine 
fi'om  the  magnitude  of  u\n  the  frequency  domain  (Ufmag)»  the  equation  is  then: 


,0.5 


Z2-14U(/)’ 


(3.44) 


Equation  3.40  was  then  used  to  determine  Ax(f)  and  Equation  3.42  was  used  to  detemine 
the  input  of  u'(f)  (where  the  magnitude  of  Uf(f)  is  used).  These  values  could  then  be 
applied  in  the  model  to  yield  a  spectrum  of  increase  in  heat  flux  (as  presented  in  Equation 
3.39  and  in  Equation  3.45  below): 


A9)(/)  =  F 
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(3.45) 


where  A,(/)  = 


3.U 

Sti/ 


and  u\f)  =  Uf^{f) 
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3.4.4.2  Application  of  MultipiB-scale  Model 

Equation  3.45  was  applied  to  the  data  for  the  pressure  surface  location  nearest  the 
leading  edge  (PSl)  in  the  current  data  set  in  order  to  test  and  continue  development  of  the 
model.  Locations  PS2  and  SSI  will  be  discussed  later.  The  intensity  and  length  scale  of 
turbulent  fluctuations  were  determined  as  a  function  of  frequency  as  input  into  the 
multiple-scale  model  and  the  resultant  increase  in  heat  flux  was  plotted  as  a  function  of 
frequency.  The  model  predicted  data  are  shown  in  Figure  3.25  below. 


10^ 


Frequency  (kHz) 


Figure  3.25:  Model  predicted  heat  flux  increase  as  a  function  of  frequency 

From  the  figure,  it  is  apparent  that  the  multiple-scale  model  grossly  over  predicts 
the  heat  transfer  at  high  frequency.  The  spectrum  of  heat  flux  in  Figure  3.13 
demonstrated  that  the  energy  in  the  unsteady  heat  flux  signal  is  attenuated  at  higher 
frequencies,  thus  it  is  expected  that  the  heat  flux  increases  at  higher  frequencies  would 
also  be  attenuated.  The  breakdown  in  the  model  may  be  a  direct  result  of  one  of  the 
primary  assumptions  of  the  model,  that  the  scale  of  the  turbulence  is  much  larger  than  the 
boundary  layer  thickness.  For  a  boundary  layer  thickness  of  approximately  1  mm,  this 
would  correspond  to  eddy  scales  with  a  frequency  of  approximately  10-15  kHz.  If  we 
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assume  “much  larger”  to  be  approximately  2-3  times  the  boundaiy  layer  thickness,  this 
corresponds  to  eddy  scales  of  5-10  kHz.  The  model  would  then  break  down  above  these 
frequencies.  Mayle,  et  al.  (1998)  state  that  the  frequency  of  eddies  most  effective  at 
increasing  laminar  heat  transfer  is  about  .  For  the  data  in  this  work,  this 

frequency  is  approximately  1  kHz.  The  coherence  presented  in  Figure  3.15  earlier 
demonstrated  high  coherence  in  this  range  for  all  three  gage  locations,  and  the  frequency 
of  the  integral  length  scale  ranged  from  approximately  550-700  Hz. 

Figure  3.26  plots  the  model  predicted  heat  flux  increases  versus  the  measured  heat 
flux  fluctuations  (Qfmag)  for  comparison.  The  figure  again  demonstrates  that  the  model 
predicts  much  more  energy  (higher  increases  in  heat  flux)  at  higher  frequency  than  the 
measured  heat  flux.  The  measured  fluctuations  (Qfmag)  are  of  similar  shape  as  the 


Figure  3.26:  Comparison  of  model  predicted  heat  flux  increases  with 
measured  fluctuations 

measured  heat  flux  PSD  (Figure  3.13),  as  expected  since  the  two  are  related  (similar  to 
the  velocity  fluctuations)  by: 
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Two  key  points  arise  from  these  results.  First,  what  can  be  done  (or  more 
appropriately  what  is  missing  from  the  model)  to  qualitatively  predict  the  correct  trend  in 
the  heat  flux  increases  in  the  frequency  domain.  The  model  predicts  large  increases  in 
heat  transfer  due  to  high  frequency  velocity  fluctuations  while  the  PSD  of  heat  flux  and 
the  coherence  disagree  with  these  predictions.  The  PSD  of  measured  heat  flux  shows 
very  low  energy  at  high  frequency,  and  the  coherence  function  indicates  little  to  no 
correlation  between  the  fluctuations  in  velocity  and  heat  flux  at  frequencies  above  5  kHz 
for  the  current  data.  Second,  how  can  the  predicted  heat  flux  at  each  frequency  be 
combined  to  obtain  similar  results  to  the  single-scale  model  (i.e.  how  do  we  combine  the 
data  in  the  frequency  domain  to  determine  the  representative  heat  flux  event)?  The 
second  point  will  be  addressed  in  the  next  section. 

In  work  by  Dullenkopf  and  Mayle  (1995)  and  Mayle,  et  al.  (1998),  the  idea  of  an 
“effective”  turbulence  level  was  introduced.  Their  ideas  were  developed  based  on 
stability  theory  for  a  laminar  boundary  layer.  Fluctuations  at  high  frequencies  (small 
scale)  are  viscously  damped,  while  disturbances  at  low  frequency  (large  scale)  are  seen  as 
quasi-steady,  similar  to  bulk  flow  pulsations,  which  in  the  mean  do  not  affect  heat 
transfer.  Only  a  small  frequency  band  around  a  dominant  frequency  is  effective  at 
increasing  heat  transfer.  This  would  indicate  that  the  model  predicted  increases  at  high 
frequency  would  not  exist.  The  ability  of  freestream  turbulence  to  affect  surface  heat 
transfer  can  be  determined  by  observation  of  the  coherence  between  the  time-resolved 
velocity  and  heat  flux  signals.  The  coherence  for  each  of  the  three  gage  locations  in  this 
work  was  presented  in  Figure  3.15  in  section  3.3.2.2.  From  the  figure  (shown  again 
below  as  Figure  3.27)  it  is  apparent  that  the  effect  of  turbulent  fluctuations  on  surface 
heat  transfer  is  almost  completely  attenuated  at  frequencies  above  5  kHz. 

Due  to  the  resolution  of  the  data  at  low  frequency,  it  is  difficult  to  obtain  accurate 
information  from  the  spectra  and  coherence  functions  at  low  frequency.  The  coherence 
function  does,  however,  agree  with  the  ideas  presented  by  Dullenkopf  and  Mayle  where 
high  frequency  fluctuations  do  not  affect  heat  transfer  and  coherence  between  velocity 
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and  heat  flux  is  only  present  in  a  narrow  frequency  range  around  the  integral  length  scale 
(as  previously  discussed).  The  dashed  lines  on  the  coherence  plot  represent  the  range  of 
the  frequencies  of  the  integral  length  scales  measured  at  the  three  blade  locations  with  an 
uncertainty  band  of  ±  5%  applied.  The  logical  step  is  then  to  make  use  of  the  knowledge 
gained  from  the  time-resolved  measurements  of  velocity  and  heat  flux  in  the  form  of  the 
coherence  function  and  incorporate  the''coherence  into  the  multiple-scale  model.  Since 
the  coherence  is  a  dimensionless  function,  it  was  decided  to  perform  a  direct 
multiplication  of  the  coherence  function  with  the  model  output,  coupling  the  two  data  sets 
in  the  frequency  domain. 


Figure  3.27:  Coherence  between  velocity  and  heat  flux  signals 


The  resultant  revised  model  equation  is  then: 


H(f)  = 


k-AT 


na 


Kif) 

Af) 


ni/2 


riif) 


(3.47) 
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where  the  value  ylgif)  is  the  measured  coherence  function  between  the  time-resolved 

velocity  and  heat  flux  data  as  a  function  of  frequency.  This  value  is  then  equal  to  the 
coherent  output  power  (COP)  or  the  energy  resultant  in  the  output  signal  (surface  heat 
flux)  that  is  correlated  with  the  input  signal  (freestream  velocity  fluctuations).  The 
resultant  model  predicted  data  are  shown  in  Figure  3.28a. 

From  the  figure,  it  is  shown  that  the  model  predicted  data  follows  the  same 
general  trend  as  the  measured  heat  flux  spectra,  although  there  is  some  disagreement  in 
the  magnitude-  of  the  heat  flux.  This  difference  in  magnitude  is  further  illustrated  in 
Figure  3.28b,  which  presents  the  same  data  as  Figure  3.28a  plotted  on  a  linear  scale.  The 
application  of  the  coherence  function  has  effectively  attenuated  the  large  predicted 
increases  in  heat  transfer  at  higher  frequencies,  which  is  expected  from  the  coherence  and 
measured  heat  flux  spectra.  Although  the  predicted  heat  flux  at  high  frequency  is 
reduced  and  the  spectrum  of  heat  flux  increase  has  a  shape  that  is  more  consistent  with 
the  measured  spectra,  the  question  of  how  to  use  this  data  to  find  a  total  increase  in  heat 
transfer  still  remains. 


Figure  3.28a:  Predicted  model  heat  flux  with  coherence  applied  and  measured  heat 

flux  fluctuations 
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Figure  3.28b:  Predicted  model  heat  flux  with  coherence  applied  and 
measured  heat  flux  fluctuations  (Linear  scale) 

The  second  key  point  discussed  earlier  is  then  important,  how  can  the  data  be 
combined  in  the  frequency  domain  to  yield  a  total  increase  in  heat  transfer?  Is  the  total 
increase  in  heat  flux  equal  to  the  sum  of  the  predicted  increase  at  all  frequencies?  From  a 
physical  standpoint,  summation  across  all  frequencies  would  indicate  that  each  individual 
eddy  acts  independently  of  eddies  at  other  frequencies  to  increase  surface  heat  transfer. 
The  next  step  is  then  to  determine  how  to  combine  the  heat  flux  increases  for  the 
frequency  domain  modeled  data  to  give  a  value  of  predicted  heat  flux.  The  following 
section  investigates  several  different  methods  for  combining  the  data. 

3.4.4.3  Combination  of  Muitipie-scaie  Modei  Data 

In  the  single-scale  model,  the  increase  in  heat  transfer  is  partially  governed  (along 
with  the  length  scale)  by  the  root  mean  square  (RMS)  of  the  velocity  fluctuations,  rather 
than  treating  each  eddy  size  with  its  associated  energy  independently.  The  RMS  of  the 
fluctuating  energy  is  packed  into  one  representative  eddy  size  in  the  single-scale  model. 
This  concept  may  be  applied  to  the  summation  of  the  predicted  heat  flux  increases,  which 
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inherently  assumes  that  if  one  eddy  of  smaller  size  (higher  frequency)  is  resident  at  the 
surface  with  a  larger  (lower  frequency)  eddy,  the  two  do  not  act  independently  of  one 
another  to  increase  heat  transfer,  rather  the  two  eddies  combine  energy  in  a  non-linear 
fashion.  Graphically,  the  combined  increases  in  heat  flux  can  be  represented  in  the  time 
domain  by  Figure  3.29. 
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Figure  3.29:  Illustration  of  increase  in  heat  transfer  by  turbulent  fluctuations  in 

heat  flux 

The  figure  demonstrates  the  idea  that  the  fluctuations  in  heat  transfer  can  be 
added  to  the  heat  flux  without  freestream  turbulence,  assuming  that  even  eddies  moving 
away  from  the  surface  (negative  velocity  fluctuations)  cause  an  increase  in  heat  transfer. 
The  value  of  Ah,  in  the  figure  would  then  be  equal  to  the  RMS  of  the  fluctuating 
component  of  measured  heat  flux.  This  also  assumes  that  ALL  frequencies  of  heat  flux 
fluctuations  are  included  in  the  total  increase.  As  the  coherence  has  demonstrated,  only 
frequencies  in  the  range  up  to  5  kHz  are  correlated  with  velocity  fluctuations.  The 
modeled  heat  flux  could  then  be  combined  by  finding  the  total  heat  flux  increase  by  using 
an  equation  of  similar  form  as  was  used  for  determining  RMS  of  velocity.  This 
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essentially  characterizes  a  representative  heat  flux  event,  similar  to  how  the  RMS  of 
fluctuating  velocity  characterizes  turbulence. 

The  combination  of  the  heat  flux  increases  for  the  multiple-scale  model  was 
investigated  in  various  forms  in  an  attempt  to  determine  how  to  apply  the  multiple-scale 
model  and  to  investigate  which  frequencies  have  the  dominant  effect  on  surface  heat 
transfer.  The  methods  used  to  combine  the  heat  flux  increases  are  discussed  below  along 
with  an  illustration  of  the  application  of  the  equations  to  the  data  for  each  method.  The 
results  of  each.diflerent  method  are  presented  in  Table  3.10.  The  methods  can  essentially 
be  categorized  into  three  forms.  The  first  (Methods  1  and  2)  uses  the  RMS  of  the 
measured  time-resolved  heat  flux  data  to  compare  these  results  to  the  average  measured 
heat  flux  increase  and  single-scale  model  data  (as  illustrated  in  Figure  3.29),  the  second 
(Methods  3  and  4)  uses  the  multiple-scale  model  developed  in  the  previous  section  and 
the  third  (Methods  6  and  7)  is  a  variation  on  the  single-scale  model,  using  the  frequency 
domain  velocity  data  to  determine  the  ratio  of  Ufmag(f)  and  Ax(f)  from  equation  3.45  and 
finding  the  RMS  of  this  data.  This  method  is  essentially  a  hybrid  of  the  single-scale  and 
multiple-scale  models,  applying  concepts  developed  for  the  multiple-scale  model  to  the 
velocity  and  length  scale  (determining  each  in  the  frequency  domain)  then  entering  an 
RMS  value  of  the  ratio  between  these  two  parameters  into  the  single-scale  model.  For 
each  of  these  three  methods,  the  results  are  presented  with  and  without  the  coherence 
function  applied.  In  addition,  the  total  increase  in  heat  transfer  was  determined  using  the 
whole  frequency  range  as  well  as  with  the  data  combined  only  over  the  frequency  range 
of  significant  coherence  (0-5  kHz).  The  single-scale  model  is  presented  for  comparison 
as  Method  5. 
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Methods  investigated  for  combining  heat  flux  increases  in  frequency  domain: 

1.  Using  measured  Qfmag,  determine  Aqt,totai  using  RMS  equation  for  all 
frequencies  and  for  data  from  0-5  kHz: 


Across  all  frequencies: 


(3.48a) 


From  0-5  kHz  (range  of  coherence): 

,toto/(0-5kHz)  ~  I  (3.48b) 

Y  f=0-5kHz 


Real  physical 
mechanism 


Measured 


2.  Using  measured  Qfmag  multiplied  by  the  coherence  (Yuq^)>  determine  Aqt,totai 
using  RMS  equation  for  all  frequencies  and  for  data  from  0-5  kHz: 


Across  all  frequencies: 

^,,,0,01  =  JX^'(Q/mag(/)’r^(/)y  (3.49a) 
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From  0-5  kHz  (range  of  coherence): 
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3.  Using  multiple-scale  modeled  Aqt(f),  determine  Aq^otai  using  RMS  equation 
for  all  frequencies  and  for  data  from  0-5  kHz: 

Across  all  frequencies: 

2  •  ,  (3.50a) 

From  0-5  kHz  (range  of  coherence): 

^r./oto/(o-5kHz)  “  /  ^2*Agf,(/)^  (3.50b) 
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4.  Using  multiple-scale  modeled  Aqt(f)  multiplied  by  the  coherence  (Vuq^), 
determine  Aqt,totai  using  RMS  equation  for  all  frequencies  and  for  data  from  0- 
5  kHz:  . 


Across  all  frequencies: 


total 
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From  0-5  kHz  (range  of  coherence): 
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5.  Using  measured  u'nns  and  Ax  (equivalent  to  single-scale  model)  across  all 
frequencies: 


^RMS  ~ 
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6.  Using  RMS  of  u'/Ax  including  data  across  all  frequencies  and  from  0-5  kHz: 


Across  all  frequencies: 


A.(/)J 


(3.54a) 


From  0-5  kHz  (range  of  coherence): 
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7.  Using  RMS  of  uVAx  multiplied  by  the  coherence  including  data  across 
all  frequencies  and  from  0-5  kHz: 

Across  all  frequencies: 
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Using  each  of  these  methods,  the  total  increase  in  heat  transfer  due  to  turbulence  was 
determined  and  the  results  are  shown  below  in  Table  3.10. 
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Table  3.11:  Multiple-scale  model  results  by  various  methods  for  location  PSl 


Method 

A(]  total 

Aqtotii  (0-5  kHz) 

(W/m^) 

(W/m^) 

1 

2785 

2698 

2 

1668 

1109 

3 

,100,400 

32,560 

4 

6513 

6467 

5 

2124 

6 

6048 

3283 

7 

1367 

1360 

The  results  from  the  table  were  analyzed  to  determine  the  method  that  most 
closely  predicts  the  measured  heat  flux  increase  and  the  value  from  the  single-scale 
model  (Method  5).  From  the  table,  it  appears  as  if  the  idea  presented  in  Figure  3.29,  that 
the  RMS  of  the  fluctuating  velocity  can  be  added  to  the  measured  low  turbulence  heat 
transfer  to  predict  the  heat  transfer  with  turbulence,  is  a  reasonable  assumption.  The 
value  of  Aqtotai  determined  by  this  method  (Method  1)  is  very  close  to  the  measured  and 
single-scale  model  predicted  increase  in  heat  transfer  (Method  5).  Applying  the 
coherence  to  the  measured  RMS  fluctuations  (Method  2)  decreases  the  predicted  value. 
Clearly,  applying  the  coherence  to  measured  heat  flux  data  is  not  necessary,  since  it  is 
data  that  is  already  a  function  of  the  coherence  between  velocity  and  heat  flux  based  on 
the  real  physical  mechanism  by  which  turbulence  affects  surface  heat  transfer. 

The  results  of  Methods  3  and  4  were  presented  earlier  graphically  in  Figure  3.26 
and  Figure  3.28,  respectively.  Figure  3.26  presents  the  multiple-scale  model  data  without 
the  coherence  function  applied,  which,  as  discussed  earlier,  does  not  take  into  account  the 
fact  that  high  frequency  increases  in  heat  transfer  are  attenuated.  The  results  in  Table 
3.10  demonstrate  that  finding  the  value  of  Aqt,totai  for  this  data  gives  extremely  high 
predictions  of  the  increase  in  heat  transfer.  The  data  with  the  coherence  applied  (Method 
4)  predicts  increases  in  heat  transfer  on  the  order  of  the  measured  increases,  but 
approximately  3  times  higher.  This  method  makes  the  most  physical  sense  and  shows 
promise  in  predicting  heat  transfer  increases  from  frequency  domain  data,  and  may 
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possibly  be  used  to  determine  the  frequency  band  which  is  most  effective  at  increasing 
heat  transfer.  This  method  will  be  revisited  shortly. 

The  results  in  Table  3.10  for  Methods  6  and  7  demonstrate  that  both  methods 
predict  an  increase  in  heat  transfer  that  are  on  the  same  order  as  the  measured  increase. 
As  discussed  earlier,  this  method  is  essentially  a  hybrid  of  the  single-scale  and  multiple- 
scale  models.  The  method  uses  the  ihput  velocity  and  length  scale  to  determine  the 
residence  time  of  each  individual  frequency  of  fluctuation  and  uses  an  equation  of  the 
form  of  the  RMS  to  determine  a  representative  residence  time.  This  residence  time  is 
then  entered  into  the  model  equation.  From  the  values  presented  for  the  data  out  to  5  kHz 
in  the  two  methods  it  can  be  seen  that  results  of  these  methods  bound  the  measured  (and 
single-scale  model)  increase  in  heat  transfer. 

Revisiting  the  results  presented  for  Method  4,  the  frequency  range  for  which  data 
were  included  was  investigated  by  centering  the  frequency  band  on  the  frequency  of  the 
integral  length  scale.  This  frequency  was  determined  in  section  3.3 .2.2  to  be  570  Hz  for 
the  data  presented  above  (location  PSl  -  pressure  surface  nearest  the  leading  edge).  The 
multiple-scale  model  results  with  the  coherence  applied  were  combined  using  equation 
3.49b,  but  rather  than  combining  data  from  0-5  kHz,  a  frequency  band  spanning  a  range 
of  fA  +-  Afb  was  used  (where  fA  is  the  frequency  corresponding  to  the  measured  integral 
length  scale  (570  Hz)  and  Afr  is  the  half-width  of  a  frequency  band  centered  on  fA).  This 
value  of  Afb  was  varied  until  the  total  heat  flux  increase  (Aq^otai)  was  approximately 
equal  to  the  measured  increase  in  heat  transfer.  The  results  of  modeling  showed  that  Afr 
was  equal  to  approximately  75  Hz.  This  result  is  demonstrated  graphically  in  Figure  3.30 
below.  In  the  flgure,  the  modeled  heat  flux  data  is  identical  to  that  presented  in  Figure 
3.28  except  die  data  has  been  frequency  averaged  to  reduce  the  noise  in  the  spectrum  for 
presentation  purposes. 

The  analysis  of  the  model  data  in  this  form  is  similar  to  the  analysis  performed  in 
Dullenkopf  and  Mayle  (1995).  In  that  work,  it  was  stated:  “The  effective  turbulence 
level,  Tueff,  is  obtained  by  integrating  the  spectral  energy  density  over  a  band  width  of 
frequencies  centered  about  the  dominant  frequency.”  The  dominant  frequency  in  their 
case  was  a  dimensionless  frequency  ((|)x)  based  on  the  frequency  of  the  integral  length 
scale  and  the  frequency  band  was  A<t).  Their  frequency  band  is  a  dimensionless 
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equivalent  of  2Afb.  The  work  of  Dullenkopf  and  Mayle  applied  this  concept  to  the 
velocity  spectra  and  then  used  this  effective  turbulence  level  in  a  correlation  with  the 
integral  length  scale  to  determine  increases  in  heat  transfer.  Their  results  showed  that 
the  frequency  band  which  is  effective  at  increasing  surface  heat  transfer  is  much  smaller 
than  the  dominant  frequency  (A(|K<(t)x).  In  this  work,  the  heat  transfer  has  been  modeled 
and  this  concept  is  applied  to  the  resultant  data  with  similar  results.  The  frequency  band 
that  yields  the  measured  increase  in  heat  transfer  is  small  (approximately  1/10  of  the 
integral  scale  frequency),  as  illustrated  in  Figure  3.30. 


Figure  3.30:  Illustration  of  frequency  band  of  multiple-scale  heat  flux  model 

With  the  data  modeled  and  analyzed  for  one  location  in  the  current  data  set,  the 
next  logical  step  is  to  apply  the  analyses  to  the  other  two  locations  to  increase  the  number 
of  data  points  for  the  analysis  of  the  different  methods.  The  methods  presented  above 
were  applied  to  the  data  sets  for  location  PS2  (on  the  pressure  surface  of  the  blade  further 
along  the  blade  chord  in  the  streamwise  direction)  and  for  location  SSI  (on  the  suction 
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surface  of  the  turbine  blade).  The  results  for  these  two  locations  are  show  in  Table  3.1 1 
in  the  same  format  as  the  data  presented  in  Table  3.10  for  location  PSl. 


Table  3.12:  Multiple-scale  model  results  by  various  methods  for  locations  PS2  and 


SSI 


Method 

>  Aqtotsi 
(W/m^) 

Aqtoiai  (0-5  kHz) 
(W/m^) 

Location  PS2 

1 

3740 

3636 

2 

2183 

1438 

3 

53,365 

16,186 

4 

3158 

3131 

5 

1116 

6 

3218 

1638 

7 

668 

663 

Location  SSI 

1 

3861 

3685 

2 

2075 

1321 

3 

37,280 

9769 

4 

1873 

1849 

5 

668 

6 

2201 

990 

7 

398 

392 

The  results  demonstrate  several  differences  from  the  results  for  location  PSl 
(Table  3.10).  The  results  from  Method  1  for  both  PS2  and  SSI  do  not  agree  well  with  the 
single-scale  model  results  (although  they  are  of  the  right  order  of  magnitude),  while  for 
PSl  the  value  for  Method  1  was  in  good  agreement  with  the  single-scale  model.  Method 
4  demonstrates  similar  results  as  the  results  for  PSl,  where  the  multiple-scale  model 
prediction  was  approximately  3  times  higher  than  the  single-scale  model  value.  Methods 
6  and  7  also  give  similar  results  for  these  two  locations,  where  the  results  of  data 
combined  out  to  5  kHz  between  the  two  methods  bound  the  single-scale  model  value. 

The  frequency  range  for  Method  4  that  matches  the  measured  and  single-scale 
model  results  was  investigated  for  both  of  these  locations,  centering  the  range  on  the 
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frequency  corresponding  to  the  measured  integral  length  scale  for  each  location,  700  Hz 
and  650  Hz  for  PS2  and  SSI,  respectively.  The  frequency  range  Avas  approximately  30 
Hz  for  both  locations.  Ibe  coherence  for  both  of  these  locations  (shown  in  Figure  3.27) 
was  lower  than  for  location  PS  1,  and  both  PS2  and  SSI  show  more  of  a  sharp  peak  near 
the  frequency  of  the  integral  length  scale.  The  coherence,  combined  with  the  smaller 
increase  in  heat  flux  results  in  a  nairower  frequency  band  to  match  the  measured 
increase. 

3.4.4.4  Application  of  Multiple-scale  Model  to  Other  Data 

The  multiple-scale  model  was  applied  to  the  cylinder  stagnation  region  data  of 
Oilier  and  Van  Fossen  (1996),  which  allows  the  model  to  be  tested  on  a  different 
geometry  and  flowfield.  The  requirements  for  application  of  the  multiple-scale  model 
were  measurements  of  steady  heat  transfer  with  and  without  freestream  turbulence,  as 
well  as  time-resolved  velocity  and  heat  flux  data.  The  model  was  applied  to  the  data  for 
grid  42  in  Oilier  and  Van  Fossen’ s  data  for  the  hot-wire  location  closest  to  the  surface  of 
the  cylinder.  The  results  for  each  of  the  methods  presented  earlier  are  listed  in  Table  3.12 
below.  Note  that  the  data  in  the  second  column  are  included  in  a  range  up  to  2  kHz  (as 
opposed  to  5  kHz  in  the  earlier  data)  since  the  coherence  range  for  this  data  was  only  2 
kHz.  It  should  be  mentioned  that  none  of  the  data  analyzed  from  the  work  of  Oilier  and 
Van  Fossen  demonstrated  significant  levels  of  coherence,  with  most  values  measured  as 
less  than  0.2. 

Table  3.13:  Multiple-scale  model  results  by  various  methods  for  data  from  Diller 


and  Van  Fossen  (1996) 


Method 

AC}total 

Aqtotai  (0-2  kHz) 

(W/m^) 

(W/m^) 

1 

974 

927 

2 

139 

24 

3 

33,770 

14,421 

4 

281 

243 

5 

1937 

6 

5697 

1118 

7 

149 

135 
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For  this  data  set,  it  appears  that  none  of  the  various  methods  are  good  predictors 
of  die  measured  heat  transfer,  except  Method  6,  which  under  predicts  the  increase,  but  is 
relatively  close  to  the  measured  value  compared  to  the  other  methods.  The  low 
coherence  between  the  time-resolved  velocity  and  heat  flux  data,  when  applied  to  the 
model  data,  greatly  reduces  the  modeled  values. 
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Chapter  4 


Conclusions 


4.i  Summary 

The  primary  goals  of  this  work  were  to  measure  and  analytically  model  the 
increases  in  heat  transfer  due  to  the  generated  high  intensity,  large  scale  freestream 
turbulence  field  based  on  measured  time-averaged  baseline  (low  turbulence)  and  high 
turbulence  surface  heat  transfer  and  on  time-resolved  measurements  of  surface  heat  flux 
and  local  fluctuating  velocity.  Along  the  way,  it  was  hoped  that  the  model  and  its 
application  would  shed  some  light  onto  the  mechanism  by  which  freestream  velocity 
fluctuations  affect  surface  heat  transfer. 

Section  3.2.3  presented  the  low  turbulence  heat  transfer  coefficients  along  with 
the  measured  high  turbulence  heat  transfer  coefficients.  The  increase  in  heat  transfer  was 
shown  to  be  approximately  two  times  higher  on  the  pressure  surface  than  on  the  suction 
surface  (17%  on  the  pressure  surface  versus  8%  on  the  suction  surface).  This  result  is 
consistent  with  observations  of  other  researchers.  This  corresponds  to  increases  in  blade 
surface  temperature  of  5-10%,  which  can  significantly  reduce  the  life  of  a  turbine  blade 

The  analytical  model  was  presented  in  two  forms,  documented  as  the  “single¬ 
scale”  and  “multiple-scale”  models.  The  results  presented  in  sections  3.4.2  and  3.4.3 
demonstrated  how  the  single-scale  model  developed  predicted  the  measured  increases  in 
heat  transfer  for  the  data  in  this  work,  as  well  as  for  the  data  from  the  work  of  various 
other  researchers  in  varying  flowfields  and  on  different  geometries  (turbine  vane,  turbine 
blade,  cylinder  stagnation  region  and  flat  plate).  The  model  was  shown  to  fairly 
accurately  predict  the  increases  in  heat  transfer,  except  in  regions  of  turbulent  flow  and 
transition.  The  model  was  shown  to  slightly  over  predict  the  heat  transfer  for  turbulent 
boundary  layers  and  demonstrated  a  deficiency  in  predicting  early  laminar  to  turbulent 
boundaiy  layer  transition  caused  by  elevated  levels  of  freestream  turbulence. 
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The  multiple-scale  model  development  provided  some  interesting  insight  into  the 
frequencies  of  turbulent  eddies  which  are  most  effective  at  increasing  surface  heat 
transfer.  It  was  shown  that  the  time-resolved  velocity  data  could  be  transformed  into  the 
frequency  domain  and  decomposed  into  the  individual  eddy  sizes  (scales  or  frequencies), 
and  that  the  energy  associated  with  each  individual  scale  of  eddy  could  be  determined. 
The  multiple-scale  model  then  applied  die  analytical  model  developed  to  each  eddy  scale 
to  determine  a  spectrum  of  heat  flux  increase.  The  time-resolved  heat  flux  data  was  used 
to  compare  with  this  spectrum  of  increases  and  also  to  determine  the  coherence  between 
the  fluctuating  velocity  and  surface  heat  transfer.  The  coherence  and  the  measured  heat 
flux  spectra  demonstrated  that  there  is  a  small  band  of  frequencies  for  which  freestream 
velocity  fluctuations  affect  surface  heat  transfer,  roughly  centered  around  the  frequency 
of  the  integral  length  scale  (most  energetic  eddies)  of  the  freestream  turbulence.  This 
result  was  in  good  agreement  with  previous  work  by  Dullenkopf  and  Mayle  (1995)  and 
Mayle,  et  al.  (1998). 

The  multiple-scale  model  attempted  to  determine  how  to  combine  the  spectrum  of 
increases  in  heat  transfer  by  various  methods  with  moderate  success.  The  resultant 
predicted  increases  were  on  the  right  order  of  magnitude  with  the  measured  increases,  but 
were  approximately  a  factor  of  3  greater. 

What  was  learned: 

1.  The  single-scale  model,  developed  based  on  a  physical  model,  demonstrated 
good  ability,  within  certain  limitations,  at  predicting  increases  in  heat  transfer 
due  to  freestream  turbulence.  Based  on  the  results,  the  mechanism  of  the 
model  seems  reasonable. 

2.  The  current  data  set  agrees  well  with  existing  correlations,  although  there  is 
some  discrepancy  with  the  suction  surface  data. 

3.  Both  measured  heat  flux  spectra  and  coherence  between  freestream  velocity 
and  surface  heat  transfer  fluctuations  are  valuable  tools  in  analyzing  the 
frequencies  (eddy  sizes)  of  turbulence  that  are  effective  at  increasing  surface 
heat  transfer. 
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4.  The  analytical  model  assumption  of  large-scale  eddies  penetrating  the 
boundary  layer  and  superimposing  their  energy  in  the  form  of  heat  flux 
increases  on  the  mean  heat  flux  appears  to  be  reasonable.  The  frequencies  in 
the  measured  spectra  and  coherence  functions  where  energy  and  coherence  are 
attenuated  correlate  well  with  the  integral  length  scale  and  effective 
turbulence  frequencies  presented  by  Mayle,  et  al.  (1998).  Above  these 
frequencies,  the  effectiveness  of  the  turbulent  freestream  eddies  at  increasing 
heat  transfer  decreases  rapidly  due  to  viscous  dissipation  of  energy  at  higher 
frequencies  (smaller  scales). 

5.  The  heat  flux  increases  in  the  multiple-scale  model  do  not  seem  to  add  as  the 
RMS  of  the  increases  of  the  individual  eddies.  The  single-scale  model,  where 
the  increase  is  a  function  of  the  combined  turbulent  energy  at  the  frequency  of 
the  integral  scale,  was  much  more  effective  at  predicting  increases  in  surface 
heat  transfer  due  to  turbulence. 

4.2  Application  to  Turbine  Design 

The  results  of  the  time-resolved  measurements  provide  a  very  limited  data  set  that 
is  applicable  to  design  of  transonic  rotor  blades.  The  measurements  were  limited  to  three 
blade  locations  that  have  been  reported  in  previous  work  in  the  transonic  facility.  The 
results  of  these  measurements  agree  well  with  existing  correlations.  The  application  of 
the  analytical  model  developed  has  reasonable  value  to  turbine  designers,  within 
limitations.  The  model  requires  knowledge  of  the  local  turbulence  intensity  (RMS  of 
freestream  velocity  fluctuations  near  the  edge  of  the  boundary  layer)  and  integral  length 
scale.  These  values  may  not  be  known  for  real  engine  conditions,  although 
computational  fluid  dynamics  predictions  may  be  able  to  determine  a  reasonable 
distribution  of  these  values  as  a  function  of  the  inlet  turbulence  conditions,  making  the 
model  moderately  applicable  to  engine  development.  The  results  of  the  multiple-scale 
model  analysis  do,  however,  provide  insight  into  the  scales  of  turbulence  that  are  most 
effective  at  affecting  surface  heat  transfer. 
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4.3  Outlook  and  Future  Work 

The  present  investigation  generated  turbulence  that  is  consistent  with  recent 
reports  of  the  turbulent  flowfield  from  modem  gas  turbine  combustors  entering  the  first 
stage  rotor  blades.  The  turbulence  intensity  and  scale  were  not  varied  to  provide  a  more 
extensive  data  set  for  the  analyses.  Previous  work  in  this  facility  varied  intensity  and 
scale;  however,  neither  was  consistent  with  what  is  currently  believed  to  consistent  with 
real  engine  conditions.  Future  work  should  vary  the  intensity  and  scale  (within 
reasonable  limjts  compared  to  the  model  of  engine  conditions)  and  should  investigate 
more  gage  locations.  Gages  currently  being  developed  provide  a  higher  density  of 
measurement  locations  with  a  smaller  spatial  resolution.  Provided  that  the  time-response 
of  these  gages  is  sufficient  to  resolve  the  scales  of  turbulence  present  in  the  flowfield, 
more  gage  locations  can  be  sampled  simultaneously.  This  provides  more  measurement 
locations  and  higher  spatial  resolution,  but  will  also  allow  for  cross  correlation  of  the 
time  signals  from  various  gages.  Cross  correlation  in  the  time  domain  adds  the  ability  to 
track  the  progression  of  eddies  in  the  boundary  layer,  providing  an  additional  analysis 
technique. 

Within  the  measurement  limitations  of  the  facility,  the  state  of  the  boundary  layer 
and  the  details  of  the  boundary  layer  thickness  and  growth  should  be  determined.  These 
details  were  not  resolved  in  this  work,  but  with  great  care  optical  techniques  such  as 
shadowgraph  and  Schleiren  photography  and  oil  flow  visualization  should  provide  decent 
approximations  of  these  parameters. 

Finally,  similar  analyses  should  be  performed  with  a  turbulent  boundary  or  film 
layer.  Film  cooling  and  turbulent  boundary  layers  have  associated  energy  and  scales  of 
turbulence  that  can  be  measured  with  techniques  similar  to  those  presented  in  this  work. 
The  analytical  model  application  could  be  tested  under  these  conditions  to  test  its 
validity. 
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Appendix  A 


Turbulence  Generation  Paper  and  Presentation 

The  following  pages  detail  the  development  and  testing  of  the  turbulence 
generation  grid  in  the  form  of  a  paper  and  presentation.  The  paper,  written  by  the  author 
of  this  dissertation  and  submitted  to  the  2002 IGTI  conference  in  Amsterdam,  NL,  details 
the  design  and  testing  of  the  turbulence  generation  grid.  The  presentation  given  by  the 
author  at  the  conference  follows  the  paper  and  includes  additional  data  (hot-wire  surveys) 
not  included  in  the  paper. 
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ABSTRACT 

Heat  transfer  predictions  in  gas  turbine  engines  have  focused  on 
cooling  techniques  and  on  the  effects  of  various  flow  phenomena.  The 
effects  of  wakes,  passing  shock  waves  and  freestream  turbulence  have 
alt  been  of  primary  interest  to  researchers.  The  focus  of  the  work 
presented  in  this  paper  is  to  develop  a  turbulence  grid  capable  of 
generating  high  intensity,  large-scale  turbulence  for  use  in 
experimental  heat  transfer  measurements  in  a  transonic  facility.  The 
grid  is  desired  to  produce  freestream  turbulence  characteristic  of  the 
flow  exiting  the  combustor  of  advanced  gas  turbine  engines.  A 
number  of  techniques  are  discussed  in  this  paper  to  generate  high 
intensity,  large  length-scale  turbulence  for  a  transonic  facility. 
Ultimately,  the  passive  grid  design  chosen  is  capable  of  producing 
freestream  turbulence  with  intensity  of  approximately  10-12%  near  the 
entrance  of  the  cascade  passages  with  an  integral  length-scale  of  2  cm. 

NOMENCLATURE 

c  Blade  chord  length 

di  Original  grid  bar  diameter  -  2.54  cm  ( 1  in.) 

d2  Modified  grid  bar  diameter-  5.08  cm  (2  in.) 

Dw  Hot-wire  probe  diameter  (5pm) 

E(f)  Energy  spectra  of  streamwise  velocity  component 
f  Frequency 

I  Momentum  flux  ratio  (pjUjVp  ja»^) 

kro  Thermal  conductivity  of  air  at  mean  temperature 

M  Mass  flux  ratio  (pjU/pcju 

Ma  Mach  number 

Nu'  Nusselt  number  of  hot-wire 

P,  Freestream  static  pressure 

P,  Freestream  total  pressure 

R , ,  Autocorrelation  of  streamwise  velocity  component 

Re  Reynolds  number 

Re».  Reynolds  number  of  hot-wire 

T  Integral  lime  scale 

T,„  Mean  temperature  ((T».+T,)/2) 


T, 

T. 

T„ 

Freestream  static  temperature 
Freestream  total  temperature 

Hot-wire  operating  temperature 

Tu 

Turbulence  intensity  l^j 

it 

Instantaneous  velocity  {U  +  «) 

U 

f  1  7  1  A' 

Mean  velocity  ly 

u 

Fluctuating  component  of  velocity  (i7 

-u) 

u' 

f  1  A 

RMS  of  fluctuating  velocity  I  —2“, 

J-l 

V 

Hot-wire  output  voltage 

Greek  Symbols 

K  Wave  number  (27tf/U) 

A*  Integral  length -scale 

p  Density  of  air 

p„,  Density  of  air  at  mean  temperature 

T  Autocorrelation  time 

Dynamic  viscosity  of  air  at  mean  temperature 

INTRODUCTION 

The  desire  for  increased  engine  efficiency  in  modem  gas  turbine 
engines  has  led  to  increases  in  thermal  load  on  high-pressure  turbine 
blading  due  to  high  combustor  exit  temperatures.  The  flowficid  in  the 
turbine  section  is  highly  turbulent  as  a  result  of  combustion  dynamics 
and  blade  row  interactions.  Prediction  of  blade  surface  temperatures 
therefore  requires  researchers  to  study  the  effects  of  unsteady  flow 
phenomena  on  turbine  blade  surface  heat  transfer.  These  unsteady 
phenomena  include  wakes  and  shock  waves  due  to  rotor-stator  blade 
passing  and  freestream  turbulence  mainly  generated  in  the  combustor. 
The  present  research  focuses  on  the  effects  of  combustor  exit 
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freeslrcam  turbulence  on  blade  surface  heat  transfer  and  film  cooling 
effectiveness. 

In  order  to  study  the  realistic  effect  turbulence  has  on  heat 
transfer  in  an  experimentally  controlled  environment,  the  turbulence 
must  be  generated  to  model  combustor  exit  turbulence.  Several 
research  groups  have  generated  high  levels  of  freestream  turbulence 
using  different  technique  that  were  investigated  specifically  for^ their 
applicability  in  a  transonic  blow  down  cascade  wind  tunnel  facility. 
TTiese  techniques  can  be  categorized  into  mock  combustor/jets  in 
crossflow  and  passive  and  active  grids.  The  majority  of  the  previous 
research  with  high  levels  of  freestream  turbulence  has  been  in  wind 
tunnel  facilities  with  relatively  low  speed  flows.  These  facilities 
typically  concentrate  on  matching  engine  representative  Reynolds 
number  (Re).  The  facility  for  this  work  is  a  transonic  facility  that 
matches  closely  with  engine  Reynolds  number,  as  well  as  Mach 
number  (Ma).  The  high  velocity,  slightly  compressible  flow  at  the 
inlet  of  the  cascade  poses  a  problem  in  generation  of  high  levels  of 
freestream  turbulence. 

This  paper  focuses  on  the  development  of  a  turbulence  grid 
designed  to  produce  high  intensity,  large  length-scale  freestream 
turbulence  in  high-speed  flow.  The  generator  was  designed  to  produce 
freestream  turbulence  with  a  turbulence  intensity  (Tu)  of 
approximately  15-20%  with  an  integral  length  scale  of  3-4  cm.  This 
desired  turbulence  level  (intensity  and  scale)  was  dctcmiincd  by 
matching  results  by  Van  Fossen  [1].  He  used  a  60®  section  of  a  GE  90 
combustor  in  a  wind  tunnel  to  model  the  exit  flow  of  a  combustor, 
Including  inlet  swirl  vanes,  film  cooling  holes  and  dilution  holes.  All 
three  of  these  combustor  components  contribute  strongly  to  the 
combustor  exit  turbulence  conditions.  Hot-wire  surveys  were 
performed  downstream  of  the  combustor  section  to  determine 
turbulence  intensity  and  length  scale.  Combustion  was  not  included, 
but  previous  research  by  Moss  and  Oldfield  [2]  showed  that  turbulence 
intensity  and  length  scale  were  not  significantly  affected  by  the 
presence  of  combustion.  Van  Fossen  measured  turbulence  intensities 
as  high  as  20%  with  length-scales  between  1  and  1.5  cm.  Using  these 
results,  the  desired  length  scale  was  based  on  matching  the  ratio  of 
length  scale  to  blade  chord  length  (Aj^/c)  of  the  transonic  turbine 
cascade  facility  to  the  blades  of  a  GE  90  engine. 

In  performing  this  research,  a  methodical  approach  was  taken. 
First,  grid  turbulence  was  generated  to  compare  with  previously 
reportwi  data  under  low  speed  conditions.  Second,  the  grid  was  tested 
in  the  transonic  wind  tunnel  to  quantify  the  turbulence  and  compare 
the  results  with  the  low  speed  testing.  Third,  based  on  inadequate 
turbulence  results  of  the  initial  design,  the  grid  was  modified  to  the 
current  design.  Results  of  measured  mean  and  turbulent  quantities  will 
be  presented. 

REVIEVV  OF  GENERATING  TECHNIQUES 
Mock  Combustor/Jets  In  Crossflow 

Ames  [3]  designed  a  mock  combustor  where  air  was  directed 
through  a  liner  into  the  mainstream  flow.  The  main  flow  is  directed 
along  the  walls  of  the  liner  and  two  rows  of  holes  were  used  to 
simulate  dilution  injection  flows.  Turbulence  levels  of  13%  were 
achieved  with  an  integral  length  scale  of  1.6  cm.  The  ratio  of  length 
scale  to  true  chord  length  was  AJe-O.  II.  A  similar  combustor 
simulator,  designed  by  Wang,  Goldstein  and  Olson  (4],  had  a  front 
panel  on  the  inlet  to  the  simulator  with  four  slots  and  guide  vanes  on 
the  inside  to  force  a  su'irling  flow  along  the  sidewalls  of  the  simulator. 
Rows  of  holes  were  placed  along  the  sidewalls  to  simulate  cooling 
flow  and  dispersion  jets.  The  highest  turbulence  parameters  achieved 


were  a  turbulence  intensity  of  18%  and  length  scale  of  8  cm 
(A/c=0.43). 

Schauer  and  Pestian  [5]  designed  a  turbulence  generator  where 
secondary  flow  was  injected  into  the  mainstream  from  two  opposing 
rows  of  holes  located  upstream  of  the  test  section.  The  jet  to 
freestream  velocity  ratio  was  14.  Turbulence  parameters  achieved 
were  1 7%  turbulence  intensity  and  length  scale  was  8  cm  [6]. 

Thole,  et  al.  [7]  used  small,  high  velocity  norma)  jets  injected  into 
cross  flow.  A  splitter  plate  was  placed  between  the  jets  to  prevent 
interaction  of  the  jets.  Further  improvements  using  a  downstream 
cylinder  were  made  by  Harrington  ct  al.  [8]  to  reduce  bulges  in  the 
mean  flow.  The  turbulence  parameters  achieved  were  up  to  20% 
intensity  and  2-3  cm  length  scale.  Later  work  by  Thole  led  to  the 
design  of  a  large-scale  combustor  simulator  in  a  low  speed,  close-loop 
wind  tunnel  facility.  The  combustor  simulator,  detailed  in  Barringer  ct 
al.  [9]  produced  turbulence  intensities  of  1 8%  and  length  scale  of  6.5 
cm  (A/c=0.1 1),  using  large  dilution  Jets. 

Grid  or  Bar  Generators 

Mchendalc  and  Han  [10]  used  an  active  grid  made  of  hollow 
brass  tubes  with  holes  drilled  at  the  grid  intersections.  Blowing  was 
done  uniformly  in  the  co-flow  (downstream)  direction  only,  with 
injection  ratio  (injection  mass  flow  ratc/frecslream  mass  flow  rate)  of 
2.5%.  Turbulence  intensity  of  12.9%  was  reported,  but  flow 
uniformity  was  reportedly  not  good. 

Sauer  and  Thole  [1 1]  and  Radomsky  and  Thole  [1 2]  developed  a 
turbulence  generator  with  simultaneous  upstream  and  dovvnstream 
blowing  from  vertical  bars.  The  generator  was  designed  to  produce 
turbulence  up  to  20%  intensity  and  a  range  of  length-scales  by 
adjusting  bar  spacing  and  blowing  rates. 

Polanka  and  Bogard  [13]  developed  a  turbulence  generator  which 
creates  turbulence  up  to  20%  in  a  passive  bar  arrangement.  The 
generator  provided  high  intensity,  large-scale  turbulence  with  good 
flow  uniformity  by  designing  a  bar  grid  with  large  flow  blockage. 

Boyle  et  al.  [14]  produced  turbulence  using  parallel  round  bars 
with  relatively  large  flow  blockage  and  the  option  of  active  blowing. 
Turbulence  of  greater  than  10%  was  produced  with  approximately  1.5 
cm  length  scale.  It  was  shown  that  active  blowing  significantly 
improved  flow  uniformity  but  reduced  overall  turbulence.  Giel  el  al. 
[15]  produced  turbulence  in  a  transonic  facility  using  an  active 
blowing  grid  of  square  bars.  The  grid  incorporated  blowing  in  the 
upstream  direction  only  at  a  rale  of  approximately  5%  of  the  cascade 
mass  flow  rate.  Blowing  air  supply  pressure  to  the  grid  was  125  psig 
(965  kPa).  Turbulence  intensity  and  length  scale  were  10%  and  2.2 
cm  (A/c=0.17),  respectively.  The  work  by  Giel  was  of  significant 
interest,  since  the  facility  is  a  transonic  tunnel  with  inlet  Mach  number 
similar  to  the  blowdown  facility  for  this  study.  This  was  the  only 
other  study  producing  high  intensity,  large-scale  turbulence  in  a 
transonic  facility. 

EXPERIMENTAL  APPARATUS  AND  INSTRUMENTATION 
Low  Speed  Facility 

The  setup  for  low  speed  testing  of  the  grid  consisted  of  a  5  HP 
blower  attached  to  several  transition  pieces,  followed  by  a  lest  section. 
The  blower  created  a  mainstream  flow  velocity  of  approximately  20 
m/s  at  the  entrance  to  the  test  section,  whereas  the  inlet  velocity  of  the 
transonic  wind  tunnel  is  approximately  120  m^s.  The  low  speed  test 
setup  was  designed  such  that  the  inlet  geometry  was  identical  to  the 
transonic  facility  from  the  grid  downstream  to  the  test  section.  The 
cross-sectional  area  at  the  grid  location  is  38.1  cm  (15  in.)  by  22.86  cm 
(9  in.)  and  is  followed  by  a  2-D  contraction  to  the  test  section  area.  A 
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Plexiglas  test  section  5 1  cm  (20  in.)  in  length,  15.2  cm  (6  in.)  wide  and 
30.5  cm  (12  in.)  high  was  fabricated  to  lake  data  downstream  of  the 
turbulence  grid.  The  cross  sectional  area  of  the  test  section  was  also 
designed  to  match  the  area  of  the  transonic  facility.  A  flow 
straightener  section  was  installed  in  the  setup  to  lower  background 
turbulence  levels  to  less  than  2%  upstream  of  the  turbulence  grid. 
Blowing  air  for  the  turbulence-generating  grid  was  supplied  through 
high-pressure  lines  to  a  plenum  tank.  The  plenum  supplied  air  to  each 
of  the  three  blowing  bars  of  the  grid. 

Transonic  Wind  Tunnel  Facility  and  Cascade  Test  Section 

The  experimental  facility  in  which  the  turbulence  grid  was 
ultimately  to  be  designed  for  use  was  a  transonic  blowdown  cascade 
wind  tunnel.  Previous  research  in  this  facility  has  been  published  by 
Holmberg  fl6],  Nix  ei  al.  [17,18),  Smith  et  al.  [19]  and  Popp  et  al. 
[20].  The  facility  allows  for  run  times  of  up  to  35  seconds  with  the 
inlet  total  pressure  maintained  by  a  feedback  control  scheme.  The 
facility  is  capable  of  heated  runs  by  way  of  a  healing  loop  built  into 
the  tunnel  upstream  of  the  cascade  lest  section.  A  schematic  of  the 
wind  tunnel  and  cascade  test  section,  including  the  turbulrmce  grid,  is 
shown  in  Figure  I. 

The  test  section  contains  a  cascade  of  four  full  and  two  half  first 
stage  rotor  blades  forming  five  passages.  The  turbine  blades  are  scaled 
up  three  times  to  allow  for  a  film  cooling  scheme  and  instrumentation 
on  both  surfaces  of  the  blades.  The  span  is  IS. 3  cm  (6  in.)  and  the 
aerodynamic  chord  is  13.7  cm  (5.4  in.).  The  Reynolds  number  based 
on  aerodynamic  chord  and  exit  conditions  is  about  610*.  The  inlet 
Mach  number  to  the  cascade  test  section  is  approximately  0.4,  with 
unheated  flow  mainstream  velocity  of  approximately  120  m/s.  For 
ease  of  installation  and  removal,  the  turbulence  grid  was  designed  to 
be  placed  directly  upstream  of  the  2-D  contraction  entering  the 
cascade  normal  to  the  inlet  flow  direction.  Again,  the  geometry  of  the 
cross-sectional  area  at  the  grid  location  downstream  to  the  test  section 
was  identical  to  the  low-speed  facility. 


Figure  1 .  Transonic  blowdown  facility  used  for  the  experiments 


Instrumentation 

Pressure,  temperature  and  velocity  were  measured  in  this  study. 
In  the  low  speed  facility,  blowing  bar  plenum  pressure  and 
temperature  were  measured.  The  bar  plenum  pressure  laps  were 
connected  to  0-100  psi  pressure  transducers.  Frccstream  velocity 
upstream  of  the  turbulence  grid  was  measured  using  a  pitot-static 
probe  connected  to  a  differential  pressure  transducer  with  a  range  of  2 
inches  of  water.  In  the  transonic  facility,  freestream  total  and  static 
pressures  were  measured  using  a  pitot-static  probe,  total  pressure 
probe  traverse  and  wall  static  pressure  taps.  AH  pressure  taps  were 
connected  to  0-15  psi  differential  pressure  transducers.  Temperature 
in  the  transonic  facility  was  measured  using  a  type  K  total  temperature 
probe  located  upstream  of  the  cascade  test  section.  Time-resolved 


velocity  was  measured  in  both  the  low  and  high-speed  facilities  using 
a  hot-wire  anemometer.  The  hot-wire  was  5  pm  in  diameter  and 
approximately  1.5  mm  in  length.  Measurements  were  performed  at  an 
overheat  ratio  of  1.5.  TTie  orientation  of  the  wire  was  such  that  the 
wire  was  parallel  to  the  span  of  the  turbine  blade.  Two  separate 
methods,  one  for  each  facility,  were  required  to  calibrate  the  hot-wire. 

The  hot-wire  anemometer  for  low  speed  testing  was  calibrated 
using  an  air  jet  calibrator.  The  hot-wire  probe  was  placed  in  a  low 
turbulence  air  jet  of  known  velocity  and  similar  operating  temperature 
as  the  low  speed  tunnel.  The  velocity  of  the  jet  was  varied  and  the 
output  voltage  of  the  wire  was  measured  and  recorded.  The  hot-wire 
voltage  vs.  velocity  data  was  then  fit  with  a  fourth  order  polynomial. 
Using  this  calibration  curve,  the  time-resolved  hot-wire  voltage 
measurements  performed  in  the  low  speed  tunnel  could  be  converted 
to  velocity  for  calculation  of  turbulent  quantities. 

Use  of  the  hol-wre  probe  to  measure  time  resolved  velocity  in 
the  transonic  facility  required  a  separate  calibration  of  the  probe.  Use 
of  the  probe  in  a  compressible  flow  field  does  not  allow  for  a  simple 
calibration  of  voltage  versus  velocity.  To  account  for  compressible 
flow,  the  wire  was  calibrated  in-situ  in  the  wind  tunnel  through  a 
blowdown  calibration.  The  wire  was  placed  in  the  wind  tunnel  and  the 
wire  output  voltage  was  recorded,  along  with  the  tunnel  total 
temperature  and  total  and  static  pressures  from  a  pitot-static  pressure 
probe  located  very  close  to  the  hot-wire  probe.  The  tunnel  was  run 
and  the  total  pressure  was  allowed  to  drop  to  very  low'  values, 
unchoking  the  cascade  and  providing  lower  Reynolds  numbers.  The 
data  from  these  measurements  was  then  used  to  plot  a  calibration 
curve  of  wire  Reynolds  number  (Re^)  and  Nusselt  number  (Nu’).  The 
data  was  then  curve  fit  with  a  power  law  to  give  an  equation  of  the 
type: 

iVw'^CRet  (i) 

where  the  wire  Reynolds  number  and  Nusselt  number  are  detenruned 
by: 

Re  =£a£A,  (2) 

‘  M 

Nu<-‘ - — -  0) 

All  wire  properties  were  constant  given  a  constant  temperature 
anemometer  that  maintains  wire  resistances  during  the  run. 

Timc-rcsolvcd  data  from  the  hot-wire  during  test  runs  could  then 
be  converted  from  hot-wire  voltage  to  Nusselt  number  and  the  wire 
Reynolds  number  could  be  determined  with  the  blowdown  calibration 
curve.  All  values  in  Equations  (2)  and  (3)  could  be  measured  during  a 
test  run,  but  the  product  (p^w )  could  not  be  decoupled  since  the  wire 
will  respond  to  fluctuations  in  velocity  as  well  as  density.  The  product 
of  density  and  velocity  {p„u  )  was  determined  for  every  hot-wire  data 
point  from  a  measured  hot-wired  voltage,  mean  static  pressure  and 
static  temperature.  Since  fluctuations  In  density  are  small  compared  to 
fluctuations  in  velocity,  the  instantaneous  velocity  could  be 
determined  by  dividing  the  Pn,M  product  by  the  average  flow  density 
determined  for  the  data-sampling  period  from  measurements  of  static 
pressure  and  temperature  to  yield  velocity: 

w=pmM/p  (4) 

Figure  2  shows  a  sample  blowdown  calibration  of  the  hot-wire.  Also 
included  in  the  plot  is  data  taken  from  a  small  blowdown  wind  tunnel 
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used  for  calibration  of  various  probes  for  use  in  compressible  flows. 
The  data  from  the  small  wind  tunnel  spans  a  larger  Reynolds  number 
range  than  that  achievable  during  blowdown  calibration.  The  data 
from  the  tw'o  different  tunnels  are  in  good  agreement. 


Figure  2.  Sample  hot-wire  blowdown  calibration  curve 


Data  AcQuIsItlon  and  Roductlon 

To  characterize  the  turbulence  from  the  turbulence  grid, 
calibration  of  the  hot-wire  was  necessary,  and  statistical  and  spectra! 
analyses  were  performed  on  the  time-resolved  velocity  data.  For  low 
speed  tests,  the  data  were  sampled  at  20  kHz,  filtered  at  10  kHz  and  a 
5V  DC  offset  was  applied  to  increase  dynamic  resolution.  All  data 
points  were  sampled  for  6  seconds  and  four  samples  were  taken  at 
each  point  for  averaging  purposes.  For  the  transonic  facility  testing, 
the  data  was  sampled  at  100  kHz  and  filtered  at  40  kHz,  providing 
shorter  sampling  periods  of  approximately  1 .3  seconds.  The  DC  offset 
was  adjust^  to  account  for  higher  flow  velocity  to  provide  better 
dynamic  resolution  of  the  signal. 

Statistical  and  spectral  analyses  included  calculation  of  the  mean 
velocity  (U),  fluctuating  component  of  velocity  (u),  root  mean  square 
(RMS)  of  the  fluctuating  component  of  velocity  («'),  autocorrelation 
(Rn),  probability  density  functions  (PDF),  and  power  spectra)  density 
(PSD). 

The  streamwisc  instantaneous  velocity  signal  («)  can  be 
decomposed  into  a  mean  and  fluctuating  component  as  follows: 

ti-U  +  u  (5) 

where  U  is  the  mean  velocity  and  u  is  the  fluctuating  component  of 
\^locity.  The  mean,  fluctuating  component  and  RMS  of  the 
discretized  velocity  data  were  calculated  according  to  the  equations  in 
the  nomenclature  section.  Using  these  values,  the  turbulence  intensity 
(Tu)  of  the  flow  at  a  given  point  was  calculated  by: 

Tu  =  ^  («) 

U 

The  autocorrelation  (Rn)  of  the  fluctuating  component  of 
velocity  (normalized  such  that  Rn  is  equal  to  !  at  zero  lag,  r=0),  was 
calculated  by: 

r  ^ 


V  J 


where  r-  j'At 

The  autocorrelation  was  used  to. determine  the  integral  time  scale  (T) 
by  integrating  under  the  curve  to  the  first  zero  crossing. 

r=]ff„(rVr  =  f;/f„,-Ar  (») 


where  No  is  the  point  of  the  first  zero  crossing. 

'fhe  streamwise  integral  length  scale,  which  is  representative  of 
the  largest  eddies  in  the  turbulent  fiowfield,  was  then  determined  by 
invoking  Taylor’s  hypothesis  of  frozen  turbulence: 

A,  =  u-r  (9) 


The  power  spectral  density  (PSD)  of  each  data  set  was 
determined  and  frequency  averaged  to  give  a  cleaner  curve  by 
averaging  PSD  data  over  a  frequency  range  of  Af=100  Hz.  The  power 
spectral  density  gives  an  indication  of  turbulent  energy  content  of  the 
flow  at  different  frequencies,  which  arc  inversely  proportional  to  the 
size  of  the  turbulent  eddies  in  the  flow  (the  wave  number,  tc,  is  equal 
to  27tfAJ).  The  PSD  was  also  used  to  check  for  vortex  shedding  from 
the  bars.  The  PSD  data  was  normalized  and  compared  to  the 


theoretical  Von  Karman  spectrum  equation  [21]: 


1  + 


(10) 


Figure  3  shows  an  example  of  measured  turbulence  data  compared  to 
the  Von  KarmAn  equation. 


Figure  3.  Low  speed  turbulence  PSD  compared  with  Von  KarmSn 
equation 

Using  these  statistical  and  spectra)  analysis  tools,  the  characteristics  of 
the  turbulence  could  be  determined. 

Uncertainty  Analysts 

Precision  uncertainty  of  measurements  of  velocity  was 
determined  using  the  computerized  uncertainty  analysis  method 
outlined  in  Moffat  [22].  Uncertainty  in  the  hot-wire  velocity 
measurements  in  the  low  speed  facility  was  estimated  to  be  ±1.5%. 
Uncertainties  in  turbulence  intensity  and  length  scale  in  the  low  speed 
facility  were  ±2%  and  ±6%,  respectively. 
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Uncertainty  in  velocity  measurements  with  the  hot-wire  in  the 
transonic  facility  was  estimated  to  be  ±2.3%.  Uncertainties  in 
turbulence  intensity  and  length  scale  in  the  transonic  facility  were 
±3.2%  and  ±6%,  respectively. 

LOW  SPEED  TESTING  OF  TURBULENCE  GRID 

As  mentioned  previously,  initial  testing  of  the  grid  was  performed 
in  a  low  speed  wind  tunnel.  The  initial  testing  served  two  purposes. 
First,  testing  was  done  to  benchmark  the  turbulence  measurement 
methodology.  Second,  results  were  used  to  obtain  baseline  numbers  to 
determine  what  levels  of  turbulence  to  expect  in  the  transonic  facility 
and  for  comparison  of  turbulence  intensity  and  length  scale  in  two 
different  flow  Helds. . 

The  initial  design  of  the  turbulence  grid  was  based  on  a 
combination  of  results  from  previous  grid  generated  turbulence  work. 
Simulation  of  turbulence  with  a  mock  combustor  or  jets  in  crossflow 
were  not  considered,  since  it  was  desired  to  employ  a  simplistic 
method.  The  initial  grid  design  for  this  work  was  an  actively  blown 
turbulence  generating  grid  incorporating  both  upstream  and 
downstream  blowing.  The  design  was  based  primarily  on  work  by 
Radomsky  and  Thole  [12]  and  by  Giel  [15].  Both  of  these  studies 
successfully  generated  large-scale,  high  intensity  turbulence  with 
actively  blown  turbulence  grids  in  low  and  high-speed  flow  facilities, 
respectively. 

The  design  that  was  studied  for  the  present  work  consisted  of 
three  2.54  cm  (1  in.)  hollow,  square  bars  spaced  7.62  cm  (3  in.)  apart 
between  bar  edges  with  no  vertical  bars.  Each  bar  has  nine  3.2mm 
(1/8  in.)  holes  spaced  2.54  cm  (1  in.)  between  centers  on  both  the 
upstream  and  downstream  sides  of  the  bars,  allowing  for  blowing  of 
high-pressure  air  in  both  directions.  The  upstream  blowing  was 
performed  to  increase  shear  stresses  upstream  of  the  grid  to  produce 
high  turbulence  levels.  Downstream  blowing  was  performed  to  fill  in 
the  wake  areas  behind  the  bars  to  create  better  flow  uniformity 
entering  the  cascade. 

Mean  and  root-mean-square  (RMS)  velocity  measurements  were 
taken  at  various  locations  throughout  the  test  section.  Flow  uniformity 
was  also  checked  (if  wakes  are  present  from  bars  on  grid)  in  the  test 
section.  This  traverse  was  performed  at  13  bar  widths  downstream  of 
the  grid.  Flow  uniformity  at  this  location  was  shown  to  be  ±2.5%. 

Four  locations  in  the  streamwise  direction  downstream  of  the  grid 
were  chosen  to  take  data  in  the  center  of  the  passage  (behind  the  center 
bar)  to  calculate  the  streamwise  variation  of  the  turbulence  intensity 
wid  length  scale.  These  locations  are  detailed  in  Table  1 . 

Table  1:  Low  Speed  Testing  Data  Acquisition  Locations 


1  Location  * 

Olstanc*  Downstream  of  Grid  1] 

tx/d)  1 

1  ... 

9.5  i 

1  2 

13.5  1 

I  3 

17.5  1 

LJ__ 

21.5  1 

The  mass  flow  ejected  from  (he  grid  holes  was  adjusted  in  an 
attempt  to  reach  the  turbulence  intensity  and  length  scale  desired. 
Momentum  flux  ratio  (I)  was  believed  to  be  the  blowing  parameter  to 
match  since  turbulence  is  generated  from  shear  stresses,  which  balance 
with  momentum.  Figure  4  -  Figure  6  plot  the  mean  velocity, 
turbulence  intensity  and  length  scale  versus  the  blow'ing  rate  (as 
momentum  flux  ratio).  The  maximum  value  of  blowing,  if  plotted  as 
mass  flux  ratio  (M)  would  be  approximately  a  mass  flux  ratio  of  14. 
At  a  momentum  flux  ratio  of  70,  the  mass  flux  ratio  is  approximately 


9.  The  mass  flux  ratio  is,  of  course,  not  linear  with  momentum  flux 
ratio.  Of  primary  importance  was  the  turbulence  intensity  and  length 
scale,  decay  of  turbulence,  flow  uniformity  and  the  effect  of  active 
blowing.  When  used  in  passive  mode  (i.c.  without  blowing  activated), 
the  grid  produced  turbulence  with  intensity  of  17%  at  a  distance  of  9.5 
bar  widths  downstream  of  the  grid  (Location  I),  which  decayed  to 
approximately  10%  at  a  streamwise  distance  of  21.5  bar  widths 
(Location  4),  with  length  scales  of  approximately  2  cm. 

The  mean  flow  in  the  test  section  without  the  grid  installed  was 
approximately  20  m/s  at  center  passage,  which  is  independent  of 
streamwise  measurement  location.  With  the  grid  installed.  Figure  4 
demonstrates  that  at  locations  near  the  grid  without  blowing,  there  is  a 
velocity  deficit  due  to  the  wake  of  the  grid  bars.  This  velocity  deficit 
decreases  in  the  streamwise  direction  and  away  from  the  influence  of 
the  bars.  With  blowing  activated,  the  mean  velocity  in  the  test  section 
will  increase  due  to  addition  of  mass  flow  from  the  blowing  bars. 
Based  on  holding  the  inlet  velocity  fixed,  a  conservation  of  mass 
analysis  predicts  the  expected  increase  to  be  from  20  m/s  with  no 
blowing  to  23  m/s  at  the  highest  momentum  flux  ratio.  Another  trend 
that  can  be  seen  from  the  data  is  the  effect  of  the  bar  wake.  Recall  that 
the  measurements  sho\vn  in  Figure  4  arc  taken  directly  behind  one  of 
the  bars  (cross-stream  surveys  were  not  taken).  At  low  momentum  flux 
ratios  there  is  clearly  a  reduced  velocity  at  the  first  position,  which  is 
actually  most  evident  at  a  momentum  flux  ratio  of  approximately  1=70. 
Even  though  there  is  blowing  from  the  bare  at  this  momentum  flux 
ratio,  it  is  not  enough  to  fill  in  the  bar  wakes  and  yet  there  is  enough 
upstream  blow'ing  to  appear  as  a  large  blockage  for  the  approaching 
flow.  As  a  result  of  this  flow  blockage,  the  velocity  deficit  of  the 
wake  appears  even  greater.  At  the  highest  momentum  flux  ratio  the 
data  show  there  is  a  fairly  constant  velocity  as  the  flow  progresses 
downstream,  indicating  the  mean  flowficld  is  likely  uniform.  The 
decrease  in  mean  velocity  at  the  highest  blowing  rate  is  due  to  the 
plenum  pressure  dropping,  most  likely  due  to  the  inlet  flow  to  the 
plenum  choking. 

The  turbulence  intensity  of  the  flow  without  the  grid  installed  was 
approximately  2%  upstream  of  the  grid.  Figure  5  demonstrates  that 
with  the  grid  in  place  without  blowing,  the  intensity  increases 
significantly  to  approximately  10-15%  depending  on  streamwise 
location.  Blowing  was  used  to  increase  turbulence  intensity  and  create 
larger  length  scales.  As  blowing  is  activated  and  increased,  there  is 
first  a  decrease  in  turbulence  intensity,  then  a  steady  increase  with 
blow  ing  rate.  The  initial  decrease  is  due  to  filling  in  the  wake  of  the 
bare,  which  essentially  reduces  the  effectiveness  (or  effective  bar 
width)  of  the  grid.  This  is  in  agreement  with  results  observed  by 
Boyle  [14]  where  active  blowing  at  low  blowing  ratios  is  seen  to 
slightly  decrease  turbulence  intensity. 

The  integral  length  scale  of  the  flow*  was  showm,  with  the 
exception  of  one  point,  to  steadily  increase  with  blowing  rate  as  can  be 
seen  in  Figure  6.  The  length  scale  at  higher  blowing  rates  is  seen  to 
settle  out  with  no  further  increase.  The  results  show  that  at 
approximately  1 8  bar  widths  downstream  of  the  grid,  the  turbulence 
intensity  is  as  high  as  17%  with  an  integral  length  scale  of  3.5  cm. 
This  corresponds  to  A^c=0.26  based  on  the  blade  chord  length  in  the 
transonic  facility. 
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Figure  4.  Plot  of  mean  velocity  at  4  streamwise  locations  at  different 
blowing  rates 


Momentum  Flux  Ratio  (I) 

Figure  5.  Plot  of  turbulence  intensity  at  4  streamwise  locations  at 
different  blowing  rates 


Figure  6.  Plot  of  integral  length  scale  at  4  streamwise  locations  at 
different  blowing  rates 

Although  the  grid  only  consisted  of  horizontal  bars,  which  are 
large  compared  to  the  size  of  the  test  section,  the  decay  of  turbulence 
intensity  without  blowing  was  shown  to  agree  well  with  published 
correlations  (Baines  and  Peterson  (23]).  The  decay  of  turbulence 
intensity  and  dispersion  of  length-scale  can  be  determined  by: 


{jJ 

where  c  and  a  arc  constants  given  as  0.80  and  0.20,  respectively,  x  is 
the  streamwise  distance,  d  is  the  bar  width,  and  n  and  m  are  exponents. 
The  value  of  n  is  typically  -5/7  and  m  is  between  0.5  and  0,56 
according  to  the  published  data. 

Figure  7  and  Figure  8  demonstrate  the  variation  of  turbulence 
intensity  and  length  scale  in  the  streamwise  direction,  respectively. 
Results  are  presented  for  both  the  low-speed  and  high-speed  testing  of 
the  grid,  although  the  high-speed  results  will  be  discussed  later  in  this 
work.  It  can  be  seen  in  Figure  7  that  the  low  speed  turbulence 
Intensity  data  without  blowing  agree  very  well  with  the  published 
correlations.  The  decay  data  with  blowing  is  seen  to  deviate  from  this 
correlation,  with  greater  deviation  with  higher  blowing  ratio. 
Dispersion  of  length  scale  does  not  agree  as  well  with  the  correlation. 


Figure  7.  Variation  of  turbulence  intensity  in  the  streamwise  direction 
with  and  without  blowing 


Figure  8.  Variation  length  scale  in  the  streamwise  direction  with  and 
without  blowing 
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TESTING  OF  TURBULENCE  GRID  IN  TRANSONIC 
FACILITY 

With  initial  testing  of  the  grid  in  the  low  speed  facility  providing 
satisfactory  results,  the  initial  grid  design  was  moved  to  the  high-speed 
facility  for  testing.  Hot-wire  data  were  taken  at  sireamwise  locations 
near  those  taken  in  the  low  speed  facility.  Data  acquisition  locations 
for  the  (lianstMik:  facility  are  iUusuated  in  Hguie  9.  The  location  are 
shown  as  T0-T4  to  avoid  confusion  with  the  locations  in  the  low-speed 
testing. 


2-D 

Contraction 


Rgurc  9.  Hot-wire  measurement  locations  in  transonic  facility 


The  turbulence  grid  was  initially  operated  in  passive  mode  in  the 
transonic  facility.  In  passive  mode,  the  turbulence  results  firom  the 
grid  did  not  match  the  levels  measured  in  the  low  speed  facility, 
generating  test  section  inlet  (location  TX))  turbulence  intensity  of  only 
5%  as  compared  to  approximately  \2%  in  the  low  speed  facility* 
which  would  decay  to  even  lower  levels  before  reaching  the  leading 
edge  of  the  turbine  blade  cascade.  Matching  active  blowing  to  the 
level  achieved  in  the  low  speed  facility  was  not  possible,  since 
matching  momentum  flux  ratio  of  the  blowing  would  require  bar 
plenum  blowing  pressures  exceeding  the  levels  achievable  within  the 
transonic  facility.  It  was  determined  that  blowing  pressure 
approaching  6(X)-700  PSIG  would  be  needed  to  match  bar  blowing 
momentum  flux  ratio  with  the  low  speed  tests  due  to  the  high  velocity 
and  higher  density  in  the  frecsircam.  Operating  in  passive  mode  (i.e. 
no  blowing),  the  grid  open  area  was  80%  and  the  bar  Reynolds  number 
was  much  higher  than  in  the  low  speed  test  rig  (approximately  8  limes 
higher),  thus  the  flow  experienced  smaller  cffcGiivc  blockage  and 
lowxr  turbulence  levels  when  compared  to  the  low  speed  tests. 

In  an  effort  to  achieve  more  desirable  turbulence  levels  in  the 
transonic  facility,  the  turbulence  generator  v«is  modified  by  increasing 
the  bar  width  to  5.08  cm  (2  in.)  in  order  to  decrease  the  open  area  of 
the  grid  to  60%.  This  modification  was  made  based  on  results 
achieved  by  Polanka  and  Bogard  (1 3)  using  large  bars.  This  modified 
grid,  operated  in  passive  mode,  generated  cascade  inlet  turbulence 
levels  approaching  15%.  with  an  integral  length-scale  of  2  cm. 
Blowing  was  performed  with  the  new  grid;  however,  the  turbulence 
intensity  decreased,  which  is  consistent  with  results  reported  by  Boyle 
[14].  Blowing  at  low  pressure  fills  in  the  wake  and  thus  reduces 
turbulence  intensity.  Also,  the  upstream  blowing  at  low  pressure  does 
not  generate  high  turbulence  intensity.  Again,  in  order  for  blowing  to 
have  a  positive  effect  on  turbulence  generation,  blowing  pressures  of 
higher  than  what  is  achievable  in  the  transonic  facility  arc  needed. 

Measurements  of  the  sireamwise  decay  of  turbulence  were 
performed  to  determine  the  variation  in  both  turbulence  intensity  and 


length  scale  as  the  turbulent  flow  progresses  from  the  inlet  of  the 
cascade  toward  the  entrance  to  the  cascade  blade  passages.  These  tests 
were  performed  with  a  coarse  spacing  (only  4  locations)  but 
demonstrated  that  the  turbulence  intensity  decays  to  approximately  10- 
12%  approaching  the  entrance  lo  the  cascade  passages  (location  T4). 
This  data  is  included  in  Figure  8  and  Figure  9.  TTie  turbulence  decay 
did  not  match  well  with  grid  generated  tuibulencc  decay  correlations. 
The  strong  decay  in  lurbwlencc  intensity  may  be  doc  to  nofi-tmiibim 
turbulence  in  the  cross-stream  diredion,  or  may  be  a  result  of  effects 
from  the  blade  (i.e.  local  velocity  is  higher  due  to  acccleratipn  of  the 
flow  on  the  suction  surface  of  the  blade).  Future  experiments  will  be 
performed  to  measure  aoss-siream  variation  in  tu^ulencc  intensity 
and  further  streamudse  decay  of  turbulence.  The  length  scale  of  the 
turbulence  docs  not  grow  to  any  extent  in  the  short  distance  measured, 
with  integral  length-scales  of  approximately  2  cm  (A^/c=0.15)  being 
measured  at  different  sireamwise  locations.  Turbulence  parameters  lor 
different  grid  designs  (vrilh  and  without  blowing)  and  measurement 
locations,  including  selected  data  from  the  low  spei^  testing,  arc 
shown  in  Tabic  2  below. 


Table  2:  Turbulence  Generation  Results 


Total  pressure  traverse  measurements  were  performed  behind  the 
turbulence  grid,  approximately  10.2  cm  (4  in.)  upstream  of  the  leading 
edge  of  the  center  cascade  passage  in  order  to  measure  flow  uniformity 
behind  the  modified  grid.  Each  traverse  covered  about  half  the 
passage  in  the  pilchwisc  direction.  The  modified  grid  has  an  open  area 
of  60%,  with  40%  flow  blockage  area,  and  thus  there  was  concern 
with  Inlet  flow  non-uniformity  due  to  large  bar  wakes.  Total  pressure 
surveys  demonstrate  that  although  there  is  a  significant  drop  in  total 
pressure  across  the  grid  (approximately  7%),  the  total  pressure  of  the 
flow  approaching  the  leading  edge  of  the  cascade  is  relatively  uniform, 
with  total  pressure  uniformity  of  ±1%.  The  main  tunnel  flow  pas^.s 
through  a  2: 1  area  contraction  downstream  of  the  grid  before  entering 
the  cascade  test  section,  which  Is  believed  to  help  in  merging  and 
mixing  of  the  relatively  large  bar  wakes,  thus  increasing  flow 
uniformity.  Results  of  the  total  pressure  traverse  for  both  the  2.54  cm 
(1  in.)  and  5.08  cm  (2  in.)  bar  grid  designs  are  shown  in  Figure  10.  It 
can  be  seen  that  the  total  pressure  with  the  2.!>4  cm  bars  is  non- 
uniform,  with  wakes  present  without  blow'ing.  With  blowing,  the  total 
pressuie  traverse  data  shows  that  there  is  a  jet  from  the  bars.  The  data 
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for  the  5.08  cm  bars  with  blowing  shows  that  there  is  a  smaller  total 
pressure  drop  across  the  grid,  but  the  flow  uniformity  is  not  greatly 
improved.  Since  active  blowing  actually  lowered  the  turbulence 
intensity  slightly,  blowing  was  not  used  once  the  grid  was 
implemented  for  wind  tunnel  heat  transfer  tests. 


CONCLUSIONS 

The  development  of  a  turbulence  generator  for  use  in  a  high* 
speed  flow  facility  has  been  documented.  The  generator  was  tested  in 
a  low'Spced  environment  with  good  results  in  turbulent  flowficid  and 
flow  uniformity.  Results  from  testing  in  a  high-speed  flowfleld  were 
in  poor  agreement  with  low-speed  results  due  to  higher  Reynolds 
number  (from  higher  velocity  and  density)  providing  inadequate  flow 
blockage.  The  momentum  flux  ratio  from  low  speed  testing  of  the  grid 
could  not  be  matched  in  transonic  facility.  The  grid  was  modified  to 
create  higher  blockage  and  used  as  a  passive  grid.  Results  from  the 
modified  grid  showed  lower  turbulence  than  originally  desired,  but 
with  reasonably  high  intensity  (10-12%  near  the  entrance  to  the 
cascade  passages),  large  length  scale  (2  cm)  compared  to  blade  chord 
and  good  total  pressure  uniformity  downstream  of  the  grid  near  the 
entrance  to  the  cascade  passages.* 
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Motivation 


Goidstein,etal.  (1983) 


Early  studies  identified  high  levels  of  freestream 
turbulence  exiting  the  combustor  ->  Effects  on  heat 
transfer? 


Kuotmos  and  McGuirk  (1989) 


Objectives 

;  The  primary  objectives  of  the  current 
research  effort  are: 

-  To  experimentally  study  the  effects  of 
combustor  exit  freestream  turbulence  on 
turbine  blade  heat  transfer 

-  Develop  and  test  a  turbulence  grid  to  produce 
high-intensity,  large-scale  freestream 
turbulence  in  a  transonic  turbine  cascade 


irelniaTech 


Turbulence  generating  techniques 


Investigator 


Turbulence  Length  Scale 


Level 


{AJO 


Method 


Flow 

Condition 


Mehendale  and  Han 
(1992) 

13% 

3 

" 

Blown  Grid 

Low  speed 

Thole,  et  al.  (1994) 

20% 

2-3 

• 

Jets  in  crossflow 

Low  speed 

Holmberg  (1996) 

5% 

0.5 

^0.1 

Grid 

Transonic 

Schauer  and 

Pestian  (1996) 

17% 

8 

• 

Jets  in  crossflow 

Low  speed 

Ames  (1997) 

13% 

1.6 

0.11 

Mock  combustor 

Low  speed 

Boyle,  etal.  (1998) 

10% 

1.5 

- 

Blown  and  unblown 
grid 

Low  speed 

Wang.etal.  (1998) 

18% 

8 

0.43 

Mock  combustor 

Low  speed 

Radomsky  and 

Thole  (1999) 

20% 

5-6 

" 

Blown  bar  array 

Low  speed 

Polanka  and 

Bogard  (1999) 

20% 

5 

- 

Passive  bar  array 
Oarge  blockage) 

Low  speed 

Giel.  etal.  (2000) 

10% 

2.2 

0.17 

Active  grid 

Transonic 

Baninger.  et  al. 
(2001) 

18% 

6-7 

0.11 

Combustor 
simulator  with  large 
jets 

Low  speed 

Current  turbufence  generator 
development 

-  Initial  design  was  an  actively  blown  bar  arrangement  based 
on  low-speed  work  by  Radomsky  and  Thole  (1998)  and  work 
in  a  transonic  facility  by  Giel,  et  al.  (2000) 


-  Modification  of  the  initial  design  was  performed  based  on 
work  of  Polanka  and  Bogard  (1999) 


Experimental  Setup 


Two  facilities  were  used  to  perform  testing  of 
the  turbulence  grid: 

-Low-speed  blower  wind 
tunnel  (Inlet  velocity  20  m/s) 


-Transonic  cascade  facility 
(inlet  Mach  number  of  M=0.4, 
approximately  120  m/s)  j 


Tuibuicncc 

-  ^  —  _  Grid  TwtSedioii 

HeetDtiJttngw  FhpwrVikfe  .  . 

.  /-  I.  l-J 


sHaab. 
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Results  of  low-speed  testing 


Streamwise  variation  in  turbulence  intensity 


mvireiniaTech 
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Testing  in  transonic  facility 

|urbulenwgrid  ^  Pitchwise  survey  streamwise  survey 

\J - /  / 


Contraction 


7^^ 


Cascade  Test  Section 


rginiaTech 


vMiMVA  ^rnicwociinrmvii  viovn^iv 


Cascade  Design: 

Blade  span  (s)  •>  16.3  cm 
Aero.  Chord  (c)  •>  13.7  cm 
Pitch  (p)-12.2  cm 


Blade  1-xA>»  7.26 
Blade  2 -xn>>  9.26 


Results  of  Transonic  Testing 


Comparison  to  Low  Speed  Testing 


-Passve  Gnd(l=0) 

1  -in  Passive  Transon  b 
1 -in  Transonic  (N10) 
-Baines  &  Peterson 


Distance  from  Grid  -  x/D 


IniaTeeJ 
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Grid  Modification 


•  Based  on  turbulence  measurements  with  the  initial 
,5  grid  design  in  the  transonic  facility,  the  grid  was 
modified 


•  Modification  of  the  initial  design  was  performed 
based  on  work  of  Polanka  and  Bogard  (1999)  by 
increasing  bar  sizeio  5.08  cm  (2  in) 


♦  T urbulence  measurements  in  the  transonic  facility 
were  performed  with  the  modified  grid 


Streamwise  turbulence  uniformity 


Streamwise  Distance  from  Grid  (x/D) 

4  5  6  7  8  9  10  11  12 


Appendix  A 


137 


Appendix  A 


138 


Conclusions 

-  The  development  of  a  turbulence  generator  for  use  in  a 
transonic  flow  facility  has  been  documented 

-  The  generator  was  tested  in  a  low-speed  environment 
with  good  results  in  turbulent  flowfield  and  flow  uniformity 

-  Results  from  testing  in  a  high-speed  flowfield  were  in 
poor  agreement  with  low-speed  results 

-  The  grid  was  modified  to  create  higher  blockage  and 
used  as  a  passive  grid 

-  Modified  grid  generates  10-12%  turbulence  near  the 
entrance  to  the  cascade  passages  with  2  cm  length  scale 

-Spanwise  and  pitchwise  uniformity  of  turbulence  was 
shown  to  be  good  with  no  vortex  shedding  measured 


Testing  in  transonic  facility 
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Measurement  Uncertainty 

The  following  section  documents  measurement  uncertainties  of  many  of  the 
reported  data  in  this  dissertation.  Results  presented  in  this  work  are  dependent  on 
experimental  measurements,  making  determination  of  the  uncertainty  in  these 
measurements  important.  The  primary  quantities  of  concern  are  the  steady  heat  transfer 
coefficients,  hot-wire  velocity  measurements  (primarily  u'nns),  the  integral  length  scale, 
and  the  model  predicted  values  of  increase  in  heat  transfer.  The  uncertainty  associated 
with  each  one  of  these  measured  and  calculated  values  will  be  presented.  The  uncertainty 
of  experimental  measurements  and  values  calculated  from  these  measurements  were 
determined  using  the  computerized  uncertainty  analysis  method  (perturbation  method) 
outlined  in  Moffat  (1988). 

B.1  Defining  Uncertainty 

Total  measurement  uncertainty  consists  of  two  components,  bias  uncertainty  (B) 
and  precision  uncertainty  (P).  Bias  uncertainty  is  error  that  is  consistent  from 
measurement  to  measurement  and  is  typically  a  function  of  errors  in  calibration  of 
transducers  and  other  measurement  devices.  Precision  or  random  uncertainty  is 
associated  with  run-to-run  variation  in  measurements.  The  total  uncertainty  (A)  is 
determined  by: 

A  =  Vb^  +  P^  (B.l) 

Many  of  the  bias  uncertainties  of  measurements  were  estimated  based  on  typical 
values  for  transducers,  thermocouples,  etc,  or  based  on  results  reported  in  previous 
research  using  similar  transducers  and  measurement  techniques. 
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B.2  Determining  Uncertainty  of  Measurements  and  Calculated 
Values 

Uncertainty  in  Heat  Transfer  Coefficients  (h) 

The  bias  uncertainty  for  the  measurements  of  heat  transfer  coefficients  is  a  result 
of  the  calibrated  sensitivity  of  each  sensor.  The  bias  uncertainty  for  the  heat  flux  sensors 
was  detailed  in  Smith,  et  al.  (1999)  and' in  Popp  (2000).  The  average  uncertainty  in  the 
heat  flux  sensors,  based  on  convection  and  radiation  heat  flux  tests  is  estimated  to  be 
approximately  ±6%.  Since  the  calculation  of  heat  transfer  coefficient  is  determined  from 
the  slope  of  a  plot  of  heat  flux  (q)  versus  Tt-Tw,  (see  section  3.2.2)  any  uncertainty  in 
thermocouple  measurements  can  only  shift  the  curve,  not  alter  the  slope  (refer  to  Figure 
3.5b).  The  bias  uncertainty  of  the  heat  flux  sensors  then  directly  translates  into  the  error 
in  heat  transfer  coefficient. 

The  precision  uncertainty  in  h  was  determined  through  repeated  measurements 
using  a  90%  confidence  interval  assuming  a  Student-t  distribution.  The  precision 
uncertainty  for  each  location  was  determined  by: 

p  =  =  (B.2) 

h  h'VN 

where  Sh  is  the  variance  of  the  heat  transfer  coefficient,  h  is  the  mean  of  the  heat  transfer 
coefficient  measurements,  N  is  the  number  of  samples  and  t  is  the  Student-t  variable. 

The  low  turbulence  precision  uncertainty  was  determined  fi’om  N=6  separate  tests 
(t=2.0I5)  and  the  high  turbulence  heat  transfer  coefficient  precision  uncertainty  was 
determined  fi'om  N=4  separate  tests  (t=2.353)  for  each  gage  location.  Table  B.l  below 
details  the  precision  uncertainty  determined  for  each  gage  location. 


Table  B.l  Precision  uncertainty  (P)  of  heat  transfer  coefficient  tests 


Gage  Location 

Low  Turbulence 

High  Turbulence 

SSI 

1.8% 

2.3% 

PSl 

2.8% 

1.7% 

PS2 

1.9% 

1.7% 
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The  uncertainty  in  measured  increase  in  heat  transfer  coefficient  due  to  turbulence  is  only 
a  function  of  the  precision  uncertainties  in  Table  B.l  since  both  low  and  high  turbulence 
measurements  are  affected  by  the  same  bias  uncertainty. 

From  equation  B.l  and  the  precision  uncertainties  in  Table  B.l  and  using  a  bias 
uncertainty  of  ±6%  for  each  gage  location,  the  total  uncertainty  for  each  measurement 
location  for  both  low  and  high  turbulence  heat  transfer  coefficients  was  calculated  and  is 
presented  in  Table  B.2. 


Table  B.2  Total  uncertainty  (A)  of  heat  transfer  coefficient  tests 


Gage  Location 

Low  Turbulence 

High  Turbulence 

SSI 

6.3% 

6.4% 

PSl 

6.6% 

6.2% 

PS2 

6.3% 

6.2% 

Uncertainty  in  hot-wire  velocity  and  turbulence  measurements 

Uncertainty  in  measurements  of  unsteady  velocity  is  of  primary  concern  in  the 
form  of  the  RMS  of  fluctuating  velocity  (u'rms).  The  RMS  of  velocity  is  used  both  in 
determination  of  the  turbulence  level  as  well  as  in  the  analytical  model  calculations.  Bias 
uncertainty  for  the  hot-wire  measurements  is  a  function  of  the  hot-wire  calibration  and 
the  variables  on  which  the  calibration  is  dependent.  The  bias  uncertainty  was  determined 
by  the  perturbation  method  outlined  in  Moffat  (1988).  Briefly,  the  perturbation  method 
combines  the  effects  on  the  dependent  variable  (quantity  of  interest)  of  perturbing  the 
values  of  the  independent  variables  by  their  respective  uncertainties.  For  the 
measurement  of  hot-wire  velocity,  this  can  be  presented  by: 


Au  = 


(B.3) 


where  for  the  hot-wire  velocity  calibration,  Xj  e{T,,T^„,P,,P5,V},  where  the  variables 
are  the  total  temperature  of  the  flow,  the  temperature  of  the  hot-wire,  the  freestream  total 
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and  static  pressures  and  the  hot-wire  voltage.  The  hot-wire  calibration  technic|ue  is 
detailed  in  Nix,  et  al.  (2002)  in  Appendix  A.  These  variables  were  perturbed  by  values 
consistent  with  uncertainty  for  each  of  the  transducers  used  to  measure  the  values  (±1°C 
for  thermocouple  measurements  and  ±0.15  psi  for  pressure  transducers)  and  an  estimate 
of  the  hot-wire  voltage  uncertainty  of  ±2  mV.  The  perturbation  analysis  was  performed 
by  the  computerized  method  presented -in  Moffat  (1988)  using  a  spreadsheet.  The  bias 
uncertainty  in  hot-wire  velocity  was  determined  to  be  ±2.3%. 

Precision  uncertainty  of  u’ms  was  determined  in  a  similar  fashion  as  for  the  steady 
heat  transfer  coefficients,  using  a  90%  confidence  level  assuming  a  Student-t  distribution: 


P 


AUnn,  _  t-S„> 

U'nns  u'^,-VN 


(BA) 


The  precision  uncertainties  were  determined  from  several  repeated  measurements 
at  each  measurement  location.  The  results  of  equation  B.4  were  combined  for  all 
measurements  and  the  precision  uncertainty  was  determined  to  be  ±4.4%.  From  equation 
B.l,  the  total  uncertainty  in  the  measurement  of  u'rms  was  determined  to  be  ±5%. 
Through  perturbation  analysis,  the  total  uncertainty  in  the  measurement  of  turbulence 
intensity  was  determined  to  be  ±5.5%. 

Uncertainty  in  integral  length  scale  (Ax) 

The  integral  length  scale  was  determined  by  integrating  the  autocorrelation  of 
streamwise  velocity  to  determine  the  integral  time-scale  and  then  invoking  Taylor’s 
hypothesis  and  multiplying  the  integral  time  scale  by  the  mean  velocity.  Error  in  the 
integral  time  scale  is  a  function  of  determining  the  zero  crossing  of  the  autocorrelation 
function.  Perturbation  of  the  integral  time  scale  was  performed  based  on  a  typical 
autocorrelation  curve  and  the  uncertainty  in  integrating  the  autocorrelation  was  estimated 
to  be  about  ±2%.  The  uncertainty  associated  with  invoking  Taylor’s  hypothesis  and 
multiplying  by  the  mean  velocity  was  estimated  by  combining  the  uncertainty  in  the 
integral  time  scale  (±2%)  and  the  mean  velocity  (±2.3%)  through  perturbation  analysis. 
The  bias  uncertainty  was  then  estimated  to  be  ±3% 
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Precision  uncertainty  in  calculating  the  integral  length  scale  was  determined  from 
repeated  measurements  of  the  integral  length  scale  at  each  gage  location.  The  precision 
uncertainty  was  determined  to  be  ±5.3%.  From  equation  B.l,  the  total  uncertainty  in  the 
measurement  of  Ax  was  determined  to  be  ±6.1%. 

Uncertainty  anal3^ical  model  predicted  heat  transfer  increase  (Aht) 

The  final  uncertainty  to  be  quantified  is  the  increase  in  heat  transfer  coefficient 
predicted  by  the  analytical  model  (Aht),  which  is  a  function  of  both  u'msand  Ax 
(neglecting  small  uncertainty  in  determining  the  thermal  conductivity  (k)  and  thermal 
diffusivity  (a)).  The  bias  uncertainty  in  the  model  predictions  was  determined  through 
perturbation  analysis  using  the  total  uncertainty  calculated  for  the  RMS  of  velocity 
fluctuations  (u'rms)  of  ±5%  and  the  integral  length  scale  (Ax)  of  ±6.1%.  The  bias 
uncertainty  on  Aht  was  determined  to  be  ±4%. 

Precision  uncertainty  in  quantifying  the  value  of  Aht  was  determined  from 
repeated  measurements  at  each  gage  location.  The  precision  uncertainty  was  determined 
to  be  ±2%.  From  equation  B.l,  the  total  uncertainty  in  calculating  Aht  was  determined  to 
be  ±4.5%. 

B.3  Uncertainty  Summary 

Table  B.3  below  summarizes  the  total  uncertainty  determined  for  all  relevant 
parameters  in  this  work. 


Table  B.3  Summary  of  total  uncertainty  (A) 


Parameter 

Total  Uncertainty 

b 

(average  of  low  and  high  turbulence) 

6.3% 

u'rms 

5.0% 

Tu 

5.5% 

Ax 

6.1% 

Ah, 

4.5% 
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Appendix  C 


Computer  Code 


The  following  code  is  a  Matlab  M-file  which  was  used  to  reduce  hotwire  voltage 
data  based  on  recent  blowdown  calibration  and  perform  analyses  in  the  time  domain  and 
frequency  domain,  as  well  as  performing  analytical  model  analyses  in  both  the  steady  and 
unsteady  form: 


HWDATA062401.M: 

%  Matlab  M-file  to  reduce  Hotwire  Data  Using  06/22  Auspex  Wire  Blowdown 

Calibration 

clear  all;  close  all; 

%Sampling  Parameters: 

fsample=l 00000;  %  Sampling  frequency  of  100  kHz 
dt=l/fsample;  %  Sampling  period  (sec) 

N=1 3 1 069;  %Number  of  Samples 
T=N*dt;  %Total  Sampling  Period  (seconds) 
df^l/T;  %Frequency  Resolution 

t=[0:dt:T];  %Define  time  T=(Number  Samples)*(dt),  dt=l /(Sampling  Frequency) 
f=[0:df:fsample];  %Defme  frequency  range  from  0  to  fsample  in  increments  of  df 

Patm=94760/6894.8; 

Twire=464; 

Dwire=5E-6; 

%  DATA  REDUCTION  USING  BLOWDOWN  CALIBRATION 

%Read  in  hotwire  data  from  LeCroy: 

load  hw062401.asc;  %Read  hotwire  voltage 

V2=hw062401; 

V2=V2+6;  %Add  on  6  V  offset 
%Mach2=0.5;  %From  Mach  number  distribution  for  SS2 
Mach2=0.25;  %From  Mach  number  distribution  for  PSl 
%Mach2=0.3;  %From  Mach  number  distribution  for  PS3 
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Pt2=9.527068893;  %From  low  speed  data 
Ps2=(Pt2+Patm)/((l+0.2*Mach2^2)^3.5)-Patm; 

Tt2=299.4831878;  %Froin  low  speed  data 
Ts2=Tt2/(l+0.2*Mach2^2); 

Tm2=(Ts2+Twire)/2; 
rho2=(Ps2+Patm)*6894.8/(287*Tm2); 
k2=0.02414*(Tni2/273.16)^1.5*(473.16/(Tm2+200)); 
Nii2=V2.^2./((Tm2-Ts2)*k2).*(Tm2/Ts2)'^(-0. 1 7); 

Re2=(Nu2./4.6267).^(l/0.2755); 

mu2=0.00001 7 1 6*(Tm2/273. 1 6^1 .5 *(383 .7 1 6/(Tm2+l  1 0.556)); 
rU2=Re2.  *mu2  ./Dwire; 

U2=rU2./rho2; 

Vniean(2)=mean(U2); 
v2=U2-Vmean(2); 
fluct2=(mean(v2.'^2)).^.5 ; 

Tu2=fluct2A^mean(2); 

%Using(rhoU)VrhoU: 

rUp2=rU2-mean(rU2); 

Um2=Mach2*(l  .4*287*Ts2)^0.5; 

Tup2=((mean(rUp2.^2)).^0.5)./(rho2*Um2); 

%Spectral  Data 

%Calculate  Statistics  for  each  data  file: 

V22=(v2).^2;  %a^2 
V23=(v2).^3;  %a^3 
V24=(v2).M;  %aM  ■ 

%Calculate  RMS  for  each  file  (RMS=  (mean(a^2))^l/2): 

RMS2=(mean(V22)).^5; 

%Calculate  Kurtosis  for  each  file  (Kurtosis  =  (mean(aM))/((mean(a^2))''2)): 
Kurt2=(mean(V24))/(tnean(V22).^2); 

%Calculate  Skewness  for  each  file  (Skewness  =  (mean(a^3))/(mean(a^2))^3/2): 
Skew2=(mean(V23))/(niean(V22).^1.5); 

%PART  b):  PLOT  AUTOCORRELATION  CURVE  AND  DETERMINE  INTEGRAL 
LENGTH  SCALE: 

% _ 


%Calculate  autocorrelation  on  each  data  file: 

Rv2=xcorr(v2,'coefF); 

Rv2=Rv2(131070:262139); 
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%PIot  autocorrelations: 
figure(l); 

plot(t,Rv2,V); 

title('Figure  1 :  Autocorrelation  of  Velocity  Signal'); 
xlabelOtau  (s)');  ylabelCRvOtau)'); 
axis([0  0.002 -0.2  1.2]); 
grid  on; 

%Perform  Numerical  Integration  to  Calculate  Integral  Time  Scale: 

%First  zero  crossing  calculated  with  find  command: 
z2=fmd(abs(Ry2)<0.02&abs(Rv2)>=0); 

T2=sum(Rv2(i:z2(l)))*dt;  %T  is  integral  time  scale 
L(2)=T2*Vmean(2); 

%PART  c):  DETERMINE  AND  PLOT  PROBABILITY  DENSITY  FUNCTION  (PDF) 
OF  DATA  SETS: 

% 


%Create  PDFs  for  each  Data  file: 

%Create  bin  file  (50  bins) 

%Hotwire  Signal: 

dV2=((Vmean(2)+3*RMS2)-(Vmean(2)-3*RMS2))/50; 

m2=[(Vmean(2)-3*RMS2):dV2:(Vmean(2)+3*RMS2)]; 


N2-hist(U2,m2);  %Create  histogram  with  50  bins  whose  centers  are  defined  in  the  vector 
m 

pdf2=N2y(sum(N2)*dV2);  %Normalize  histogram  to  create  PDF  (divide  by  Ntotal  and 


%PlotPDFs: 

figure(2); 

bar(m2,pdf2,'r'); 

xlabel('Hotwire  Voltage  (V)');  ylabel('PDF  (NJ/(N_t_o_t_aJ*\DeltaV)'); 
title('Figure  2:  Probability  Density  Function  for  (PDF)  Velocity  Signal'); 
grid  on; 

%PART  d):  PLOT  NON_DIMENSIONAL  POWER  SPECTRAL  DENSITY  OF  DATA 
SETS: 

% 


%Perform  Discrete  Fourier  Transform  of  each  data  set: 
%Perform  FFT  and  scaling  of  each  data  set: 
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%Complete  Data  set: 

Vf2=ffl(v2)./(N+1);  %Scale  1/N  FFT  and  correct  for  double  sided  FFT  (multiply  by  2) 
Vf2m=2*Vf2.*conj(Vf2)./df; 

%  Do  frequency  AVGing: 
fory=l:(N+l)/100 
fa(y)=mean(f(y*  1 00-99:y*  1 00)); 

Vfma2(y)=mean(V  f2m(y*  1 00-99:y  *  1 00)); 
end 

%Plot  Power  Spectra: 

figure(3); 

loglog(f,Vf2m); 

xlabel('Frequency  (Hz)*); 

ylabel([Trequency  Averaged  PSD']); 

title('Figure  3:  PSD  of  Velocity  Signal*); 

grid  on; 

%PIot  Frequency  Averaged  Power  Spectra: 

figure(4); 

loglog(fa,Vfrna2); 

xlabel(Trequency  (Hz)'); 

ylabel(['Frequency  Averaged  PSD']); 

title('Figure  4:  Frequency  Averaged  PSD  of  Velocity  Signal'); 

grid  on; 

fa=transpose(fa); 

Vfrna2=transpose(Vfrna2); 

%  Normalize  with  Von  Karman: 

Vfmn=(Vfma2.*Vmean(2))/(RMS2.^2.*L(2));  %Non-dimensional  PSD 
E=4./(l+(8*pi.*fa.*L(2)./(3*Vmean(2))).^2).^(5/6); 

K=L(2).*faA^mean(2); 

%Plot  Power  Spectra  along  with  Von  Karman  equation: 
figure(5); 

loglog(K,Vfmn.K,E); 

xlabel('\Lambda_xf/U'); 

ylabel(['PSD  (E_v(\kappa)*U)/(u^2\Lambda_x)']); 
title('Figure  5:  Normalized  Power  Spectral  Density*); 
legendCPower  Spectral  Density','VonKarman  Spectra'); 
grid  on; 

%Model  Analysis 
%Load  heat  flux  data: 
load  HF062401.asc 

T  =  277.4791737;  %From  low  speed  data 
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%q=((HF062401*1000000)/100)/(105);  %For  SS2  location 

q=((Iff062401*1000000)/100)/(0.236426*T+249.4697)*  10000;  %For  PSl  location 
W/cm^2 

%q=((HF062401  ♦  1 000000)/!  00)/(0. 1 48635*7+236.5432);  %For  PS3  location 

meanq=niean(q); 

qp=(q-meanq); 

%Steady  model: 
urms  =  RMS2; 

Lx  =  L(2); 

Te  =  Ts2;  %Static  temp  of  flow 
ke  =  0.02414*(Te/273.16)^1.5*(473.16/(Te+200)); 
rho  =  (Ps2+Patin)*6894.8/(287*Ts2); 
ae  =  ke/(rho*  1004.5); 

qptsteady  =  abs(ke.*(Te  -  T)./(pi.*ae.*Lx./urms).'^0.5); 
qptsteadymax  =  abs(ke.*(Te  -  T)./(pi.*ae.*Lx./max(v2)).^0.5); 
dhtsteady  =  qptsteady/(Ts2-T); 

%Unsteady  Model 
%Heat  Flux: 

QlNfft(qp)./(N+l);  %Scale  1/N  FFT  and  correct  for  double  sided  FFT  (multiply  by  2) 
Qfm=2*Qf.*conj(Qf)./df;  %PSD  of  q  fluctuations 

%  Determine  magnitude  of  both  velocity  and  heat  flux  FFTs: 

Ufmag=(imag(Vf2).^2+real(Vf2).^2).^0.5;  %Magnitude  of  u  fluctuations  in  frequency 
domain 

Qfmag=(imag(Q0.^2+real(Qf).'^2).'^0.5;  %Magnitude  of  q  fluctuations  in  frequency 
domain 

%Apply  model  to  Ufmag: 

Lxu  =  (3.*Vmean(2))./(8.*pi.*f); 

Lxu  =  transpose(Lxu); 

qptu  =  abs(ke.*(Te  -  T)y(pi.*ae.*(Lxu./Ufmag)).^0.5); 

%Apply  Coherence  function  to  data 

%coherence  file  determined  with  HFMPHW2  code  fi’om  Holmberg  (1996): 

load  next062401.coh; 

fcoh  =  next062401(l,:); 

fcoh  =  transpose(fcoh(l  :((N+l)/2))); 

coh  =  next062401(2,:); 

coh  =  transpose(coh(l  :((N+l)/2))); 

qptcoh  =  qptu(l  :((N+l)/2)).*coh; 

%Plot  measured  magnitude  of  heat  flux  fluctuations  vs.  model  with  coherence: 
figure(6) 

loglog(f,Qfmag,'b',fcoh, qptcoh, V) 
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xlabel('Frequency  (kHz)'); 

ylabel(['Heat  Flux  Fluctuation  Magnitude,  Predicted  Heat  Flux  (W/cm''^2)']); 
title('Figure  6:  Model  Heat  Flux  Comparison  with  Magnitude  of  Measured  Heat  Flux'); 
axis([10e-l  40000  lOe-14  1]); 

legendCMeasured  Heat  Flux  Fluctuation  Magnitude', 'Unsteady  Model  Results'); 

%Plot  measured  heat  flux  PSD  vs.  model  with  coherence: 
figure(7) 

loglog(f,Qfm,'b',fcoh,qptcoh,'r') 
xlabel('Frequency  (kHz)'); 

ylabel(['Heat  Flux  Spectra,  Predicted  Heat  Flux  (W/cm'''2)']); 

title('Figure  7:  Model  Heat  Flux  Comparison  with  PSD  of  Measured  Heat  Flux'); 

axis([10e-l  40000  lOe-14  1]); 

legend('Measured  Heat  Flux  PSD','Unsteady  Model  Results'); 

%Methods  of  combining  unsteady  model  heat  flux  (section  3.4 .4.3): 

%1.  Plot  qtmeas  and  calculate  qtrms  and  dqt  total  from  measured  data  (no  modeling): 

qrms  =  sqrt(mean(qp.'^2));  %RMS  of  measured  q  fluctuations  calculated  in  time  domain 
qrmsf  =  sqrt(sum(Qfm(l:60000))*df);  %RMS  of  measured  q  fluctuations  calculated  in 
frequency  domain 

dqtotall  =  sqrt(2.*sum(Qfmag(l:6555).'"‘2)); 

%  Plot  measured  values,  both  PSD  and  Magnitude 
figure(8) 

loglog(f,Qfm,'b',f,Qfmag,'r’) 

xlabel('Frequency  (kHz)'); 

ylabel(['PSD  and  Magnitude  of  HF  Fluctuations']); 

title('Figure  D1 :  PSD  and  Magnitude  of  Fluctuations  of  Measured  Heat  Flux'); 
axis([10e-l  40000  lOe-14  1]); 

legendCMeasured  Heat  Flux  PSD','Magnitude  of  Heat  Flux  Fluctuations'); 

text(10,10E-10,'(W/cm^2)^2/Hz') 

text(l  00,1  OE-1 0,’W/cm''2') 

%2.  Plot  qtmeas  x  coherence  and  calculate  qtrms  and  dqt  total  from  measured  data  (no 
modeling): 

Qfmcoh  =  Qfm(l  :65535).*coh; 

qrmsf2  =  sqrt(sum(Qfincoh)*df);  %RMS  of  measured  q  fluctuations  multiplied  by 

coherence  calculated  in  frequency  domain 

Qfmagcoh  =  Qfrnag( 1 :65535).*coh; 

dqtotalcoh2  =  sqrt(2.*sum(Qfmagcoh(l:6555).^2)); 

%  Plot  measured  values,  both  PSD  and  Magnitude  with  coherence  applied: 
figure(9) 
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loglog(fcoh,Qfhicoh,'b’,fcoh,Qfinagcoh,'r') 

xlabeI(Trequency  (kHz)'); 

yIabeI(['PSD  and  Magnitude  of  HF  Fluctuations']); 

title('Figure  D2:  PSD  and  Magnitude  of  Fluctuations  of  Measured  Heat  Flux  Multiplied 
by  Coherence'); 

axis([10e-l  40000  lOe-14  1]); 

legendCMeasured  Heat  Flux  PSD','Magnitude  of  Heat  Flux  Fluctuations')* 
text(l  0, 1 OE-1 0,'(W/cin^2)^2/Hz') 
text(  1 00, 1 OE-1 0,'W/cni^2') 

%3.  Plot  model  predicted  dqt  and  calculate  qtrms  and  dqt  total  from  modeled  data: 

qpturms3  =  sqrt(sum(2.’»qptu(l  :60000).^2));  %Modeled  heat  flux  -  all  frequencies 
dqtotal3  =  sqrt(sum(2.*qptu(l:6555).''2));  %Modeled  heat  flux  out  to  5  kHz 

%  Plot  modeled  values 

figure(lO) 

loglog(f,qptu,'b') 

xlabel('Frequency  (kHz)'); 

ylabel(['Modeled  Heat  Flux  (W/m''2)']); 

title('Figure  D3:  Modeled  Heat  Flux'); 

axis([10e-l  40000  lOe-4  1]); 

%4.  Plot  model  predicted  dqt  and  calculate  qtrms  and  dqt  total  from  modeled  data  with 
coherence: 

qpturms4  =  sqrt(2.*sum(qptcoh(l:60000).^2));  %Modeled  heat  flux  -  all  frequencies 
dqtotal4  =  sqrt(2.*sum(qptcoh(l  :6555).^2));  %Modeled  heat  flux  out  to  5  kHz 

%  Plot  modeled  values  with  coherence 
figure(l  1) 

logIog(fcoh,qptcoh,'b') 
xlabel('Frequency  (kHz)'); 

ylabe^fModeled  Heat  Flux  x  \gamma_u_q^2  (W/cm^2)']); 
titIe('Figure  D4:  Modeled  Heat  Flux  x  \gamma  u_q'''2'); 
axis([10e-l  40000  lOe-14  1]); 

figure(12) 

loglog(f,qptu,'b',fcoh,qptcoh,'r') 
xlabel('Frequency  (kHz)'); 
ylabel(['Modeled  Heat  Flux  (W/cm^2)']); 

titleCFigure  D5:  Comparison  of  Modeled  Heat  Flux  With  and  Without  Coherence')* 
axis([10e-l  40000  lOe-14  1]); 
legend('without  Coherence', 'with  Coherence'); 

%5.  Steady  model  using  u'rms  and  Lx: 
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figure(13) 

loglog(f,Ufmag,’b’,f,Vf2m,'r’) 
xIabel('Frequency  (kHz)'); 
yIabel(['Magnitude  of  u  and  PSD  of  u']); 

title('Figure  D6:  Magnitude  and  PSD  of  Measured  Velocity  Fluctuations'); 
axis([10e-l  40000  lOe-14  1]); 
legendCMagnitude  of  u','PSD  of  u'); 

%6.  Similar  to  steady  model  using  u'/Lx: 
uLxu  =  Ufmag./Lxu; 

uLxurms6  =  sqrt(sum(2.*uLxu(l  :60000).^2));  %A11  frequencies 
qpturms6  =  abs(ke.*(Te  -  T)V(pi.*ae).'^.5).*sqrt(uLxurms6); 
uLxusumb  =  sqrt(sum(2.*uLxu(l:6555).‘'^2));%Out  to  5  kHz 
qptusumb  =  abs(ke.*(Te  -  T)y(pi.*ae).'^0.5).*sqrt(uLxusum6); 

%6.  Plot  u'/Lxu  values: 
figure(14) 
loglog(f,uLxu,'b') 
xlabelCFrequency  (kHz)'); 
ylabel(['u/Lx(l/s)']); 

title('Figure  D7:  Plot  of  umag/Lxu  versus  Frequency'); 
axis([10e-l  40000  lOe-14  1]); 

%6.  Similar  to  steady  model  using  u'/Lx  with  Coherence: 
uLxucoh  =  uLxu(l:65535).*coh; 

uLxurms?  =  sqrt(sum(2.*uLxucoh(l  :60000).^2));  %A11  frequencies 
qpturms?  =  abs(ke.*(Te  -  T)./(pi.*ae).''0.5).*sqrt(uLxurms7); 
uLxusum7  =  sqrt(sum(2.*uLxucoh(l  :6555).^2));  %Out  to  5  kHz 
qptusum7  =  abs(ke.*(Te  -  T)./(pi.*ae).^0.5).*sqrt(uLxusum7); 

%6.  Plot  u'/Lxu  values  with  coherence: 
figure(15) 

loglog(fcoh,uLxucoh,'b') 
xlabelCFrequency  (kHz^); 
ylabel(['u/Lx  x  \gamma_u_q^2  (1/s)']); 

title('Figure  D8:  Plot  of  umag/Lxu  x  \gamma_u_q^2  versus  Frequency'); 
axis([10e-l  40000  lOe-14  Ij); 

%7  Autocorrelation  of  u'  and  calc  of  Lx...,Not  necessary,  already  done  earlier: 

%8  Autocorrelation  ofq’ and  calc  of  Lq: 

%Calculate  autocorrelation  on  each  data  file: 

Rq8=xcorr(qp,'coefF); 

Rq8=Rq8(131070:262139); 

%Plot  autocorrelations: 
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figure(16); 

plot(t.Rq8,V); 

title('Figure  DIO:  Autocorrelation  of  Heat  Flux  Signal*); 
xIabelOtau  (s)');  ylabel('R_q(\tau)’); 
axis([0  0.002  -0.2  1.2]); 
grid  on; 

%Perform  Numerical  Integration  to  Calpulate  Integral  Time  Scale 
%First  zero  crossing  calculated  with  find  command: 
zq=find(abs(Rq8)<0.02&abs(Rq8)>=0); 
Tq=sum(Rq8(l:zq(l)))*dt;  %T  is  integral  time  scale 
Lq=Tq*meanq; 
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